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Abstract 

A review of new developments in theoretical and experimental electronic structure 
investigations of half-metallic ferromagnets (HMF) is presented. Being semicon- 
ductors for one spin projection and metals for another ones, these substances are 
promising magnetic materials for applications in spintronics (i.e., spin-dependent 
electronics). Classification of HMF by the peculiarities of their electronic structure 
and chemical bonding is discussed. Effects of clectron-magnon interaction in HMF 
and their manifestations in magnetic, spectral, thermodynamic, and transport prop- 
erties are considered. Especial attention is paid to appearance of non-quasiparticle 
states in the energy gap, which provide an instructive example of essentially many- 
body features in the electronic structure. State-of-art electronic calculations for 
correlated d-systems is discussed, and results for specific HMF (Heusler alloys, 
zinc-blende structure compounds, Cr02, Fe304) are reviewed. 
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I. INTRODUCTION 

Twenty-five years ago the unusual magnetooptical properties of several Heusler alloys 
motivated the study of their electronic structure. It yielded an unexpected result: some of 
these alloys showed the properties of metals as well as insulators at the same time in the 
same mater ial depending on the s pin direction. This property was baptized half-metallic 



magnetism (Ide Groot et a/.l . Il983bl ). Although it is not exactly clear how many half-metals 
are known to exist at this moment, half-metallic magnetism as phenomenon has been gen- 
erally accepted. Formally the expected 100% spin polarization of the charge carriers in a 
half-metallic ferromagnet (HMF) is a hypothetical situation that can only be approached 
in the limit of vanishing temperature and neglecting spin-orbital interactions. However, at 
low temperatures (as compared with the Curie temperature which exceeds lOOOK for some 
HMF's) and minor spin-orbit interactions a half-metal deviates so markedly from a normal 
material, that the treatment as a special category of materials is justified. The confusion on 
the number of well-established half-metals originates from the fact that there is no "smok- 
ing gun" experiment to prove or dis prove half-metallicity. The most direct measurement is 



spin-resolved positron annihilation (JHanssen and Mijnarendsi . Il986l ). but this is a tedious. 



expensive technique requiring very dedicated equipment. The only proven HMF so far is 
NiMnSb to the precision of t he experiment, which was better than one hundredth of an 



electron ( Hanssen et al. 



19901 ). This number also sets the scale for concerns of temperature- 
induced depolarization and spin-orbit effects, detrimental for half-metallicity. 

The half-metallicity in a specific compound should not be confused with the ability to 
pick up 100% polarized electrons from a HMF. The latter process involves electrons crossing 
a surface or interface into some medium where their degree of spin polarization is analyzed. 
This is clearly not an intrinsic materials property. The richness but also the complications 
of surfaces and interfaces are still not fully appreciated. 

Because of all these experimental complications, it is not surprising that electronic struc- 
ture calculations continue to play an important role in the search for new HMF as well as 
the introduction of new concepts like half-metallic antiferromagnetism. However, electronic 
structure calculations have weaknesses as well. Most of the calculations are based on density 
functional theory in the LDA or GGA approximation. It is well known that these methods 
underestimate the band gap for many semiconductors and insulators, typically by 30%. It 



has been assumed that these problems do not occur in half-metals since their dielectric re- 
sponse is that of a metal. This assumption was disproved recently. A calculation on the 
HMF Lao.rSro.sMnOs employing the GW approximation (that gives a correct description of 
band gaps i n many sern i -cond uctors) leads to a half-metallic band gap 2 eV in excess to the 
DFT value (JKino et all |2003| ). The consequences of this result are possibly dramatic: if it 
were valid in half-metallic magnetism in general, it would imply that many of the materials, 
showing band gaps in DFT based calculations of insufficient size to encompass the Fermi 
energy, are actually in reality bona fide half-metals. Clearly much more work is needed in 
this area. 

The strength of a computational approach is that it does not need samples: even nonex- 
istent materials can be calculated! But in such an endeavour a clear goal should be kept 
in mind. Certainly, computational studies can help in the design of new materials, but the 
challenge is not so much in finding exotic physics in materials that have no chance of ever 
being realized. Such studies can serve didactical purposes, in which case they will be in- 
cluded in this review. However, the main attention will be devoted to materials that either 
exist or are (meta) stable enough to have a fair chance of realization. 

This review will cover half-metals and will not discuss the thriving area of magnetic semi- 
conductors. Some overlap exists, however. The much older field of magnetic semiconductors 
started with se miconductors like the europium monochalcogenides and cadmium-chromium 



chalcogenides (JNagaevL 



19831). Later, the att ention changed to the so-called diluted mag- 



netic semiconductors ( JDelves and Lewisl . ll963h . These are regular (i.e. III-V or II- VI) semi- 
conductors, where magnetism is introduced by partial substitution of the cation by some 
(magnetic) 3(i transition element. The resulting Curie temperatures remained unsatisfac- 
tory, however. The next step in the development was the elimination of the non-magnetic 
transition element altogether. This way HMF's can be realized provided that the remaining 
transition-metal pnictides could be stabilized in the zinc-blende or related structures. The 
review will treat not the (diluted) magnetic semiconductors as such, but some aspects of 
metastable zinc-blende HMF's. 

HMF's form a quite diverse collection of materials with very different chemical and phys- 
ical properties and even the origin of the half-metallicity can be quite distinct. For this 
reason it is required to discuss the origin of the band gap in terms of two ingredients that 
define a solid: the crystal structure and the chemical composition. Two aspects are of im- 



portance in this context. The first one is "strong magnetism" versus "weak magnetism". In 
a strong magnet, the magnetic moment will not increase if the exchange splitting is increased 
hypothetically. Thus the size of the magnetic moment is not determined by the strength 
of the exchange interaction, but is limited instead by the availability of electron states. In 
practice this implies that either the minority spin sub shell(s) responsible for the magnetism 
is (are) empty or the relevant majority channel(s) is (are) completely filled. In the case of 
weak magnetism the magnetic moment is determined by a subtle compromise between the 
energy gain of an increase in magnetic moment (the exchange energy) and the (band) en- 
ergy the increase of the magnetic moment costs. To avoid misunderstanding, we emphasize 
that this defin ition of "weak" and "strong" magnets differs from that used in Moriya's book 



(iMoriyal . Il985l ) and most theoretical works on itinerant-electron magnetism. According to 
Moriya, "strong" magnets are those with well-defined magnetic moments which means, e.g., 
the Curie- Weiss behavior of the wavevector-dependent magnetic susceptibility x(q; T) in the 
whole Brillouin zone. In this sense, all HMF's (e.g., containing Mn ions) are strong magnets. 
However, within this group of materials we may introduce a more fine classification based 
on the sensitivity of magnetic moment with respect to small variation of parameters. 

Any combination of "weak" or "strong" magnetism with majority or minority spin band 
gaps is known today. Thus "weak" magnets with minority spin band gaps are found in 
the Heusler alloys and artificial zinc-blendes, examples of weak magnets with majority spin 
gaps are the double perovskites and magnetite. The colossal magnetoresistance materials 
are examples of strongly magnetic half-metals with minority spin band gaps as is for ex- 
ample, chromium dioxide, while the anionogenic ferromagnets like rubidium sesquioxide are 
examples of a strongly magnetic half-metal with a majority spin band gap. 

An interesting and relatively new development is the work on half-metallic sulphides. 
The HMF state in oxides with the spinel structure is relatively rare. The prime example 
of course is magnetite. However any substitution into the transition-metal sublattices leads 
practically invariable to a Mott insulating state, as magnetite itself possesses below the 
Verwey transition at 120K. On the other hand the electrons in sulphides, are substantially 
less correlated. Hence a wealth of substitutions is possible in order to optimize properties, 
design half-metallic ferro- or antiferromagnets etc. without the risk of losing the metallic 
properties for the second spin direction as well. There is a prize to be paid however: since the 
cation-cation distances are larger in sulphides, the Curie and Neel temperatures are lower 
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in comparison with the oxides. Nevertheless, the work on half-metalhc sulphides deserves 
strong support. 

In all the metallic ferromagnets, interaction between conduction electrons and spin fluctu- 
ation is of crucial importance for physical properties. In particular, the scattering of charge 
carriers by magnetic excitations determines transport properties of itinerant magnets (tem- 
perature dependences of resistivity, magnetoresistivity, thermoelectric power, anomalous 
Hall etc.). From this point of view, HMF, as well as ferromagnetic semiconductors, differ 
from "normal" metallic ferromagnets by the absence of spin-flip (one-magnon) scattering 
processes. This difference is also important for magnetic excitations since there no Stoner 



damping, and spin waves are wel 



insulators (jAuslender and Irkhin 



defined in the whole Brillouin zone, as well as in magnetic 



1984al : 



Irkhin and Katsnelsonl . 



19941 ). 



Electron-magnon interaction modifies also considerably electron energy spectrum in 
HMF. These effects take place both in usual ferromagnets and in HMF. However, peculiar 
band structure of HMF (the energy gap for one spin projection) results in important conse- 
quences. In generic itinerant ferromagnets the states near the Fermi level are quasiparticles 
for both spin projections. On the contrary, in HMF and important role belongs to inco- 
herent (non-quasipartic l e, NQ P) states which occur near the Fermi le yel in the energy gap 



Irkhin and Katsnelson. 



1973 



1994) 



Irkhin and Katsnelson 



. The appearance of the NQP states in (JEdwards and Hertzl . 



19831 ) is one of the most interesting correlation effects typi- 



cal for HMF. The origin of these states is connected with "spin-polaron" processes: the 
spin-down low-energy electron excitations, which are forbidden for HMF in the one-particle 
picture, turn out to be possible as a superposition of spin-up electron excitations and virtual 
magnons. The density of the NQP states vanishes at the Fermi level but increases drastically 
at the energy scale of the order of a characteristic magnon frequency ID. These states are im- 



portant for spin-polarized elect ron spectroscopy (jlrkhin and Katsnelsonl . 



tions (Andreev reflection) ( jTkachov et al. 



2001 



2006|), NMR 



(jlrkhin and Katsnelsonl . l2005al). and subgap tran sport in ferromagnet-superconductor junc 



200ll ). Recently, the dei isity of NQP states ha s 



been calculated from first principle s for a prototype HMF , NiMnSb ( IChioncel et a/.l . l2003al ). 



as well as for other Heusler a lloys (IChioncel et al 



(IChioncel et al. 



2005 



2006al ) , zinc -blend structure compounds 



2006bl ) and Cr02 (IChioncel et all 120071 ). Fig [T| shows in a pictorial 



way the NQP contribution to the density of states. 

Therefore, HMF are very interesting conceptually as a class of materials which may be 



convenient to treat many-body solid state physics that is essentially beyond band theory. It 
is accepted that usually many-body effects lead only to renormalization of the quasiparticle 
parameters in the sense of Landau's Fermi liquid ( FL) theo r y, the electronic liquid being 
qualitatively similar to the electron gas (see, e.g., ( iNoziered . Il964j )). On the other hand, 
NQP states in HMF are not described by the FL theory. As an example of highly unusual 
properties of the NQP stat es, we mention that they can contribute to the T-linea r term in 



1990 



Irkhin et al. 



1989 



19M), despite 



the electron heat capacity fllrkhin and Katsnelsonl . 

their density at Ep is zero at temperature T = K. Some developments concerning physical 

effects of NQP states in HMF are considered in the present review. 



II. CLASSES OF HALF-METALLIC FERROMAGNETS 



A. Heusler alloys and zinc-blende structure compounds 

In this chapter we will treat HMF with the Heusler Clt and L2i structures. Although 
being not Heusler alloys in the strict sense, artificial half-metals in the zinc-blende structure 
will be also discussed because of their close relation with the Heusler CI;, ones. Zinc blende 
has a face centered cubic {fee) Bravais lattice with a basis of (0,0,0) and (1/4,1/4,1/4), 
both species coordinating each other tetrahedrally. The Heusler 01^ structure consists of 
the zinc-blende structure with an additional occupation of the (1/2, 1/2, 1/2) site. Atoms 
at the latter position, as well as those in the origin, are tetrahedrally coordinated by the 
third constituent, which itself has a cube coordination consisting of two tetrahedra. The 
Heusler L2i structure is obtained by an additional occupation of the (3/4, 3/4, 3/4) by the 
same element already present on the (1/4, 1/4, 1/4) site. This results in occurrence of an 
inversion center which is not present in the zinc-blende and Heusler CI;, structures. This 
difference has important consequences for the half-metallic band gaps. Electronic structure 
of the Heusler allo ys has been reviewed recently (frorn a bit different positions) by Galanakis 
and Mavropoulos (JGalanakis and Mavropoulod . 



20071 ). 



1. Heusler Clf, alloys 

Interest in fast, non- volatile mass storage memory sparked many activities in the area of 
magnetooptics in general and the magnetooptic Kerr effect specifically in the beginning of 
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the 80-s of last century. All existing magnetic solids were i nvestigated leading to a record 



MOKE rotation of 1.27° for PtMnSb (Ivan Engen et al. 



19831 ) . The origin of these properties 



remained an unsolved problem. This formed the motivation for the study of the electronic 
structure of the isoelectronic Heusler CI;, compounds NiMnSb, PdMnSb and PtMnSb and 
the subsequent discovery of half-metallic magnetism. Interestingly enough, there seems to 



be still no consensus on t he origin o 
and intuitive explanation (jde Groot , 



the magnetooptical properties. The original simple 



199ll ) was complementary to the production of spin- 



polarized electrons by optical excitation in III-V semiconductors. In that case, the top of the 
valence band is split by the the spin-orbit coupling, and the photoexcitation of electrons from 
the very top of the band by circularly polarized light leads to 50% spin polarization. Vice 
versa, excitations from a lOOvalenceband is possible for only one of the two components of 
circular light, as in the case of PtMnSb, this should result in a strong difference of refraction 
and absorption for two opposite polarizations. In PtMnSb this difference is maximal for the 
visible light, and for NiMnSb the maximum of off-diagonal optical conductivity is shifted 
to the ultraviolet region. The main contribution to th is shift comes from scalar relativistic 
interactions in the final state (jWijngaard et alV Il989l ). which are much weaker for Ni than 
for Pt due to difference in nuclear charges. Fur ther the magnetoopt ical properties of the 



Heusler alloys were calculated by Antonov et al. (jAntonov et al. 



19971 ) in a good agreement 



with experim ental data, but 



;he p hysical explanation was not given in this paper. Recently, 



Chadov et al ( IChadov et all 



20061 ) demonstrated that an agreement between calculated and 



experimental values of the Kerr rotation and ellipticity in NiMnSb can be further improved 
by taking into account correlation effects within so-called LDA-I-DMFT approach (see below 

SectUVZl). 

Since NiMnSb is the most studied HMF (at least within the Heusler Gift's) we will 
concentrate on it here. The origin of its half-metallic properties has an analogy with the 
electronic structure of III-V zinc-blende semiconductors. Given the magnetic moment of 4/^b 
manganese is trivalent for the minority spin direction, and antimony is pentavalent. The 
Heusler Gl^ structure is the zinc-blende one with an additional site (1/2, 1/2, 1/2) being 
occupied. The role of nickel is to supply both Mn and Sb with the essential tetrahedral 
coordination and to stabilize MnSb in the cubic structure (MnSb in the zinc-blende structure 
is half-metallic, but not stable). Thus a proper sit e occupancy is e s sentia l : nickel has to 



occupy the double tetrahedrally coordinated site (jHelmholdt et al 



1984 



Orgassa et al. 
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19991 ). The similarity in chemical bonding between NiMnSb and zinc-blende semiconductors 
also explains why it is a "weak" magnet, in the sense discussed in the Introduction: the 
presence of occupied manganese minority d states is essential for the band gap. These states 
play the same role as the metal p states in zinc-blende semiconductors, a situation which 
is possible only because of the absence of inversion symmetry . The similarity of che mical 
bonding in Heusler's and zinc-blende in the original paper (Ide Groot et all Il983bl ) was 
illustrated by "removing the nickel d states from the Hamiltonian" . This phrasing has lead 
to considerable confusion. Actually, the coupling of manganese states and antimony states 
through non-diagonal matrix elements of nickel d states was maintained in this calculation. 
Se veral explanations of th e band gap have been given in terms of a Ni-Mn interaction 



only (JGalanakis et al 



2002al ). While this interaction is certainly present, it is not sufficient 



to explain the band gap in NiMnSb. These analyses are based on calculations of NiMnSb 
excluding the antimony, but keeping the volume fixed. This is a highly inf lated situation with 



a volume more than twice the equilibrium one (JEgorushkin et al. 



19831 ). Under expansion 



bandwidths in metals decrease, leading eventually to a Mott insulating state. But even 
before this transition a band gap appear, simply due to the infiation itself. This is not 



a hypothetical scenario: a so lid as simp 



ferromagnet under expansion (iMin et al. 



e as elemental lithium becomes a half-metallic 



19861 ). yet there is no evidence for half- metallic 



magnetism under equilibrium conditions for this element. Also, it is not clear from these 
considerations why NiMnSb is half-metallic only in the case of tetrahedrally coordinated 



manganese . Probab 
by in Ref. ( JKiiblerl . 



y the chemical bonding in relation to the band gap is best summarized 



20001 ): a nickel induced Mn-Sb covalent interaction. 



Surfaces of NiMnSb d o not show the 100% spin polar i zation as determined by positron 



annihilation for the bulk (JBona et al. 



1985 : 



Soulen et al 



tendency of showing surface segregation of manganese ( 



1999). Part of the reason is their 



Ristoiu et al. 



2000) • Also, surfaces of 



NiMnSb are quite reactive and are easily oxidized. But even without contaminations none o 



2003; 



de Wijs and de Groot 



the s urfaces of NiMnSb are genuinely half-metallic (iGalanakis 
20011 ). This is just another example of the sensitivity of the half-metallic properties in 
NiMnSb on the correct crystal structure. But this does not necessarily imply that interfaces 
of NiMnSb with, for example, semiconductors cannot be completely spin-polarized. For 
example, it was shown that at the 1116 interface of NiMnSb with CdS or InP the HMF 
properties are completely conserved, if the semiconductors are anion terminated at the 
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interface (jde Wijs and de Grootl . l200ll ). This anion-antimony bond may look exotic, but 
such a coordination is quite common in minerals like costobite and paracostobite (minerals 
are stable on a geological timescale) No experimental verification is at hand at this moment, 
partially because experimentalists tend to prefer the easier 100 surfaces in spite of the fact 
that calculations show that no half-metallic properties are possible here^ 



Several photoemission measurements have been reported on NiMnSb (jCorrea et allUOOw 
as well as the closely related PtMnSb (JKisker et all 119871 ). We concentrate on the latter 
one, because it is the first angular-resolved measurement using single crystalline samples. A 
very good agreement with the calculated band structure was obtained. This is a remark- 
able result. In calculations based on density functional theory eigenvalues depend on the 
occupation. These occupations deviate from the ground state in a photoemission experi- 
ment. To very good precision the dependence of the eigenvalues on the occupation numbers 
is given by the Hubbard parameter U. The consequence is that the effective U value in 
alloys like NiMnSb and PtMnSb is much smaller then, e.g., in Ni metal where photoe- 
mission experiments indic ate the satellite structure related to the so-called Hubbard bands 
(JLichtenstein et aliuOQu ). 



The transpor t properties of NiMnS b were studied extensively (JBorca et al 



Hordequin et al. 



2000 



Otto et al. 



2001 



19891 ) (a theoretical discussion of transport properties 



in HMF is given in Sect. IIII.I|) . At low temperatures, the temperature dependence of the 
resisti vity follows a T^ law. However, the T^ law at low temperatures is absent in thin 
films ( JMoodera and Mootod . Il994l ). At around 90K a transition takes place, beyond which 
the temperature dependence is T^-^^ , The nature of this phase transition is unknown. One 
possibility is the effect of thermal excitations provided the Fermi energy is positioned close 
to the top of the valence band or the bottom of the conduction band. For example, in 
the latter case thermal excitations are possible from the metallic majority spin direction to 
the empty states in the conduction band of the minority spin direction. Such excitation 
reduces the magnetic moment, which in its turn will reduce the exchange splitting resulting 
in an even further reduction of the energy difference between Fermi energy and bottom of 
the conduction band. This positive feed-back will lead to a collapse of the half-metallic 
properties at a certain temperature. An analogous situation exists for Ep close to the va- 
lence band maximum. Numerical simulations indicate that this scenario is highly unlikely 
in the case of NiMnSb because of its unusual low density of states at the Fermi energy. 
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Another explanation has been put forward b ased on the crossing of a magnon and a phonon 



branch at an energy corresponding to 80K (JHordequin et al 



1997a 



b|). It is unclear how 



this phonon-magnon interaction influences the electronic properties of NiMnSb. 

Local magnetic moments were studied experimentally as a function of temperature with 
polarized neutron scattering. The manganese moment decreases slightly with temperature 
from 3.79 ^b at 15K to 3.55 fiB at 260K, while the nickel moment remains constant at 
0.19 fiB in the same temperature range. On the other hand, magnetic circular dichroism 
shows a reduction of both the manganese and nickel moments around 80K. Borca et al. 
(JBorca et al\ . l200lh conclude that at the phase tra nsition the couplin g of the manganese 
and nickel moments is lost. A computational study (iLezaic et all |2006| ) shows vanishing of 
the moment of nickel at the transi tion temperature. None of these anomalies are reflected 



in the spontaneous magnetization (jOtto et al 



19891) of bulk NiMnSb. 



Two Heusler 01^ alloys exist, isoelectronic with NiMnSb: PdMnSb and PtMnSb. Their 
electronic structures are very similar, but we will discuss the differences. The calculated DFT 
band structure for PdMnSb is not half-metallic. The minority spin direction does show a 
band gap very similar as NiMnSb but the Fermi energy intersects the top of the valence 
band. Reliable calculations (e.g., based on the GW approximation) are needed here to 
settle the issue whether PdMnSb is half-metallic or not. PtMnSb is very similar to NiMnSb, 
the largest differences being in the empty states just above the Fermi energy. The direct 
band gap (at the gamma point) is between the triplet top of the valence band (neglecting 
spin-orbit interactions) and a total symmetrical singlet state in PtMnSb. This singlet state 
is positioned at much higher energy in NiMnSb. These differences have been attributed t o 
the much stronger mass- velocity and Darwin terms in platinum ( Wijngaard et all Il989l ). 



Platinum does not carry a magnetic moment in PtMnSb. Consequently, no 90K anomaly 
like in NiMnSb is to be expected and none has been reported so far. 

Let us consider whether half-metals in the Heusler CU structure exist when substitut- 
ing other than iso-electronic elements for Ni in NiMnSb. Since NiMnSb is a weak magnet, 
substitutions with elements reducing the total magnetic moment are possible only while 
maintaining the half- metallic properties. Thus cobalt, ir on maiiganese and chromium will 
be considered. The case of Co was studied by Kiibler (JKubleiJ, Il984l ). The half- metallic 



properties are conse rved, consequeii t 



lations on FeMnSb (Ide Groot et al. 



y the magnetic n aoment is reduced to 3 /iR. Calcu- 



19861 ). MnMnSb flWijngaard et all Il992l ) as well as 
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CrMnSb (jde Grootl . Il99ll ) all show the preservation of the half-metallic properties. In the 
case of FeMnSb this implies a reduction of the total magnetic moment per formula unit to 
2fiB, which is an unusually small moment to be shared by iron and manganese. The way 
out is that FeMnSb orders antzferromagnetically. Thus the 2fiB total magnetic moment 
corresponds to the difference in moments of iron and manganese rather than to their sum 
implied by ferromagnetic ordering. This way to preserve a band gap (energetically favorable 
from a chemical-bonding point of view) together with the maintenance of sizable magnetic 
moments (favorable for the exchange energy) determine the magnetic ordering here. Both 
these effects are usually larger than the exchange-coupling energies. The antiferromagnetic 
ordering is maintained in MnMnSb with a total moment of IfiB- In the case of CrMnSb 
the antiferromagnetic coupling leads to a half- metallic solution with a zero net moment. 
This is a really exotic state of matter. It is genuinely half-metallic, im plying 100% spi n 



polarization of the conduction-electrons, yet it lacks a net magnetization (jde Grootl . 1 19911 ). 
The stability of such a solution depends sensitively on the balance between the energy gain 
of the band gap and the energy gain due to the existence of magnetic moments: if the first 
one dominates a nonmagnetic semiconducting solution will be more stable (remember that 
both spin directions are isoelectronic here). 

In reality, the situation is more complex. CoMnSb does exist, but i t crystallizes in 



a tet ragonal superstructure with Co partially occupying the empty sites (ISenateur et al. . 



19721 ). The magnetic moments deviate from the ones expected for a half-metallic solution. 



FeMnSb does not exist, but part of the nickel can be substituted by iron. Up to 10% 
the Heusler CI;, structure is maintained, from 75% to 95% a stru cture comparab l e wit h 



CoMnSb is stable and between 10% and 95% both phases coexist (jde Groot et all Il986l ). 
MnMnSb exists, orders antiferromagnetically and has a net moment of I^b- It does not 
crystallize in the Heusler Clf, structure and consequently is not half-metallic. CrMnSb exists, 
is antiferromagnetic at low temperatures, and shows a transition at room temperature to a 
ferromagnetic phase. 

A different substitution is the replacement of Mn by another transition metal. An in- 
teresting substitution is a rare earth element R. Because of the analogy of half-metals with 
Clfo structure and III-V semiconductors one expects NiRSb compounds to be nonmagnetic 
semiconductors. 

Several of these compounds do exist in the Cl^ structure, examples are Sc, Y and heavy 
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rare-earth elements from the second half of the lanthanide series. All of them are semi- 



conductors indeed (jPierre and Karla 



2000 



Pierre et al. 



1999f ). Doping of NiMnSb by rare 



earth elements has been suggested as a way to improve the finite-temperature spin polariza- 
tion in NiMnSb. These substitutions do not influence the electronic band structure much, 
(see also Sect. IV.A.2p . the band gap for the minority spin direction remaining completely 
intact. However, the random substitution of nonmagnetic (Y, Sc) or v ery different (Ho-Lu) 



magnetic elements for manganese will modify the magnon spectrum (jAttema et al 



Chioncel 



2004 



20041 ). This could be beneficial to increase spin polarization in some temperature 



range. 



2. Half-metals with zinc-blende structure 



The Curie temperatures of diluted magnetic semiconductors remain somewhat disap- 
pointing. A solution is to replace all the main group metals by transition metals. But there 
is a heavy prize to be paid: These systems can only be prepared as metastable states - 
if at all - on a suitable chosen substrate. An alternative to come to the same conclusion 
is to consider Heusler Clf,'s with larger band gaps. This is most easily accomplished by 
replacement of the antimony by arsenic or phosphorous. No stable Heusler 01^ alloys exist 
with these lighter pnictides, however. An alternative is to try to grow them as metastable 
systems on a suitable chosen substrate. This makes the nickel superfiuous: since it fails in 
the case of lighter pnictides to play the role it does so well in the NiMnSb. The bottleneck 
in this quest is not so much in predicting systems that are good half-metals, but to design 
combinations of half-metals and substrates that are meta-stable enough to have a chance of 
being realized e xperimentally . 



Shirai et al. (IShirai et ali 119981 ) were the first to relate concentrated magnetic semicon- 
ductor with half-metallic magnets in their study of MnAs in the zinc-blende structure. The 
experimental realization showed an increase of th e Curie temperat ure indeed: 400K was 
reported on for CrAs grown on GaAs. Xie et al. (jXie et all l2003bl ) calculated the stabil- 
ity of all 3d transition-metal chalcogenides in the zinc-blende with respect to the ground 
state structure. Chromium telluride and selenide, as well as vanadium telluride, are good 
half-metals which are stable towards a tetrahedral and rhombohedral distortions. Zhao and 



Zunger (IZhao and Zungerl . l2005l ) consider the stability of an epitaxial layer as a function of 
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the lattice parameter of the substrate allowing for relaxation in the growth direction. The 
result is that while the bulk zinc-blende phase is always unstable with respect to the (equi- 
librium) NiAs structure, there exist lattice constants where the epitaxial zinc-blende phase 
is more stable as compared with the epitaxial nickel arsenide structure. This is realized 
(computationally) for half-metallic CrSe. 

An alternative to the concentrated III-V magnetic semiconductors is given by delta doped 
IIl-V semiconductors. Here the magnetic properties are not introduced by a more or less 
homogeneous replacement of main group metals by magnetic transition metals. Instead, a 
very thin transition- meta l layer is sandwiched between undoped III-V s emiconductor mate - 



rial ( Nazmul et al. 



20021 ). A clear increase in Curie temperature results ( Chiba e^ a/.l.l2003r) . 



This is not unrelated to the interface-half-metallicty introduced before (Ide Groot 



199l|). 



3. Heusler L2i alloys 

The crystal structure of the Heusler L2i alloys is closely related with that of the Clf, alloys. 
In the L2i structure the (1/2, 1/2, 1/2) position, empty in the 01^ structure, is occupied by 
the same element that occupies the (0, 0, 0) position. The similarity in structure suggests 
a similarity in interactions and physical properties, but on the contrary, the interactions 
and the physical properties of the two classes are actually quite distinct. The introduction 
of the fourth atom in the unit cell introduces an inversion centre in the crystal structure. 
The bandgap in the 01^ compounds resulted from an interaction very similar to that in 
III-V semiconductors, where the manganese t2g d electrons play the role of the p electrons 
in the III-V semiconductor. This is no longer possible in the presence of an inversion 
centre. Consequently bandwidths are reduced and usually Van Hove singularities occur in 
the vicinity of the Fermi energy. The smaller bandwidth leads to several (pseudo) gaps. 
Correlation effects are expected to become better observable here. 

Another difference is the occurrence of defects. Experimentally it was noted that "The 
strong effect of cold work on Heusler alloys (L2i structure) contras ts with almost unobserv- 



able effects in the CI?, structure alloy NiMnSb" (iSchaf et ami983l ). But also here there are 
indications that defects that destroy the bandgap are energetically less favourable. 

Experimental work goes back to Heuler in the beginning of the last century. The motiva- 
tion of his work was the possibility of preparing magnetic alloys out of non- magnetic elements 
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(JHeuslerl . Il903l ) (A material was only considered magnetic in that period if it possessed a 
spontaneous net magnetisation). More recently the landmark work of Ziebeck and Webster 
on neutron-diffractions investigations (jZiebeck and Websterl . 119741 ) deserves mentioning as 
well as the NMR work in the Orsay group of Campbell. 



The f irst b andstructure calculations were by Ishida an d coworkers (llshida et all 1 19821 . 



1976a| ]bl. 119801 ). as well as Kuebler, Williams and Sommers (JKiibler et al. 



3). 



19831). The latter 



paper contains a clue to half-metallic properties in the L2i compounds: the authors remark 
that "The minority state densities at the Fermi energy for ferromagnetic Co2MnAl and 
Co2MnSn nearly vanish. This should lead to peculiar transport properties in these two 
Heusler alloys". 

Calculations that explicitly addressed the question o f half- metallic properties i n the full 



heuslers appeared not earlier than in 1995 (JFujii et al 



1995 : 



Ishida et al 



19951). A sys- 



tematic study of the electronic structur e of Heusler L2i comp ounds was undertaken by 
Galanakis, Dederichs and Papanikolaou (JGalanakis et all l2002bl ). This paper also reviews 
the work on half-metallic magnetism in full heusler compounds till 2002. For this reason we 
refer to it for details and concentrate on subsequent developments here. 

The heusler L2i compounds take a unique position in the spectrum of halfmetals because 
of their Curie temperatures. High Curie temperatures are important in the application of 
halfmetals at finite temperature, since many of the depolarisation mechanisms scale with the 
reduced temperature T /T^. Curie temperatures approach lOOOK: Co2MnSn shows a Curie 
temperature of 829K, the Germanium analogue 9 05K while Co ^ MnSi was a record holder 
for some time with a Curie temperature of 985K (jBrown et all l2000l ). A further increase 
was realized in Co2FeSi. Experi mentally it shows an i nteger magnetic moment of 6 /ig and 



a Curie temperature of HOOK (jWurmehl et al 



20051 ) . This result was not reproduced in 



calculations employing the LDA approximation. The magnetic moment of 6 Bohr-magnetons 
could only be reproduced by the application of U in excess of 7.5 electron-volt. This is an 
unusual high number and alternative explanations should also be considered. The question 
of lattice defects has been studied. On the basis of neutron-diffraction, Co-Fe disorder could 
be excluded, but no data are available for the degree of Fe-Si interchange. A calculation of 
the magnetic saturation moment as function of the iron-silicon disorder seems a logical next 
step in the understanding of this fascinating compound. 

Whereas the investigations of the bulk electronic structures of full heuslers has advanced 
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comparable as with the half-heuslers, the situation with respect of the preservation of 
halfmetaUic properties at surfaces and interfaces clearly still lacks behind. Two impor- 
tant results were obtained recently. One result is the preservation of halfmetaUic properties 
of an Co2MnSi (001) surface provided it is purel y manganese terrainated . This is the only 
surface of this half-metal showing this property (JHashemifar et all |2005| ). 

No genuine halfmetaUic interfaces between full-heuslers and semiconductors are reported 
yet. But the results for Co2CrAl/G a.As look proraising . For an (110) interface a spin- 
polarization of ~ 90% was obtained (JNagao et all |2004J ). Although this is clearly not a 
genuine half-metallic interface yet it should provide a good basis for analyses why half- 
metallic behaviour is lost at in interface in analogy with the successful work for the Clf, 
case. 

An interesting development in half-metallic-magnetism is in electron-deficient full-Heusler 
alloys. Reduction of the number of valence electrons to 24 per formula unit leads to either 
a non-magnetic semiconductor or a halfmetaUic anti-ferromagnet. But remarkable enough, 
the reduction or the number of valence-electrons can be continued here, re-entering a range 
of half-metals but with a bandgap for the majority spin-direction now. This is best ex- 
emplified for the case of Mn2VAl. It is a halfmeta Uic ferrimagnet of calcula ted with the 
generalized gradient exchange-correlation potential (jWeht and Pickettl . Il999l ). Halfmetals 
with a bandgap for the majority spin-direction hardly occur. The search for new candidates 
should strongly be supported. 



B. Strongly magnetic half-metals w^ith minority spin gap 



1. Chromium dioxide 



Chromium dioxide is the only metallic oxide of chromium. It orders ferromagnetically 



with a Curie tempera ture of about 390 K. Its ha l 



structure calculations (JMatar et al . 



1992 



Schwarz 



metallic state was discovered by band 



19861 ). The origin of the half-metallicity 



is straightforward: in an ionic picture the chromium is in the form of a Cr^+ ion. The two 
remaining d electrons occupy the majority d states. The crystal field splitting is that of 
a (slightly) deformed octahedron. The valence band for the majority-spin direction is 2/3 
filled, hence the metallic properties. The minority-spin d states are at a significant higher 
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energy due to the exchange sphtting. For this reason the Fermi level falls in a band gap 
between the (filled) oxygen 2p states and the (empty) chromium d states. Thus the HMF 
properties of chromium dioxide are basically a property of chromium and its valence and, as 
long as the crystal-field splitting is not changed too drastically, the half-metallic properties 
are conserved. This implies that the infiuence of impurities should not be dramatic and 
a number of surfaces retain the half-metallicity of the bulk. As a matter of fact, all the 
surfaces of low index are half- metallic with a possible exception of one of the (101) surfaces 



(Attema et al. 



2006 



van Leuken and de Groot 



not confirm these expectat ions f|Kamper et al. 



1995 



. Although initial measurements did 



19871 ). they were confirmed late r by the 



experiments like tunneling (JBratkovskyl . 119971 ). Andreev refiection (jji et all 1200 ll ) on well- 
characterized surfaces. Recently, the fiow of a triplet-spin s uper current has bee n realized in 
Cr02 sandwiched employing two superconducting contacts (JKeizer et a/.l . l2006l ). 

As mentioned before, an interesting question is the origin of the metallic ferromagnetism 
in Cr O^. This was explained in terms of the double exchange (Zener) model by Korotin 



et al. ( jKorotin et al. 



1998 : 



Schlottmann 



20031 ). The octahedral coordination in the rutile 
structure is slightly distorted. This leads to splitting of the degenerate t2g state into a more 
localized d^y state and more delocalized d^z and dy^ states (or linear combinations of these). 
The localized filled dxy state plays the same role as the filled t2g majority spin state in 
the Zener double-exchange model, while the partially occupied d^y ± dyz majority states in 
Cr02 the role of the p artially occupied e ^ states. The transport properties of Cr02 were 
investigated in detail (jWatts et all 12000) and interpreted in terms of a two-band model, 
very much in line with the double exchange model for Cr02. 

The importance of explicit electron-electron interactions in C rO^ remains a subj ect of 
active research. On one hand, Mazin, Singh and Ambrosch-Draxl (JMazin et all Il999l ) com- 
pared LSDA calculations with experimental optical conductivities and found no indications 



for strong correl ations relatec 



Miiller-Hartman (ICraco et al. 



exotic phenome i ia. O n the other hand, Craco, Laad and 



2003 



Laad et al 



200ll ) considered photoemission results and 



conductivity (both DC and optical) and concluded the importance of dynamical correlation 
effects. The ferrom agnetic correlated stat e was investigated also in a combined local and 
non-local approach (jChioncel et all 120071 ) which demonstrates that the d^ y orbital is not 



comp 



1998 



etely filled and loca 



Schlottmann 



2003; 



ized as described by L DA-I-U or model calculations (JKorotin et al. 



Toropova et al 



20051 ). More recently, Toropova, Kotliar, Savrasov 
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and Oudovenko (JToropova et ali l2005l ) concluded that the low-temperature experimental 
data are best fitted without taking into account the Hubbard U corrections. Chromium 
dioxide will clearly remain an area of active research. 



2. The colossal magnetoresistance materials 

The interest in ternary oxides of manganese with d i - or t rivalent main group metals goes 
back to Van Santen and Jonker (jJonker and Santenl . Il950l ). The occurrence of ferromag- 
netism in transit i on m etal oxides, being considered unusual at that time, was explained 
by Zener ( IZenerl . 



1951 



J) by the introduction of "double exchange" mechanism. In 70th 
and 80th these systems were investigated theoretically in connection with the problem of 
phas e separation and "ferron" ( r nagne t ic polaron) for mation in ferromagnetic semiconduc- 



tors (jAuslender and Katsnelson 



1982 



Nagaevl . Il983l ). The interest in spintronics fifteen 



years ago revived the interest in the ternary manganese perovskites, generally referred to as 
colossal magnetoresistance materials. A wealth of interesting physics is combined in a single 
phase diagram of, for example, Lai.^Sr^MnOa. From a "traditional" antiferromagnetic in- 
sulator for X = 1, the reduction of x results in a ferromagnetic metallic state, while finally at 
X = a Mott insulating antiferromagnet is found. Some of the transitions are accompanied 
with charge and or orbital ordering. Finite temperatures and applied magnetic fields com- 
plicate the phase diagram substantially. The ferromagnetic metallic phase for intermediate 
values of x is presumably half- metallic (jPickett and Singhl . Il996l ). We will concentrate on 
this phase here and refe r to other reviews for a more complete overvi ew of the manganites 



(JDagotto 



2003 



Nagaev 



2001 



Salamon and Jaime 



2001 



Ziese 



2OO2I). 



Once the occurrence of a ferromagnetic magnetic ordering is explained, the occurrence 
of half-metallic magnetism is rather straightforward. Manganese possesses around 3.5 d 
electrons in the metallic high-spin state; its rather localized majority spin t2g state is filled, 
the majority, much more dispersive, Cg state is partially occupied and the minority d states 
are positioned at higher energy, thus being empty. Hence a rather large band gap exists for 
the minority spin at the Fermi energy and the manganites are strong magnets. Correlation 
effects are expected to be much stronger here. Notice that no reference has been made 
to the actual crystal structure: subtleties like in the Heusler structure are absent here. 
The half-metallic properties are basically a property of the valence of the manganese alone. 
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Surface sensitivity of the HMF properties is not to be expected as long as the valence 
of the manganese is main tained: This is easily accomplished in the layered manganites 



(Ide Boer and Groot 



19991). 



Experimental verification of the half-metallic properties has not been without debate. 
The origin of the controversy is that the calculated position of the Fermi energy in the 
energy gap is invariably very close to the bottom of the conduction band. The experimental 



conformation of the HMF behavior by photoemission (jPark et al 



19981 ). was contested on 



the basis of Andreev ref l ection measurements, that did show minority-spin d states at the 



Fermi energy (JNadgorny 



2007; 



Nadgorny et al. 



casted doubt on the half-metallic properties (jJo et al. 



20011). Also, tu n nelling exper i ments initially 



2000 



IViret et al. 



19971 ). Mazin 



subsequently introduced the concept of transport half-metal: the Fermi energy may straddle 
the bottom of the minority spin t2g band, but since these states are loc alized this does not 



influence 



: he ha lf-metallic properties, as far as transport is concerned ( iMazin et al. . 



Nadgorny 



1999 



20071 ). Recent magnetotr ansport measuremen ts on better samples support the 



Kino et al. (JKino et al. 



HMF pictur e of the CMR ma terials (jBowen et al\ . l2003l ). The recent GW calculations by 



20031 ) shed a different light on this matter. In these calculations 



the half-metallic band gap is increased by as much as 2 eV with respect to the DFT value. 
This implies that the minority spin d band is not even close to the Fermi energy and the 
CMR materials should be considered as genuine, real HMF's. 

C. Weakly magnetic half-metals with majority spin gap 



1. The double perovskites 

The double perovskites have a unit cell twice the size of the regular perovskite structure. 
The two transition-metal sites are occupied by different elements. Double perovskites are 
interesting f or two reasons . Half-metallic antiferromagnetism has been predicted to occur for 



LaaVMnOg flPickett 



19981 ) (we will return to this question later in Sect JII.ET2l) . The second 



reason the double perovskites are important is that the high Curie temperatures can be 
obtained in them as compared with the regular perovskites. SraFeMoOe w as the first exam 



pie to be studied in this respect by means of band structure calculations ( jKobayashi et al. 



19981 ). The density of states is shown in figure [3l In the majority spin direction the valence 
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band consists of filled oxygen 2s and 2p states, as well as a completely filled Fe 3d band, 
showing the usual crystal field splitting. A band gap separates the conduction band, which 
is primarily formed by molybdenum d states. The minority spin direction shows an occu- 
pied oxygen derived valence band and a hybridized d band of mixed iron and molybdenum 
character. It intersects the Fermi energy. 



The Curie temperature is in the range of 



was found for Sr2FeRe06 (JKobayashi et al. 



10 to 450K. More recently, a similar behavior 



19991 ). Optical measurement s did show exci 



tations across the half-metallic band gap of 0.5 eV (JTomioka et all 12000 ) . A substantial 



higher Curie temperature is found in Sr2CrRe05, Tc = 635K (JKato et al. 



2002), but band 



structure calculations sho w that the band gap for th e majority spin direction is closed by 



the spin-orbit interaction (jVaitheeswaran et al. 



I2OO5I) 



2. Magnetite 

Magnetite Fe304 is one of most wide-spread natural iron compounds and the most ancient 
magnetic material known to humanity. Surprisingly, we still have no complete explanation 
of its magnetic, electronic and even structural properties, many issues about this substance 
remaining controversial. At room temperature magnetite has inverted cubic spinel structure 
with tetrahedral A-sites occupied by Fe'^"^ ions, whereas octahedral B-sites are randomly 
occupied by Fe^"^ and Fe^+ ions with equal concentrations. FeoO^jsa ferrimagnet with a 
high Curie temperature Tq — 860 K. As discovered by Verwey (IVerweyl . ll939l ). at Ty ~ 120 
K magnetite undergoes a structural distortion and metal-insulator transition. Usually the 
Verwey transition is tr eated as a charge ordering of Fe^"*" and Fe^^ states in octahedral 



sites (for a review, see (JMott 



1974 



19801) )• The nature of the Verwey transition and low- 



temperature phase of Fe304 is a subject of numerous investigations which are beyond our 



topic, see, e.g., recent reviews fiGarsia and Sabia; 



iad. 



2004 



Walz 



2OO2I ). As demonstrated by 



the band-structure calculation (lYanase and Siratoril . 119841 ). magnetite in the cubic spinel 
structure is a rather rare example of HMF with majority-spin gap. This means a satu- 
rated state of itinerant 3 d electrons propagating over B-sites, the magnetic moment being 
close to 4/iB per formula unit. Recentl y this pictur e was questioned by the x-ray mag- 



netic circular dichroism (XMCD) data (JHuang et all |200J) which were interpreted as an 



evidence of large orbital contribution to the magnetization and non-saturated spin state. 
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However, later XMCD experiments (JGoering et all |2006| ) confirm the purely spin saturated 
magnetic stat e. Direct measureme nts of spin polarization by spin-polarized photoemission 
spectroscopy (JMortonx et a/.l . l2002l ) yield the value about —40% (instead of —100% predicted 
by naive band picture), which can be due to both surface effects and electron correlations 
in the bulk (see Sect JIII.CIl. Transport properties of Fe,s04-based fi lms are now intensively 



studied (see, e.g., Refs. ( JEerenstein et al. 



2002 



Zhao et al. 



20051 )). In particular, a large 



magnetoresistance owin g to electron propagation through antiphase boundaries was found 
(JEerenstein et alimm . 

Unlike the Heusler alloys, magnetite is a system with a narrow 3d band and therefore 
strong correlation effects. The fact of the metal-insulator transition itself can be already 



considered as an evidence of strong electron-electron interaction (JMott 



1974 ). The influence 



of these effects on the electronic structure of Fe304 has been recently considered in Refs. 



( ICraco et al. 



2006 



Leonov et al. 



20061). 



D. Strongly magnetic half-metals with majority spin gap 



1. Anionogenic ferromagnets 



Until recently, strongly magnetic HMF with a majority spin band gap were absent. The 
chemical composition of the compou nds calculated to be half-metallic in this category was 



quite unexpected: heavy alkali oxides (JAttema et al. 



20051 ) . The magnetic moment is carried 



by complex oxygen ions, hence the above name. Besides the oxygen molecule, that has 
two unpaired electrons, the 0^~ ion occurs in the so-called hyperoxides like Rb02 and 
Cs02. These are antiferromagnetic insulators with rather low Neel temperatures. Another 
molecular ion of interest is the non- magnetic peroxide ion 0\~ . In the series molecular 
oxygen - hyperoxide ion - peroxide ion the antibonding tt orbital is progressively fllled, 
leading to the vanishing of the magnetic moment for the peroxides. Also sesquioxides exist 
which have composition between peroxide and hyperoxide. They are rather thermally stable, 
but do react with atmospheric water and carbondioxide. The analogy between the holes 
in the antibonding double-degenerate tt level and the electrons in the double- degenerate 
antibonding Cg level of the colossal magnetoresistance materials motivated a computational 
study. This yielded a HMF state with surprisingly high Curie temperatures (300K). Partial 
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explanation is the absence of superexchange in these oxides, since the mediator for it, the 
alkah ions, do not possess the required electron states in the vicinity of the Fermi level. Direct 
experimental evidences are unfortunately lacking. An indirect evidence is the cubic crystal 
structure measured down to 5K (unlike peroxides and hyperoxides) , the crystallographic 



equivalence o 



(IJansen et al. 



;he m olecular oxygen ions, the occurrence of charge fluctuations down to 5K 



1999I ). the opaque optical properties and indications of unusual widths of 



the stabili ty regions of th e sequioxides in the oxygen-rubidium and oxygen-caesium phase 



diagrams (JRengade 



190 



I 



E. Sulphides 

The spectacular developments in the area of high temperature superconductivity suc- 
ceeded by the interest in colossal magnetoresistance materials have pushed the interest in 
sulphides and selenides somewhat to the background. These materials have some advantages 
over oxides, however. Two main differences, both due to the increased metal-anion covalence 
as compared with oxides, are of importance here: the more correlated behavior of the oxides 
as well as their preference for a high-spin configuration. Sulphides often prefer a low-spin 
configuration, which make their behavior less predictable without careful computation. So, 
the sulphur analogue of magnetite, the mineral Greigite, has a magnetic moment of 2^b 
only, compared with the 4/1^ of magnetite. Consequently, it is not half-metallic. In the 
widespread mineral pyrite FeS2, iron has a non- magnetic d^ configuration, unimaginable in 
oxides. As mentioned before, magnetite shows half-metallic properties, but is at the brink 
of Mott localization: cooling down below 120K suffices to accomplish this. On the other 
hand, the much less correlated behavior of sulphospinels allows the occupation of a broad 
range of different transition metals on the octahedral and tetrahedral cation sites without 
the risk of a Mott insulating state. This does not hold for all the pyrites, however. Thus 



FeS2 is a non-magnetic semiconductor. The excellent agreement bet ween LDA calcu 



ations 



and the photoemission spectra indicate negligible correlation effects (iFolkerts et a/.l . 119871 ). 
C0S2 is a ferromagnetic metal with a Curie temperature of 122K. The magnetic moments 
were calculated function of the Hubbard U and comparison with experimental data indicated 
the importance of U (of less than 1 eV). NiS2 is a Mott insulator. NiSe2 is metallic while 
in NiSe2xS2(i-i:) the strength of the correlation effects can be adjusted by variations in the 
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composition. 

Magnetic ordering temperatures, important for maintaining the polarization of charge 
carriers at finite temperature, of oxides are usually superior to those of sulphides and se- 
lenides. 

1. Pyrites 



Saturated itinerant ferrom agnetism in the pyri te-structure system Fei_2;Coa;S2 was discov- 



ered experimentally in Ref. (IJarrett et al 



point of view in Ref. (JAuslender et al. 



19681) and discussed from the theoretical-model 



19881 ). Half- metallic ferromagnetis ni in pyrites 



was fi rst considered in band calculations by Zhao, Callaway and Hayashibaran (jZhao et al. 



19931 ). Their results for C0S2 show near the Fermi energy two completely filled t2g bands 
for the two spin directions, a partial filled eg majority band as well as a minority eg band 
just overlapping the Fermi energy. At slightly higher energy the antibonding sulphur 3p 
states are found. Clearly, cobalt disulphide is an almost half- metallic ferromagnet. Also it is 
suggested that half-metallic magnetism can_be obtained in the ternary system Fe2;Coi_2;S2, 
an idea further worked out by Mazin (JMazinl . I2OOOI ). He calculated the expected HMF re- 
gion in the phase diagram to extend from 0.2 to 0.9. A detailed study, both computational 



and experimental (IWang et al 



2OO5I ). reveals a strong dependence of the spin polarization 



at the Fermi level on the composition. Theoretically, 100% spin polarization is obtained 
for X = 0.25, whereas the maximal polarization (85%) determined with Andreev reflec- 
tion at 4.2K is obtained at x = 0.15. The polarization drops for higher concentrations of 
iron. The Fermi level is located very close to the bottom of the conduction band. This 
can lead to thermal instabilities of the half-metallicity as discussed for NiMnSb. Recently 
half -metallic properties o f pyrite-structure compounds have been reviewed by Leighton et. 



al. (JLeighton et al. 



20071 ). 



2. Spinels 

The activities in the area of half-metallicity are somewhat underrepresented. A com- 
plication in this class of compounds is that of cation ordering. The application of high 
temperatures leads to disproportionation, so long annealing at lower temperatures may be 
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required. The type of cation ordering depends on the preparation conditions. On the other 
hand, once being controlled, the cation occupancy can form a degree of freedom to achieve 
HMF materials. 

One of the compounds considered in a study on chromium chalcogenides is of interest 
here CuCr2S4. It shows an almost HM F band structure: the Fermi level is positioned 50 
meV below the top of the valence band (JAntonov et ali Il999l ). 

Sulphospi nels were a l so considered in detail in the quest for the elusive half-metallic anti- 



ferromagnet (JMin et al. 



2004 



Park et al. 



2OOII ). Mn(CrV)S4, with chromium and vanadium 



occupying the octahedral sites is calculated to fulfill all the requirements. It shows a band 
gap of approximately 2eV, while the Fermi level intersects a band of primarily vanadium 
character. The Mn moment is compensated by the moments of chromium and vanadium 
on the octahedral sites. Another sulphospinel with predicted half-metallic properties is 
(Feo.5Cuo.5)(Vo.5Tii 5)284. In this case, the metallic behavior is attributed to the atoms at 
the tetrahedral site; their magnetic moments are exactly canceled by those at the octahedral 
site. 

F. Miscellaneous 

1. Ruthenates 



The ?)d transition-elements and their compounds have been studied in much more detail 
as compared with their Ad and bd analogues. Part of the reason is that magnetism is 
expected to be favored more in the 3(i series where no d core is present. Ruthenium is a 
perfect example of the contrary. The binary and ternary oxides of this 4d transition metal 
show a rich variety of physical properties like ferr omagnetism in SrRuOs, unconventional 
superconductivity in Sr2Ru04 ( JMaeno et ali Il994j ). Here we consider the case of SrRuOa. 
Ruthenium is tetravalent in this compound, just as in RUO2. The latter compound is a non- 
magnetic metal with Ad electrons in the slightly split t2g subband. In SrRuOa, a magnetic, 
low-spin state occurs with a filled t2g majority spin band and a partially filled t2g minority 
spin band. Thus all the ingredients seem to be present for a half-metal. Calculations show 
that the exchange and crystal field splitting is not sufficient to create a band gap large 
enough to encompass the Fermi energy. Recently, it was shown that the application of the 
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LDA+U method leads to a substantial increase in band gap in conjunction with orbital 
ordering. Thus a half-metallic solution is obtained. Comparison with experiment does not 
lead to a definite conclusion. No experimental determination of U is available. The measured 
magnetic moment is more in line with the LDA results (0.8 ^ 1.6 /i^), but extrapolation 
to the high- field limit could lead to an integer magnetic moment. There is no experimental 
evidence for orbital ordering, however. 

The research on ruthenates is relatively recent and especially ternary compounds have 
not yet been investigated exhaustively. 

2. Organic half-metals 

Conducting organic materials have been an area of active re s earch since the discovery of 



electrical conduction in doped polyacethylene (IShirakawa et all 119771 ) . A surge of activities 



has resulted in applications generally referred to as plastic electronics. Recently, attempts 
have started to develop HMF suitable for these applications. Originally, the focus was 
on carbon nanotubes where magnetism was achieved by the introduction of 3d metals. 
Calcul ations were performed for 3,3 single wall carbon nanotubes with a linear iron nanowire 



inside ( iRahman et ali . |200J). Structure optimization resulted in a slightly asymmetrical 
position of the iron wire in the nanotube. The results were somewhat disappointing: the 
iron looses its magnetic moment and the overall system is semiconducting. A subsequent 
investigation of the 3,3 single wall carbon nanotube with a linear cobalt wire inside resulted 
in a HMF band structure with the band gap for the minority spin direction of order 1 eV. 
This band structure is very much like that of the iron system. The metallic properties 
are caused by the extra electron of the cobalt system, that is completely absorbed by the 
majority spin band structure. 

Another series of materials investigated is inspired by the molecule ferrocene. This is 
a so-called sandwich complex with an Fe ion between two cyclop entadienyl anions. Fer- 
rocene can be considered to be the first member of a series of so-called multiple decker 
sandwhich structures. They are formed by adding additional pairs of iron atoms and cy- 
clop entadienyl molecules. Thus the chemical structure is fundamentally different from the 
nanotubes discussed above: The latter can be thought of as two interacting wires in parallel, 
one of organic and one of metal nature. The former is characterized by a parallel stacking 
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of cyclopentadienyl (or benzene) rings coupled together by transition-metal atoms. The 
syntheses of these sys tems was shown to be possible for various vanadium benzene clusters 



(IHoshino et al. 



19951 ). The m ost promising cand idate at the moment is the one- dimensional 



manganese-benzene polymer (jXiang et all |2005| ). It is a genuine HMF with a moment of 1 



^B- The ferromagnetic ordering is much more stable than the antiferro magnetic one. This 
large difference (0.25 eV) can be traced back to the coexistence of a rather narrow and a 
rather dispersive band for the metallic spin direction, a scenario very reminiscent of the 
double-exchange model. 

III. MODEL THEORETICAL APPROACHES 

A. Electron spectrum and strong itinerant ferromagnetism in the Hubbard model 

To investigate the spectrum of single-particle and spin-wave excitations in metallic mag- 
nets we use many-electron models which permit to describe effects of inter-electron cor- 
relations. The simplest model of such a type is the Hubbard model. In the case of a 
non-degenerate band its Hamiltonian reads 



ko- 
^int = UY,C^i]Ci^c\^Cii (1) 



with U being the on-site repulsion parameter, tk the bare electron spectrum. The Hubbard 
model was widely used to consider itinerant electron ferromagnetism since this takes into 
account the largest term of the Coulomb interaction - the intra-atomic one. Despite apparent 
simplicity, this model contains a very reach physics, and its rigorous investigation is a very 
difficult problem. 

The simplest Hartree-Fock (Stoner) approximation in the Hubbard model ([T]), which 
corresponds formally to first-order perturbation theory in U , yields the electron spectrum 
of the form 



Tl 

£^k. = tk + Un_, = tk + U{- - a{S')) = tk. (2) 
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so that we have for the spin sphtting A = U{n^ — rii) = 2U{S^) and U plays the role of the 
Stoner parameter. 

Consider more strictly the case of a half-metallic (saturated) ferromagnet where n_|_ = 
n = 1 — Uq, n_ = (note that for realistic HMF the saturated ferromagnetic behavior is 
described by the generalized Slater-Pauling rule, see Sect JIILFj) . Then the spin up electrons 
behave at T =0 K as free ones, 



G^^iE) = iE-t^)-' 



(3) 



For spin-down states the situation is non-trivial. Writing down the sequence of equations of 
motion for G]f_i{E) and for the Green's function 



*kp(^) = {{S;c^-p^\cI))e,S^ = E^TCk+ci 



(4) 



and performing decoupling in spirit of a ladder approximation we obtain for the self-energy 

Uri: 



'H 



UR^{E) 



where n^ = /(^k) is the Fermi function, 



I — n 



k q 



E — tl,_r 



UJc 



(5) 



(6) 



q -- (-k-q u^q 

describes the el ectron-magnon scat t ering . This result corresponds to the Edwards-Hertz 



approximation (JEdwards and Hertz 



19731 ) 



In a more general non-saturated situation one obtains for the self-energy to second order 
in U (jlrkhin and Katsnelsonl . Il990l ) 

where Nb{uj) is the Bose function. Retaining only the magnon pole contribution to the spin 
spectral density (i.e. neglecting the spin-wave damping) we have 

SkT(i^) = ^AE J^q + ^'^+q^ (8) 

ci E - tk+ql + CJq 



Ski(i?) = f/AE 



1 + iVq - ?^k-qT 



d E - tk_q| 



UJn 



(9) 



where Uq is the magnon energy, A'"q = NB{uJq). These results are valid in the s — d model 
{U —>■ I, see below) to first order in the small parameter 1/25. Taking into account the 
relation 

{S')=So-Y.Np (10) 
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where 5*0 is the saturation magnetization one obtains for the spin-wave correction to the 
electron energy 

2' 






dkl 



a 



dh 



U{S') \ dk 



where 






aU- 



^k+q ~ ^k 



;ii) 



(12) 



The T^/^-dependence of the electron spectrum owing to magnons is weaker than the T^/^- 
dependence of the magnetization. This fact is due to vanishing of electron-magnon inter- 
action amplitude A at zero magnon wavevector, which is connected with the symmetry of 
exchange interac tion. Such a weak ening of temperature dependence of the spin splitting was 
observed in iron (ISpringfordl. 



was calculated in Ref.(Auslender and Irkhin 



1980). Theone-e 



ectron damping in the half-metallic situation 



1984al ). a Fermi-liquid-type behavior (small 



damping near the Fermi level containing high powers of temperature) being obtained. 

The problem of ferromagnetic ordering in narrow energy bands is up to now extensively 
discussed. To stabilize the ferromagnetic solution within the Hubbard model is yet another 
difficult problem. It was proved recently, that the necessary conditi ons for ferrom agnetism 



is a density of state with large spectral weigh t near the banc 



Hund's rule coupling for the degenerate case (IVoUhardt et al 



edges (lUlmkd . Il998l ) and the 



19991 ) . Real examples of sat- 



urated ferromagnetic ordering are provided by pyrite structure systems 



itinerant-electron 



1988 



Jarrett et al 



' erroinagnetism in double-de generate narrow Cg band (jAuslender et al. 



^ei_.T^Co.T,S';i with 



1968; 



Ramesha et al. 



20041 ). CMR manganites, magnetite Fe304 above 



the Verwey transition temperature, "anionic" half-metallic ferromagnets are another ex- 
amples (see Section |V]). Recently, a model of sp electro n magnetism in narrow impurity 
bands has been proposed (JEdwards and Katsnelsonl . 120061 ) which may be applicable to some 
carbon- or boron-based systems such as doped CaBg. In this model, the magnon excitations 
turn out to have higher energy than the Stoner ones. Also, T-matrix renormalization of 
the Stoner exchange parameter which decreases its value essentially in a typical itinerant- 
electron magnets is much less relevant. For these reasons, the narrow-band sp systems can 
provide an example of real "Stoner" m agnets which can have rather high Curie temperatures 
at small enough magnetization value (JEdwards and Katsnelsonl . |2006| ). According to that 
model, these ferromagnets also should be saturated. 
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Systems with strong inter-electron correlations are the most difficult case for standard 
approaches in the itinerant electron magnetism theory (band calculations, spin-fluctuation 
theories). Physically, the magnetism pict ure in this case differs substantially from the Stoner 



picture of a weak itinerant magnetism (JMoriyal . 



19851 ) since correlations lead to a radical 



recon struction of the electron spectrum — formation of the Hubbard's sub bands (IHubbard , 



1963 ) which are intimately connected with the local magnetic moments (jAuslender et al. 
1988k 



In the limit [/ — *> oo, considering the case where the number of electron n = 1~ 5 <1 {5 
is the hole concentration), the Hubbard Hamiltonian reads 



'H = Y^ ek^° k^i 



trO 
k ; 



ko- 



^o/3 



where X^f is the Fourier transform o: 



' the Hubbard operators X^ 



al3 



\ia){ip\, Ek 



(13) 



-tk. 



According to Nagaoka (JNagaokal . Il966l ). the ground state for simple lattices is a saturated 



ferromagnetic state for a low density S of current carriers ( "doubles" or "holes" in an almost 



half-fllled band). Roth flRoth . 



1969 



applied a variational principle to this problem and 



obtained two critical concentrations. The flrst one, Sc, corresponds to instability of saturated 
ferromagnetic state, and the second one, 6'^, to the transition from non-saturated ferromag- 
netic into paramagnetic state. For the simple cubic (sc) lattice, the values 6c = 0.37 and 
6'^ = 0.64 were obtained. Next, the stability of ferromagnetism was investigated within vari- 
ous approximations and methods. Most calculations for a numbe r of lattices yie l d the value 



of 6r which is close to 



expansion (jZhao et al 



.3. In particular, the Gutzwiller method (JFazekas et al\ . \l99w . t/U 



19871 ). density matrix renorma lization group approach and Quantum 



Monte-Carlo (QMC) method flLiang and Pangl . Il99,^ yielded 6c = 0.2 - 0.35 



At the same time, for the critical concentrations 6' the interval o: 



" values is broader 



and varies from 0.38 to 0.64. Irkhin and Zarubin (llrkhin and Zarubinl . I2004J . |2006| ) have 
obtained the density-of-states (DOS) pictures in a Hubbard ferromagnet with account of the 
"Kondo" scattering and spin-polaron contributions and calculated the values of the critical 
concentrations of current carriers. This approach yields a rather simple interpolational 
description of saturated and non-saturated ferromagnetism. 



The simplest "Hubbard-I" approximation for the electron spectrum (jHubbardl . Il963h 
corresponds to the zeroth order in the inverse nearest-neighbor number 1/z ("mean-fleld" 
approximation in the electron hopping). This approximation is quite non-satisfactory at 
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describing ferromagnetisin (in particular, ferromagnetic solutions are absent, except for pe- 
culiar models of bare density of states). Therefore, to treat the problems connected with 
the ferromagnetism formation in narrow bands the one-particle Green's functions were cal- 
culated to first order in \j z and in the corresponding self-consistent approximations. 



The retarded anticommutator Green's functions G'ko-(-E') = ((^^°l^°1))f; can 



3e cal- 



culated by using the equation- of motion-approach of Refs. (Ilrkhin and Katsnelson 



Irkhin and Zarubiru . 



tation one obtains ( Irkhin and Zarubin 



2004 



1988 : 



20061) with account of spin fluctuations. In the locator represen- 



2004h 



G^,,{E) = \F^,{E)-e^\-\ F^„{E) 






(14) 



with 






Xc 



+ ^k-q-a 



no + n_^ 



Gk^d-aiE 



crujn 



E^k-c 



Xc 



no + n, 



■G'k-qa(-S) 



b^^{E) = E + J2e 



^k- 



k-q 



no + n_^ 






auo, 



q/' 



(15) 



where x^q^" = ('S'qSz^) = {X"^'" Xz"^") and x^"'" = ('^^q''"''^^rq"'') are the correlation 
functions for spin and charge densities, nk^ = (X'^'^X^^). To simplify numerical calcula- 
tions, the long-wavelength dispersion law cjq = Vq^ {T> is the spin-wave stiffness constant) 
was us ed with the ma gnon spectral function K^^ioj) being average this in q. Following to 



Ref. (JNagaoka 



19661 ) the value T) = 0.666\t\ was taken for the cubic lattice and the same 



K{uj) adopted for other lattices (the choice of T> influences weakly the critical concentra- 
tion). Then a{E) and b{E) do not depend on k and can be expressed in terms of the bare 
electron density of states. In the case of the saturated ferromagnetism the expressions (I14p 
reduce approximately to the result ([6]) for [/ — > oo, 

Ski(^) = -(1 - ^o) (e r. ?'"",, . 1 ' • (16) 



E 



£^k-q + ^q, 



To write down the self-consistent approximation one has to replace in (1151) G^^^{E) — >■ 
Gka{E) and calculate n-ko- via spectral representation. In such an approach, large electron 
damping is present which smears the "Kondo" peak. 

The l/2;-corrections lead to a non-trivial structure of the one-particle density of states. In 
the non-self-consistent approach the integral with the Fermi functions yields, similar to the 
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Kondo problem, the logarithmic singularity (jlrkhin and Zarubinl . 120001 ) . For very low S a 



significant logarithmic singularity exists only in the imagina ry part of the Green's 



which corresponds to a finite jump in the density of states (jlrkhin and Katsnelson 



unction, 



1985ah . 



However, when 6 increases, it is necessary to take into account the resolvents in both the 
numerator and denominator of the Green's function, so that the real and imaginary parts are 
"mixed" and a logarithmic singularity does appear in DOS. The magnon frequencies in the 
denominators of Eqs. flTSl) result in that the singularity is spread out over the interval ciJmax 
and the peak becomes smoothed. In the self-consistent approximation the form of N^{E) 
approaches the bare density of states and the peak is smeared, even neglecting spin dynamics. 
There are no poles of the Green's function for a =| above the Fermi level at small 
6, i.e. the saturated ferroma gnetic state is pres e rved. Unlike most other analytical ap- 



2004, 



20061 ) for the one-particle Green's 



proaches, the results of Ref. (jlrkhin and Zarubinl . 
function describe formation of non-saturated ferromagnetism too, the account of longitudi- 
nal spin fluctuations Xn'^'^ being decisive for obtaining the non-saturated solution and 
calculating the second critical concentration 6[. where the magnetization vanishes. For 
6 > 6c this dependence deviates from the linear one, (S^) = (1 — 6)/2. The calculations of 



Ref. (jlrkhin and Zarubinl . 



2004 



20061) yield th e 6'^ va : 



ues w. 



1969a 



lich are considerably smaller than 



In the non-saturated state a 



the results of the spin- wave approximation (jRoth . 

spin-polaron pole occurs, so that quasiparticle states with a =| occur above the Fermi level 

(Fig.SD. 

The finite-f/ case can be also treated within the Green's function methods discussed. The 
Edwards- Hertz approximation (JSj) enables one to investigate stability of the saturated ferro- 
magnetic state only, i.e. calculate 6c. The c orresponding results are presented in Fig. [51 For 
comparison, the variational results of Refs. (jvon der Linden and Edwards! . |l99l|) are shown 
which yield a strict upper boundary for the saturated state. An agreement takes place for 
large U (far from antiferromagnetic or phase-separation instability which are not taken into 
account in the calculations) . It should be noted that DMFT yields qualitatively similar 
results ( jObermeier et ali . 119971 ). One can see that saturated ferromagnetism can occur for 
large U/\t\, and its existence at realistic U is, generally speaking, a not too simple problem. 

Now we treat the orbital-degenerate case which is more realistic for transition-metal 
compounds. Consider the many-electron system with two ground terms of the d" and d""*"^ 
configurations, r„ = {SL} and r„_|_i = {S'L'}. It is suitable to use the representation of the 
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Ferm i oper ators in terms of the many-body atomic quantum numbers ( llrkhin and Irkhinl . 



1994 



2004) 



+ ^Y" E G^L cZ^Im'^.^^ {S'L'M'^,', SLM^^) (17) 






where Gf/^ are the fractional parentage coefficients. We can introduce a further simphfica- 
tion assuming that only one of the competing configurations has non-zero orbital moment 
L = I. This assumption holds for the d^ and d^ ground-state configurations of Fe^"*" and Fe^+ 
respectively, the first configuration having zero orbital moment. A similar situation takes 
place for the CMR manganites (with d^ and d^ configuration for Mn^+ and Mn^+: due to 
the relevance of t2g — Cg crystal-field splitting the former configuration corresponds to the 
completely filled t2g band with L = 0. 

We treat the narrow-band case which should be described by a two-configuration Hubbard 
model where both conduction electrons and local moments belong to the same d-band, the 
states with n + 1 electrons playing the role of current-carrier states. After performing the 
procedure of mapping onto the corresponding state space, the one-electron Fermi operators 
for the strongly correlated states Cn^m ^^^ replaced by many-electron operators according to 
Eq. (TT7I) . Taking into account the values of the Clebsh-Gordan coefficients which correspond 
to the coupling of momenta S and 1/2 we obtain the "double-exchange" Hamiltonian 

'^ = E ti<.mgkam9kam- (18) 

kt7m 

Here we have redefined the band energy by including the many-electron renormalization 
factor, t^rnin + l){Gl'^f/{2l + 1) ^ tkm, and 



9lm= E ^^|^X,(5- 1/2, M + |;5Mm), S' = S - 1/2, 



9lm= J: ^|?^f^X,iS + 1/2, M+^;SMm), S' = S + 1/2 (19) 

where 15* Mm) are the empty states with the orbital index m, \S'M') are the singly-occupied 
states with the total on-site spin S' = S ± 1/2 and its projection M', o" = ±. We 
see that the two-configuration Hamiltonian is a generaliza tion of the narrow-band s — d 



exchange model with | /| — » cxj (double-exchange model) (jlrkhin and Irkhin 



1994 



Irkhin and Katsnelson 



2007 : 



1985a[ ): in the case where the configuration c?""*"^ has larger spin 



than the configuration d^, we have the effective "s — d exchange model" with ferromagnetic 
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coupling, and in the opposite case with antiferromagnetic couphng. In the absence of orbital 
degeneracy the model flTSj) is reduced to the narrow-band Hubbard model. 

For S = 1/2 the narrow-band s — d exchange model with |/| -^ — oo is equivalent to 
the Hubbard model with the replacement t^ -^ ik/2, so that the ferromagnetism picture 
corresponds to that described above. In a general case the criteria for spin and orbital 



instabilities (llrkhin and Katsnelson 



2005bl ) are different. It turns out that the saturated 
spin ferromagnetism is more stable than the orbital one in the realistic case S > 1/2 (e.g., 
for magnetite and for colossal magnetoresistance manganites). This means that the half- 
metallic ferromagnetic phases both with saturated and non-saturate d orbital mome nts can 
arise. The phase diagram at finite temperatures is discussed in Ref . (JEdwardsi . |2002| ) . 

In contrast with usual itinerant-electron ferromagnets, additional collective excitation 
branches (orbitons) occur in the model. Also, mixed excitations with the simultaneous 
change of spin and orbital projections exist ("optical magnons"). All these excitations can 
be well defined in the whole Brillouin zone, the daniping du e to the interaction with current 



carriers being small e nough fjirkhin 



and 



Katsnelson 



2005b|). 



The XMCD data fJHuang et al\ . l2004l ) suggest large orb ital contributions to magnetism 



in Fe304. However, more recent experimental XMCD data (JGoering et al 



20061 ) yield very 



small orbital moments in Fe304 and confirr n HMF behavior of magnetite. In any case, the 
model of orbital itinerant ferromagnetism (llrkhin and Katsnelsonl . l2005bl ) is of a general 
physical interest and can be applied, e.g., to CMR manganites. 



B. Electron spectrum in the s-d exchange model: The non-quasiparticle density of 
states 



Besides the Hubbard model, it is often convenient to use for theoretical description of 
magnetic metals the s — d{f) exchange model. The s — d exchange model was first proposed 
for tr ansition d r netals to consider peculiarities of their electrical resistivity (Vonsovsky 
1971) (I Vonsovskyl . Il974l ). This model postulates the existence of two electron subsystems: 
itinerant s electrons which play the role of current carriers, and localized d electrons which 
give the main contribution to magnetic moments. Such an assumption can be hardly justi- 
fied quantitatively for d metals, but is useful at qualitative consideration of some physical 
properties, especially of transport phenomena. At the same time, the s — d model provides a 
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good description of magnetism in rare-earth metals and their compounds with well-locahzed 
4f states. Now this model is widely used in the theory of anomalous / s ystems f i nterrn ediate 
valence compounds, heavy fermions...) as the "Kondo-lattice" model (JHewsoru . Il993l ). 

The Hamiltonian of the s — d exchange model in the case of arbitrary inhomogeneous 
potential reads 






«5 = -|;iV^ + W 



(20) 



where K-(r) is the potential energy (with account of the electron-electron interaction in 
the mean field approximation) which is supposed to be spin dependent, '^aij) is the field 
operator for the spin projection a, S(r) is the spin density of the localized-moment system, 
(5S(r) = S(r)— (S(r)) is its fluctuating part, the effect of the average spin polarization (S(r)) 
being included into V^(r). We use an approximation of contact electron- niagnon interaction 
described by the s — d exchange parameter /, 



'Hd — — 2^ JqSqS- 



(21) 



(for simplicity we neglect the inhomogeneity effects for the magnon subsystem), Sq are 
operators for localized spins, Jq are the Fourier transforms of the exchange parameters 
between localized spins. In rare earth metals the latter interaction is usually the indirect 
RKKY exchange via conduction electrons which is due to the same s-d interaction. However, 
at constructing perturbation theory, it is convenient to include this interaction in the zero- 
order Hamiltonian. 

Although being more complicated in its form, the s — d model turns out to be in some 
respect simpler than the Hubbard model ([T]) since it permits to construct the quasiclas- 
sical expansion in the small parameter \J2S. Within simple approximations, the results 
in the s — d{f) and Hubbard models differ as a rule by the replacement I —>■ U only. To 
describe e ffects of electron-magnon int e ractio n we use the formalism of the exact eigen- 



f unctions fllrkhin and Katsnelson 



1984 



20061 ). Passing to the representation of the exact 



eigenf unctions of the Hamiltonian Ti^, 






r c„ 



(22) 
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one can rewrite the Hamiltonian (120|) in the following form: 

'H = Y1 ^va(^y„Cy„ - I Yl ('^"' Ai/5|q) ^Sq4„cr„^c^^ + Hd (23) 

where 

(z/a,/ia'|q) = (/ia'|e''''k^)- (24) 

We take into account again the electron-spectrum spin splitting in the mean-field approxi- 
mation by keeping the dependence of the eigenf unctions on the spin projection. 

In the spin- wave region one can use for the spin operators the magnon (e.g., Dyson- 
Maleev) representation. Then we have for the one-electron Green's function 

Gl{E) = [E-e,^-^l{E)]-\ (25) 

with the self-energy S^(£') describing correlation effects. 

We start with the perturbation expansion in the electron-magnon interaction. To second 
order in I one has 

SUE) = 2/^5E I (- T,/^ i |q) r^A^^, Si(i^) = 2/^5E I (- i,/^ T |q) r^.^^^^^^ 

(26) 
where n^ = /(e^o-) (discussion of a more general "ladder" approximation is given below). 
Using the expansion of the Dyson equation fl25l) we obtain for the spectral density 

A.(^) = --lmGl{E) = 6{E - e,,) 

TT 

-5\E - e..)ReS^(E) - l lf^"j% (27) 



The second term in the right-hand side of Eg. (1271) gives the shift of quasiparticle energies. 
The third term, which arises from the branch cut of the self-energy, describes the incoherent 
(non-quasiparticle) contribution owing to scattering by magnons. One can see that this does 
not vanish in the energy region, corresponding to the "alien" spin subband with the opposite 
projection —a. 

Neglecting temporarily in Eq. fl26|) the magnon energy uj^ in comparison with typical 
electron energies and using the identities 

fiq E — Sim 



E ^^F (.,) = -i / dE'l^lm (H (E' -no + ^0r W) (28) 
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one derives at zero temperature 



., fiE' 



KiE) = 2l'sjdE'f-^^{v]\5{E-E'-H'^\v]) 
EUE) = 2PS ldE ^~^^^} {v[\5(E-E'-Hl)\vi) 



(29) 

Using the tight-binding model for the ideal-crystal Hamiltonian one obtains in the real-space 
representation 



Sj,^^,(E) = 2PS j dE'f{E')[--hmG\^^^{E')y^^^> 
where R, R' are lattice site indices, and therefore 

s^(E)=5:i^,.(R)|2s^,i,(i?). 

R 



(31) 

(32) 



(33) 



One can generalize the above results to the case of arbitr ary s — d exchange paramete r. 
Simplifying the sequence of equations of motion (cf. Ref. (jlrkhin and Katsnelsonl . Il984j )) 
we have for the operator Green' function 



G"(E) 



E-HI + (yl{ni - Ho'') 



1 



1 -1 



-R^iE) 



(34) 



l + cxIR''{Ey 

If we consider spin dependence of electron spectrum in the simplest rigid-splitting approxi- 
mation e^^ = El, — al {S^) and thus neglect spin- dependence of the eigenfunctions ipua (R) 
the expressions (!29l) . (!30l) are drastically simplified. Then the self-energy does not depend 
on z/: 

2PSR''{E) 



E^(E) 



l + aIR-{E)' 



RHe) = E 



nt 



1 — n^ 

R'iE)=j: ^' 



(35) 
(36) 



n E - e^l 1^ E - e^l 

If TYq is just the crystal Hamiltonian (z/ = k, Sj^o- = t]^„, tko- being the band energy), the 
expression flM|) coincides with the result for the Hubbard model after the replacement I —>■ U 
(see Sect JIILAj) . The imaginary part of S'^ (E) determines the NQP states. Description of 
such states in the Hubbard model with arbitrary U by the dynamical mean-field theory will 
be presented in Section IIV.AI 
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The expression (IMj) can be also represented in the form 



C'iE) 



E 






{E) 



-1 



(37) 



The equation (1371) is convenient in the narrow-band case. In this hmit where spin sphtting 
is large in comparison with the bandwidth of conduction electrons we have TiJ — "Hq = —2IS 
and we obtain for the "lower" spin subband with a = —sign/ 

1 -1 



C'iE) 



E-m" + 



2S 



For a periodic crystal Eq. fl38|) takes the form 



G^iE) 



E — tw_. 



This expression yields exact result in the limit / 



Gi {E) = 

with e = E + IS. In the limit / 
quasiparticles, 



2S 
-oo Eq. 



R''{E) 

2S 
^W(E) 

-^ +00, 

m = E 



(38) 



(39) 



/(^k) 



e - tb 



(40) 



k " "k 

gives correctly the spectrum of spin-down 



gUe) 



2S 



tt] 



(41) 



25 + 1 

with e = E - I{S + l),tl = [25/(25 + l)]tk- However, it does not describe the NQP states 
quite correctly, so that more accurat e expre ssions can be obtained by using the atomic 



representation (jlrkhin and Katsnelsonl . 

25 



2005b|) 



GiiE) 



25+1 



e-tl 



25 
R*(e) 



-1 



R*i^) = E 






(42) 



The Green's functions (1391) . (HUl) . fH2]) have no poles, at least for small current carrier 
concentration, and the whole spectral weight of minority states is proy ided by the branch 



cut (non-quasiparticle states) (llrkhin and Katsnelsonl. 



1983, 



1985b 



iligd) 



19901 ). For surface states 



this result was obtained in Ref . (iKatsnelson and Edwards! . 119921 ) in a narrow-band Hubbard 
model. Now we see that this result can be derived in an arbitrary inhomogeneous case. For 
a HMF with the gap in the minority-spin subband NQP states occur above the Fermi level, 
and for the gap in the majority-spin subband below the Fermi level. 

In the absence of spin dynamics (i.e., neglecting the magnon frequencies) the NQP density 
of states has a jump at the Fermi level. However, the magnon frequencies can be restored 
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in the final result, in analogy with the case of ideal crystal, which leads to a smearing of 
the jump on the energy scale of a characteristic magnon energy uJ. It should be mentioned 
once more that we restrict ourselves to the case of the usual three-dimensional magnon 
spectrum and do not consider the influence of surface states on the spin-wave subsystem. 
The expressions obtained enable us to investigate the energy dependence of the spectral 
density. 

An analysis of the electron-spin coupling yields different DOS pictures for two possible 
signs of the s — d exchange parameter /. For / > the spin-down NQP scattering states 
form a "tail" of the upper spin-down band, which starts from Ep (Fig.l) since the Pauli 
principle prevents electron scattering into occupied states. 

For J < spin-up NQP states are present below the Fermi level as an isolated region 
(Fig. 2): occupied states with the total spin S — 1 are a superposition of the states \S)\ [) 
and IS* — 1)1 t)- The entanglement of the states of electron and spin subsystems which 
is necessary to form the NQP states is a purely quantum effect formally disappearing at 
S ^ oo. To understand better why the NQP states are formed only below Ep in this case 
we can treat the limit / = — oo. Then the current carrier is really a many-body state of 
the occupied site as a whole with total spin S — 1/2, which propagates in the ferromagnetic 
medium with spin S at any site. The fractions of the states \S)\ [) and jS* — 1)| t) in the 
current carrier state are 1/(25* -|- 1) and 2S/{2S + 1), respectively, so that the first number 
is just a spectral weight of occupied spin-up electron NQP states. At the same time, the 
density of empty states is measured by the number of electrons with a given spin projection 
which one can add to the system. It is obvious that one cannot put any spin-up electrons in 
the spin-up site at / = — oo. Therefore the density of NQP states should vanish above Ep. 

It is worthwhile to note that in th e most of known HMF the gap exists for minority-spin 



states ( Irkhin and Katsnelson 



19941 ). This is similar to the case J > 0, so that the NQP 



states should arise above the Fermi energ y. For exceptiona l cases with the majority-spin gap 



such as a double perovskite Sr2FeMo06 (jKobayashi et al\ . Il998l ) and magnetite one should 



expect formation of the NQP states below the Fermi energy. 

The presence of space inhomogeneity (e.g., surfaces, interfaces, impurities) does not 
change qualitatively the spectral density picture, except smooth matrix elements. Further 
in this section we consider, for simplicity, the case of clean infinite crystal; all the tempera- 
ture and energy dependences of the spectral density will be basically the same, e.g., for the 
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surface layer. 

The second term in the right-hand side of Eq. (|27|) describes the renormahzation of quasi- 
particle energies. The third term, which arises from the branch cut of the self-energy Ti„{E), 
describes the incoherent (non-quasiparticle) contribution owing to scattering by magnons. 
One can see that this does not vanish in the energy region, corresponding to the "alien" spin 
subband with the opposite projection —a. Consider for definiteness the case / > (the case 
/ < differs, roughly speaking, by a particle- hole transformation). On summing up Eq.([2j 
to find the total DOS A^^- (E) and neglecting the quasiparticle shift we get 



N^{E) 
NliE) 



E 

kq L 



1 



it 



k+qi 



tkt)^ 



S{E - tkT 



21^ sY. 



kq 



1 + A^q - nk| 
(^k+qi — ^kT — ^c 



-5{E - tkT - ^q) 



(43) 



The T^/^- dependence of the magnon contribution to the residue of the Green's function, i.e. 
of the effective electron mass in the lower spin subband, and an increase with temperature of 
the incoherent tail from the upper spin subband result in a strong temperature dependence 
of partial densities of states N„{E), the corrections being of opposite sign. At the same 
time, the temperature shift of the band edge for the quasipartic l e states is pr oportional to 



1983 



1984 



1985bh 



T^/^ rather than to magnetization ( Irkhin and Katsnelsonl . 

It is worthwhile to note that there exists a purely single-particle mechanism of the gap 
filling in HMF which is due to relativistic interactions. Specifically, one should take into ac- 
count spin-orbit coupling effects which connect the spin-up and spin-down channels through 
the angular momentum 1. The strength of this interaction is proportional to the spatial 
derivatives of the crystal potential V{j): Vso oc W (1 ■ s), off-diagonal elements V^q being 
non-zero. For HMF with a gap in the minority-spin (spin-down) channel one could construct 
the wave function for spin-down electrons within perturbation theory, so that the DOS in 
the gap has a square dep endence on the spin-orbit coupling strength, 5n^^{E) oc (V^J)^ 
(JMavropoulos et all 120041 ). There is an obvious qualitative distinction between the many- 
body and spin-orbit contribution in the minority spin channel; besides that, that the latter 
is orders in magnitude smaller and very weakly temperature-dependent. Fo r further discus- 
sions of the spin-orbit effects in HMF, see Ref. ( iPickett and Eschrigl . 120071 ). 



The exact solution in the atomic limit (for one conduction electron), which is valid not 
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only in the spin-wave region, but for arbitrary temperatures, reads (jAuslender et all Il983l ) 



^ ' 2S+1 E + IS 2S+1 E-I{S+1) ^ ^ 

In this case the energy levels are not temperature dependent at all, whereas the residues are 
strongly temperature dependent via the magnetization. 

Now we consider the case T = K for a finite band filling. The picture of N{E) in 
HMF (or degenerate ferromagnetic semiconductor) demonstrates strong energy dependence 
near the Fermi level (Figs. 1,2). If we neglect magnon frequencies in the denominators of 
Eg. (1431) . the partial density of incoherent states should occur by a jump above or below the 
Fermi energy Ep for / > and / < respectively owing to the Fermi distribution functions. 
An account of finite magnon frequencies cUq = Vq"^ {T> is the spin wave stiffness constant) 
leads to smearing of these singularities, N^ai^Ep) being equal to zero. For \E — Ef\ <^ cJ we 
obtain 

N.^{E) 1 ^-^^'/%(^(^_^^))^^^,ig^j (45) 



N^{E) 2S 



UJ 



(a = ± corresponds to the spin projections t, |). With increasing \E — Ep\, N_a/Na tends 
to a constant value which is of order of P within the perturbation theory. 
In the strong coupling limit where |/| ^ oo we have from (H3l) 



^^ = ^0{a{E - Ep)\ \E - Ep\ » ^ (46) 

In fact, this expression is valid only in the framework of the l/2S'-expansion, and in the 
narrow-band quantum case we have to use more exact expressions (1401) . (1421) . 

To investigate details of the energy dependence of A^ {E) in the broad-band case we 
assume the simplest isotropic approximation for the majority-spin electrons. 

Provided that we use the rigid splitting approximation tkj, = tkt + ^ (^ = 21 S, / > 0), the 
half-metallic situation (or, more precisely, the situation of degenerate ferromagnetic semi- 
conductor) takes place for A > Ep- Then qualitatively the equation fj45|) works to accuracy 
of a prefactor. It is worthwhile to note that, formally speaking, the NQP contribution to 
DOS occurs also for an "usual" metal where A < Ep. In the case of small A there is a 
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crossover energy (or temperature) scale 



T* =V{m*A/kFf 

which is the magnon energy at the boundary of Stoner continuum, T* 
At \E — Ef\ ^ uj the equation (H3l) for the NQP contribution reads 



(48) 



uJiA/Epf <^uJ. 



6N^{E) oc 



1 



hi 



1 + J{E-Ef)/T* 



{E - Ep) /T* 



eiyE-Ep). 



(49) 



\E - Ep) IT* 

At |ii^ — Ep\ <^ T* this gives the same results as above. However, at T* ^ \E ^ Ep\ <^ uJ 
this contribution is proportional to 



'{E — Ep) /T* and is negative (of course, the total 
DOS is always positive). This demonstrates that one should be very careful when discussing 
the NQP states for the systems which are not half-metallic. 



The model of rigid spin splitting used a bove is in fact not app 



wher e the gap has a hybridization origin (jde Groot et all Il983bl : llrkhin and Katsnelsonl . 



icable for the real HMF 



19941 ). The simplest model for HMF is as follows: a "normal" metallic spectrum for majority 
electrons fH7j) and the hybridization gap for minority ones. 



tki - El 



i ( a + sign (a) ^a + A^ 



(50) 



Here we assume for simplicity that the Fermi energy lies exactly in the middle of the hy- 
bridization gap (otherwise one needs to shift ^k ^ ^]s. + Eq — Ep in the last equation, Eq being 
the middle of the gap). One can replace in Eq. (l43!) ^k+q by Vkq, Vk = k/m*. Integrating 
over the angle between the vectors k and q we derive 



X^ - (X^ + 1)'^/^ + 1 +x 



(51) 



\ V^k+qi -^kT -^q/ / VpqAK?, 

where angular brackets stand for the average over the angles of the vector k, X = kpq/m*A. 
Here we do have the crossover with the energy scale T* which can be small for small enough 
hybridization gap. For example, in NiMnSb the conduction band width is about 5 eV and 
the distance from the Fermi level to the nearest gap edge (i.e. indirect energy gap which is 
proportional to A^) is smaller than 0.5 eV, so that {A/Ep)'^ < 0.1. 
For the case < E — Ep -C u one has 



mE)^b 



E-Ep 



b{y) 



2 r 
5 



y'^'-ii + yr^ + i 



5/2 



y + y 



3/2 



y3/2, y < 1 

y, y > 1 



(52) 
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Thus the behavior N^(E) oc (£" — Ep) ' takes place only for very small excitation energies 
E — Ep <^ T* , whereas in a broad interval T* <^ E — Ep ^ uJ one has the linear dependence 
Ni{E) (xE-Ep. 

C. The problem of spin polarization 



Functionality of devices that exploit charge as well as spin degrees of freedo m de; 



pend s in a crucial way on the behavior of the spin polarization of current carriers (jPring . 
19981 ). Unfortunately, many potentially promising half- m etallic systems exh ibit dramatic 



decrea se in the spin polarization. Crystal imperfe ctions (jEbert and Schutd . Il99ll ). inter- 



faces (Ide Wijs and de Groot , 



2OOII ). and surfaces (JGalanakid . 120031 ) constitute important 



examples of static perturbations of the ideal periodic potential which affect the states in the 
half-metallic gap. 



In addition, several other depolari zation mechanisms w e re su 



on magnon and phonon ex citations ( JDowben and Skomski 



y gested that are 



2003 



2004 



Skomski 



Skomski and Dowben 



3ased 
2OO7I : 



2002). These papers extend the view of spin-disorder as random 



inter-atomic exchange fields and claim that disorder rotates locally t he spin di r ection 
and thus modifies the local magnetic moment a nd spin polarization (jitoh et all 12000 : 



MacDonald et al. 



1998 : 



Orgassa et al 



1999 



2OOOI ) . The coupling between atomic moments 



can be treated in terms of the Heisenberg-type interactions (see Sect JVI. Al and |VI.B|) . The 
sign and magnit ude of the e xchan ge constants determines whether the spin structure is 
collinear on not (jSandratskii 1200 ll ) . 

Simple q ualitative considerations (Edwards! . 119831 ) . as well as direct Green's functions 



calculations (jAuslender and Irkhinl . 



1984a 



3 



I) for ferromagnetic semiconductors, demonstrate 



that spin polarization of conduction electrons in spin-wave region is proportional to magne- 
tization 



P = 



N^-N 



i 



Po{S^)/S 



(53) 



A weak ground-state depolarization 1 — Pq occurs in the case for the empty conduction 
band where / > 0. As discussed in the previous section, in the case of the Fermi statistics 
of charge carriers (degenerate ferromagnetic semiconductor and HMF) the NQP states at 
T = exist only below Ep for majority-spin gap (/ < for the case of semiconductors) and 
only above Ep for minority-spin gap (/ > for the case of semiconductors). 
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Spin-resolved photoelectron s pectra for raagnetite slightly above the Verwey transition 



point have been measured in Ref . (iMortonx et a/.U2002l ). a negative polarization about —40% 
being found near the Fermi energy. The strong deviation from —100% polarization can be, at 
l east partially, related to NQ P states. Since according to the electronic structure calculations 



(^Yanase and Siratori 



19841 ) magnetite is a HMF with majority-spin gap, the NQP should 



exist below the Fermi energy and thus be relevant for photoelectron spectroscopy. Since 
electron correlations in Fe304 are quite strong, the spectral weight of the NQP states should 
be considerable. Of course, the photoemission is a surface-sensitive method, and it is not 
quite clear to what degree these data characterize the electronic structure of bulk Fe304. 

An instructive limit is the Hubbard ferromagnet with infinitely str ong correlations U = oo 
(IT511 and ele ctron concentration n < 1 The DOS calculations yield (jlrkhin and Katsnelsonl . 



1983 



1985bf ) 



N^{E)=J2f{t^+^)6{E-t^^^ 



UJr, 



kcr 




E — Ep ^ ^max 

E < Ep 



(54) 



A schematic density of states is shown in Fig. [29k (a more realistic picture is presented in 
Fig. m see also Ref. drkhin and Katsnelsonl . l2005bl )). The result (15^ has a simple physical 
meaning. Since the current carriers are spinless doubles (doubly occupied sites), the electrons 
with spins up and down may be picked up with an equal probability from the states below 
the Fermi level of doubles, so that these states are fully depolarized. On the other hand, 
according to the Pauli principle, only the spin down electrons may be added in the singly 
occupied states in the saturated ferromagnet. 



The behavior P{T) ~ (5*^) is qualitativ ely confirmee 



sion from ferromagnetic se miconductors (Kisker et al 



half-metallic Heusler allovs ( Otto et al. 



2003 



Skomski and Dowben 



jy ex perimental data on field emis- 



19781 ) and t ransport properties of 



19891 ) . An attempt was used (iDowben and Skomskil . 



20021 ) to generalize the result (155]) on the HMF case (in fact, 
using qualitative arguments which are valid only in the atomic limit, see Eq. (H^ ). However, 



the situation for HMF is more complicated. We focus on the magnon contribution to DOS 
( l43l) and calculate the function 



A = E 



-5{Ef - tkt) 



(55) 



kq (^k+qi - tkT - ^q) 

Using the parabolic electron spectrum tk| = k"^ /2m* and averaging over the angles of the 
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vector k we obtain 

% q (g*) -g^ 

where p = N-^{Ef,T = 0), we have used the condition q <^ kp, q* = m*A/kF = A/fi?, 
where A = 2 |/| S* is the spin sphtting. In the ferromagnetic semiconductor we have, in 
agreement with the quahtative considerations presented above: 

A = ^^poc(-) p (57) 

Further on we consider the spectrum model (jSj), (1501) where the gap has a hybridiza- 
tion origin. At T ^ T* we reproduce the result ( 1571) which is actually universal for this 
temperature region. At T* ^ T ^ a; we derive 

A = 5: 2PSNM^^) -\ cxq^Y.— ^ -J^T In - (58) 

^ ^ 3vFqA ^ q Tc' ^ 

This result distinguishes HMF like the Heusler alloys from ferromagnetic semiconductors 
and narrow-band saturated ferromagnets. In the narrow-band case the spin polarization 
follows the magnetization up to the Curie temperature Tc- 

For finite temperatures the density of NQP states at the Fermi energy is proportional to 



oo 



N{Ef)^I du-^P— (59) 

JO smh(a;/i ) 

The filling of the energy gap is very important for possible applications of HMF in spintronics: 
in fact HMF have deciding advantages only provided that T <^Tc- Since a single-particle 
Stoner-like theory leads to much less restrictive (but unfortunately completely wrong) in- 
equality T ^ A, the many-body treatment of the spin-polarization problem (inclusion of 
collective spin-wave excitations) is crucial. Generally, for temperatures which are compa- 
rable with the Curie temperature Tc there are no essential difference between half-metallic 
and "ordinary" ferromagnets since the gap is filled. 



1989. 


1994) 



The corresponding symmetry analysis was performed in Ref. (llrkhin et al. 
for a model of conduction electrons interacting with "pseudospin" excitations in ferroelectric 
semiconductors. The symmetrical (with respect to Ep) part of N{E) in the gap can be 
attributed to smearing of electron states by electron-magnon scattering; the asymmetrical 
( "Kondo-like" ) one is the density of NQP states owing to the Fermi distribution function. 
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Skomski and Dowben ( Dowben and Skomski 



2003 




2004 



Skomski 



Skomski and Dowben 



2007 : 



20021 ) investigated spin mixing effects for NiMnSb by using a 
simple tiglit-binding approximation. 

Fig. [H] sliows also a schematic comparison between this approximation and many-body 
results. In the tight-binding approach, the distortion of the spin-up and spin-down DOS is 
presented by the dark regions. The spin mixing gives a non-zero symmetric Ni{E), DOS 
being weakly m odified by therma l fluctuations. 



Itoh et al. ( Itoh et al. 



20001 ) calculated the polarization for a ferromagnet-insulator 



magnetic tunnel junction with and without spin fluctuations in a thermally randomized 
atomic potential. The results indicate that the effect of spin fluctu ations is sign i ficant . 



The idea of spin fluctuations was further developed by Lezaic et al. (JLezaic et al 



20061) 



by considering the competition between hybridization and thermal spin fluctuation in the 
prototype HMF NiMnSb. 

Figure [9] shows the sublattice susceptibilities computed within a generalized Heisenberg- 
type Hamiltonian. These results demonstrate that the Ni-sublattice magnetic order is lost 
already at 50K (this effect is a consequence of weakness of the couplin g between the Ni 



mome nts and the neighbor atoms), however, neutron diffraction data (JHordequin et al. 



1997a 



b|)does not support this scenario. The right hand side of Fig. [9] presents the polariza 



tion P{T) calculated in a disordered local moment (DLM) approach (JGyorffy et a/.l . ll985l ). 
representing the system at finite temperatures in a mean-field way. It is concluded that the 
thermal collapse of the polarization is connected with a change in hybridization due to the 
moment fluctuation, the effect of non-coUinearity being much milder. 

The issue of finite-temperature spin polarization in half-metals remains an open question. 
Magnons play a crucial role, independently of the theoretical approach. Besides that, a role 
of phonon modes in the the many-sublattice half-metals is not excluded. Nevertheless, a non- 
trivial contributions to physics of half-metallic ferromagnets comes from genuine many-body 
effects. The corresponding first-principle calculations will be presented in SectJVl 



D. The tunneling conductance and spin-polarized STM 



An useful tool to probe the s pin polarization and iion-quasiparticle states in HMF is pro- 
vided by tunneling phenomena ( jAuslender and Irkhinl . Il985bl : 



Auth et al 



20031 : iBratkovskvl . 
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19981 : iGercsi et all l2006l : ISakuraba et all l2006al . 120071 . l2006bl ) , i n particular, by Andre ev re- 
flection spectroscopy for a HMF-superconductor tunnel junction (JTkachov et a/.l . l200ll ). The 
most direct way is probably the measurement of a tunnel current between two pieces of HMF 
with the opposite magnetization directions. 

Let us explain in a simple qualitative way why the NQP states are important for the 
tunneling transport. To this aim we consider the above-discussed narrow-band saturated 
Hubbard ferromagnet where the current carriers are the holes in the lowest Hubbard band 
and the NQP states provide all the spectral weight for the minority spin projection . Suppose 
we have a tunnel junction with two pieces of this ferromagnet with either parallel (Fig. [29b) 
or antiparallel (Fig. 12^ ) magnetization directions. From the one-particle point of view, 
the spin-conserving tunneling is forbidden in the latter case. However, in the framework of 
many-particle picture the charge current is a transfer process between an empty site and 
a single-occupied site rather than the motion of the electron irrespective to the site like 
in the band theory and therefore the distinction between these two cases (see Fig. [29l) . is 
due only to the difference in the densities of states. It means that the estimations of the 
tunneling magnetoresistance based on simple one-electron picture is too optimistic; even for 
antiparallel spin orientation of two pieces of the half-metallic ferromagnets in the junction for 
zero temperature the current is not zero, due to the non-quasiparticle states. More exactly, 
it vanishes for zero bias since the density of NQP states at the Fermi energy equals to zero. 
However, it grows with the bias sharply, having the scale of typical magnon energies, i.e., 
millivolts. 



Formally, we can consider a standard tunneling Hamiltonian (see, e.g., Ref. (jMahanl . 



1990)) 



U = HL + nR + E(^kp4tCpi + h.c.) (60) 

kp 

where 'Hl,r are the Hamiltonians of the left (right) half-spaces, respectively, k and p are 
the corresponding quasimomenta, and spin projections are defined with respect to the mag- 
netization direction of a given half-space (spin is supposed to be conserving in the "global" 
coordinate system). Carrying ou t standard ca lculations of the tunneling current 1 in the 
second order in Tkp one has (cf. (JMahanl . Il990l )) 



J oc 5] iTkpHl + iVq - /(Vq)][/(^k) - /(tk + eVmeV + t^- Vq + ^q) 

kqp 



(61) 
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where V is the bias voltage. 

A ver y efficient new experimental method is sp i n-polarized scanning tunne ling microscopy 



(STM) flHeinze et al. 



200G 



Kleiber et al. 



2000 



Wiesendanger et al. 



1990l ) which enables 



one to probe directly the spectral density with spin resolution in magnetic systems. The 
spin-polarized STM should be able to probe the NQP states via their contribution to the 
differential tunneling conductivity dla/dV oc A^^r {gV) (jlrkhin and Katsnelsonl . l2006l ). Note 



that the value N„{ eV) vanishes for \eV\ < hu n, where huJo is an anisotropy gap in the 
magnon spectrum (llrkhin and Katsnelson . |2002| ). which is small, but could be changed by 
suitable substitution (jAttema et all l2004l ). Keeping in mind that ferromagnet i c sern icon- 
ductors can be considered as a peculiar case of HFM fjrkhin and Katsnelsonl . Il983l ). an 
acco unt of NQP states can be important for proper description of spin diodes and transis- 



tors fJFlatte and Vignale 



2001 



Tkachov et al 



I 



200l|) 



The above formulas are derived for the usual one-electron density of states at Ep, which 
is observed, say, in photoemission measurements. However, the factors which are present in 
the expression for the tunneling current do not influence the temperature dependence, and 



therefore thes e results are valid for spin po 



bias in STM (jlrkhin and Katsnelsonl . 



2006 



arization from tu nneling conductance at zero 



Ukraintsev 



19961 ). Unlike the photoemission 



spectroscopy which probes only occupied electron states, STM detects the states both above 
and below Ep, depending on the sign of bias. 

One should keep in mind that sometimes the surface of HMF is not half-inetallic : in 
particular, this is the case of a prototype HMF, NiMnSb ( Ide Wijs and de Grootl . l200ll ). In 
such a situation, the tunneling current for minority electrons is due to the surface states only. 
However, the NQP states can be still visible in the tunneling current via the hybridization 
of the bulk states with the surface one. This leads to the Fano antiresonance picture which 
is usually observed in STM investigations of the Kondo effect at metallic surfaces. In such 
cases the tunneling conductance will be proportional to a mixture of A^^ and the real part 
of the on-site Green's function, L^-. Surprisingly, in this case the effect of NQP states on 
the tunneling current can be even more pronounced in comparison with the ideal crystal. 
The reason is that the analytical continuation of the jump in N„ {E) is logarithm; both 
singularities are cut at the energy uJ ; neve rtheless, the energy dependence of L„ (E) can be 



pronounced (llrkhin and Katsnelson 



20061) 



STM measurements of electron DOS give also an opportunity to probe bosonic excitations 
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interacting with the conduction electrons. Due to electron-phonon couphng, the derivative 
dN(j {E) /dE and thus d'^X^r {V) /dV"^ at eV = E have peaks at the energies E = ±uJi 
corresponding to the peaks in the phonon DOS. According to above treatment (see, e.g., 
Eq. fHSl) . the same effect should be observable for the case of electron- magnon interaction. 
However, in the latter case these peaks are essentially asymmetric with respect to the Fermi 
energy (zero bias) due to asymmetry of the non-quasiparticle contributions. This asymmetry 
can be used to distinguish phonon and magnon peaks. 



T hermoelectric power in 



Refs. ( jMcCann and Fal'ko 



: he tunnel situation was theoretically investigated in 



2002 



McCann and Fal'ko 



20031 ). The relative polarizations of 



ferromagnetic layers can be manipulated by an external magnetic field, and a large difference 
occurs for a junction between two ferromagnets with antiparallel and parallel polarizations. 
This magnetothermopower effect becomes giant in the extreme case of a junction between 
two half-metallic ferromagnets, since the thermopower is inversely proportional to the area 
of the maximal cross section of the Fermi surface of minority electrons in the plane parallel 
to the interface. One has a strong polarization dependence of the thermopower 

Qap = OMkB/e, Qp oc A;|T/(eE^) (62) 

This result is independent of temperature and of the specific half-metallic material, and it 
represents a giant magnetothermopower effect, AQ ~ Qap = — 55/iV/K. 

E. Spin waves 

Unlike the Stoner theory, the Hubbard model and other model with electron correlations 
enable one to describe spin- wave excitations in an itinerant ferromagnet. This was already 



done in the old approaches based on the random phase approximation (RPA) (JHerringl . 



19661 ) . To discuss related approaches we present the interaction Hamiltonian in terms of the 
spin density operators 

^int = ^ E cLck. - ^ E(^-c,^q + S^SZ^) (63) 

^ ka -^ q 

where S^ are the Fourier components of spin density operators. The first term in (1631) yields 
a renormalization of the chemical potential and may be omitted. Writing down the sequence 
of equations of motion for the spin Green's function 

G^{u) = {{S^\SZ^)). (64) 
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one derives (jlrkhin and Katsnelsonl . Il990l ) 



where 



*-^ql^J — o f \ TT 

00 — \lq[UJ) — TTq^U;) 



fiq(^) = UY,- ^^"^ ^ ^ -{n^^ - rik+qj 



k Wq — t]i + A — to 

and the self-energy vr describes corrections to RPA. Unhke the standard form of RPA, 



Gq(c.) = ^-^ 



UU^iu) 



(65) 



(66) 



(67) 



nq(^) = E ''"^"^"^"^ , 

the equivalent representation ( l65l) yields explicitly the magnon (spin-wave) pole 



^q - 



^q(0)=E^kq^k. 



ko- 



where ^^q is given by (TT2!) . Expanding in g we get cUq = Va/sqaQp where 



^ k 



<9^ik . , 1 dtk dtk 



-(rikt +^ki) - 



(nkt - nki) 



(68) 



(69) 



(70) 



dkadkf^ A 9fca dkf^ 

Eq. [70] are spin- wave stiffness tensor components. For a weak ferromagnet (A <^ Ep, U) we 
have V (X A. The magnon damping in the RPA is given by 

9nkt^ 



I'^K^) 



-Im f2q(c<j) = TT 



f/Aa;X:(- 



9t 



kT 



6{UJ - tk+qi + t 



kT> 



(71) 



'q 



iSK^n) 



7TUAuj^N^{EF)Ni{EF)e{uj^ - 00-) (72) 

with 6{x) being the step function. Here oo^ = oo{q*) is the threshold energy which is de- 
termined by the condition of entering into the Stoner continuum (decay into the Stoner 
excitations, i.e. electron- hole pairs), q* being the minimal (in k) solution to the equation 



^k+q*i — ^kT — Ep 



(73) 



The quantity oo^ determines a characteristic energy scale separating two temperature regions: 
the contributions of spin waves (poles of the Green's function [65] dominate at T < c<j_ , and 
those of Stoner excitations (its branch cut) at T > a;_. 
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Although the formal expressions in the s — d exchange model are similar, presence of 
two electron subsystems leads to some new effects, in particular, to possible occurrence 
of the "optical mode" pole uo ^ 2|/|S'. The problem of the optical mode formation and 



i ts damping was investigated in app 



(lAuslender and Irkhin 



et al (JKaplan et al. 



1984al . Il985a : 



ication to degenerate ferromagnetic semiconductors 



Irkhin 



1987 



Irkhin and Katsnelsonl . Il985al ) . Kaplan 



2OOII ) performed exact diagonalization studies of the double exchange 



model which indicate the existence of continuum states in the single-spin-flip channel that 
overlap the magnons at very low energies (of order 10~^ eV) and extend to high energies. 
This picture differs dramatically from the prevalent view, where there are the magnons, plus 
the Stoner continuum at the high-energy scale, with nothing in between. Peculiarities of 
magnons in HMF, especially, in the c oUossal magnetoresistanse materials, have been recently 
reviewed in Ref . ( IZhang et all 120071 ) . 

In the case of weak ferromagnets, the contribution of the branch cut of the spin Green's 
function may be approximately treated as that of a paramagnon pole at imaginary cj, and 
we obtain 



q* = kpi — kpi , u_ = V{k,p^ — kpif ~ A^ ~ T^jEp 

Since g* is small, we have at small g > g*, instead of Eq. fl72l) . 



(74) 



7i^H^q) 



f/Acjfio,._ =.^2 



q Att 



{my = A/q 



(75) 



with Qq the lattice cell volume. The estimation ( 17^ holds also for the s — d{f) exchange 
model with the indirect RKKY-interaction where V ~ Tc/S ~ PS/ Ep- 

The damping at very small q < q* (where ( I7T]) vanishes) is due to the two-magnon 
scattering processes. To consider these we have to calculate the function vr to leading order 
in the fluctuating part of the Coulomb interaction. Writing down the equation of motion 
for the Green's function ( !64l) we obtain 



7rq(^) = E(4q)'[5(k T, k + q - p T, ^p - ^) + 5(k + p i, k + q i, o^p - a;) 
pk 

- 5(k + pi,kT,a;p)-5(k + qi,k + q-pi,p,a;p)] (76) 

where 
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B[k. a , ko", cijj = 

The magnon damping needed is given by the imaginary part of (ff6l) . 



(77) 



7^ (cu) = vr ^(A^)2(nk. - rik+q-p.) [iVp - Nb{uJp - tu)] 5{tu + tk - tk+q-p - cUp) (78) 



kptr 



Integration for the isotropic electron spectrum gives (jlrkhin and Katsnelson 



1990|) 



7?)(-) 



"o 



1 VA-2 X 

127r3 4(5^)2 4^^'^ 



J tUq/35 

I (T/4) (ln(T/c.q) 



+ 






(79) 



These resul ts were obtained by Silin and Solontsov within the ph enomenological Fermi- 



1984a 



liquid theory (jSilin and Solontsovl . Il984j ) and Auslender an d Irkhin (lAuslender and Irkhinl . 



1985al ) within the s — d exchange model. Golosov ( 



results of Refs. ( lAuslender and Irkhin 



1984al . 



1985a : 



Golosovl. 



20001) r eprodu ced the 



Irkhin and Katsnelson 



1985aJ) within 



the l/2S'-expansion and performed numerical investigations of the magnon spectrum and 
damping in the limit of large |/| (double-exchange situation) in application to colossal mag- 
netoresistance compounds. 

Real part of (176!) describes the temperature dependence of the spin stiffness owing to 
two-magnon processes (besides the simplest T^ -contribution which occurs from the temper- 
ature dependence of the Fermi distribution functions in (165|) . The spin- wave contribution 
connected with the magnon distribution functions is proportional to T. More interesting is 
the non-analytical many-electron contribution owing to the Fermi functions: 



6n 



aP 



A{S')'- 



E 

pk 



dka dk, 



P 



^ki(l-rik-pT) 
^k — ik-p — "^p 



?^k+pi ( 1 - ^ki ) ^k+pi ( 1 - ^kT ) '^ki ( 1 - ^k-pi ) 



U 



h 



h 



h 



U 



U 



UJr- 



i-k+p — i-k — ^p i-k+pj, — i-kt ~ ^p f-ki — ''k-pt ^p 

Performing integration for parabolic spectra of electrons and magnons yields 



(80) 



5V{T) = 


( nQoT \ 
V12(S^)m*; 


' 1 
D 


with 







Y: NI{Ef) In — - 2N,{Ep)N,{Ef) In "^^^^-^T) (gl) 



uj^ = V{kF] ± kpif , N„{Ef) = m*nokF/27r^ (82) 
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It should be noted that the correction (IHTj) dominates at low temperatures over the above- 
mentioned T^-correction, which demonstrates an important role of corrections to the RPA 
approximation. Unfortunately, the T^ In T-term has not yet to be considered at analyzing 
magnon spectra of ferromagnetic metals. We see that temperature dependences of spin- 
wave characteristics in conducting magnets differ considerably from those in the Heisenberg 
model. 

F. Magnetization and local moments 

To treat the problem of magnetic moments in the Hubbard model we consider corrections 
to the magnetization {S^). We have 

(^^) = ^-E(^q^q)-(^^T^a) (83) 

^ q 

The first average involved in fl83p is calculated from the spectral representation of the RPA 
Green's function |65l 

{SZ^S+) = 2SoN^ (q < q*) (84) 

1 /■- , iVB(u;)7«(cu)(A-cu)/f/ 



^ •^-'^ [uj - Re fiq(cu)]2 + [7V (^)]^ 
In contradiction with fl83l) . fl8^ . in the true Bloch spin- wave contribution to magnetization 
every magnon should decrease (S*^) by unity. The agreemen t may be restored by a l lowin g 



not only the magnon pole, but also branch cut contributions (llrkhin and Katsnelson 



1990). 



In the semi-phenomenological manner, it is convenient to introduce "magnon" operators 
which satisfy on the average the Bose commutation relations: 

6q = (25o)-^/^S+ , fet = (25o)-^/^5lq (86) 



Then we have 



^m = - E( KK) = j^ Eis-^s^) (87) 



q V — u/ q 

Performing integration over u in (l85l) at T=0 K we obtain 



1 ^_ 7^^) W 
S{S') =_i^ia_in_ (88) 



TT V ^q ^q 
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with W being the bandwidth. This contribution of the order of U^\n(W/uj+) describes 
zero-point decrease of the magnetization due to the ground-state magnon damping which is 
owing to the Stoner excitations. For parabohc electron and magnon spectra, neglecting the 
damping in the denominator of (jHSD we obtain at low temperatures T < uj^ the dependence 
6{S^)ci oc U'^{T /uj^Y- For a weak ferromagnet, the temperature correcti on is proportio nal 
to {T/TcY, in agreement with the self-consistent renormalization theory (JMoriyal . Il985l ). It 
should be stressed that the T^-correction obtained is much larger than the Stoner contri- 
bution of the order of {T / EpY. The spin-wave corrections to the local magnetic moment 
at a site (S^) = (3/4')('n — 2( ni^ni\)) at low T <ti uo^ have a weak dependence —{T/TcY^'^ 
(llrkhin and Katsnelsonl . [19901) • This justifies their neglecting in the above discussion of the 
magnetization (1551) . 

At high T > uj^ the damping in the denominator dominates at small q in the case of a 
weak ferromagnet. Taking into account (!75|) we obtain from (185|) 

A /•<» r°^ uAqdq / T ^^/^ 



Hsz^s^) 



dujNB{uj) 



oc 



(89) 



ttU J -oo "'"'"'"' Jo {Vq^y + A^uj^q^ " KEp 

Thus we get from (!83l) the T^'^-contribution to the magnetization, which agrees with the 
result of the phase transition scaling theory near T = Tc- For a ferromagnet with well- 
localized magnetic moment s the damping may be neglect ed and we derive a small correction 
5{S')ei OC -Puj_\n{T/uj^) Jlrkhin and Katsnelsonl . Il996[ ). 

Now we discuss a more realistic situation in HMF, in particular in the Heusler 
alloys. These compounds d emonstrate h i gh va l ues of the saturation magnet i zation 



and Curie temperatu r e (see ( 



Irkhin and Katsnelsonl . 



1994 



Fecher et al. 



Ozdogan et al. 



2006 



Galanakis and Dederichs(eds.) 



2005 



20061 )). The strong ferromagnetism is mainly 



due to local moments of well-separated Mn atoms. On the other hand, the highest magnetic 
moment (6/2^) and Curie-temperature (HOOK) in the classes of Heusler cora pounds as well 
as half-metallic ferromagnets was revealed for Co2FeSi ( IWurmehl et am2005l ). It was found 
empirically that the Curie temperature of Co2-based Heusl er compounds can be estimated 



from a nearly linear dependence on the magnetic moment (jFecher et al. 



20061). 



High spin polarization and magnetic mo ment of half-rnetal. 



within the generalized Slater-Pauling rule (jFecher et al. 



2006 



i c ferromagnets can be tr eated 



Galanakis et al. 



2002bl ). Ac 



cording to the original formulation by Slater and Pauling, the magnetic moments m of 3d 
elements and their binary compounds can be described by the mean number of valence elec- 
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trons ny per atom. A plot of m versus maK netic va l ence r n(nM) is called the generalized 



198J). According to Hund's rule it 



Slater-Pauling rule, as described by Kiibler (iKiiblerl . 
is often favorable that the majority d states are fully occupied (n^f = 5). Starting from 
m = 2n^ — ny, this leads to the definition of the magnetic valence as um = 10 — ny, so that 
the magnetic moment per atom is given by ?Ti = um + 2nsp|. 

In the case of localized moments, the Fermi energy is pinned in a deep valley of the 
minority electron density. This constrains n^j to be approximately 3, and m = ny — 6 — 2rasp|. 
HMF are supposed to exhibit a real gap in the minority density of states where the Fermi 
energy is pinned. Then the number of occupied minority states has to be an integer. Thus, 
the Slater-Pauling rule will be strictly fulfilled with the spin magnetic moment per atom 
m = ny — 6. The situation for the HMF and non-HMF full Heusler alloys is shown in FigJTUl 

For ordered compounds with different kinds of atoms it may be more convenient to 
consider the total spin magnetic moment Mf of all atoms of the unit cell. This quantity scales 
with the number of valence electron Zf. Mt = Zt — 18 for the half-Heusler and Mt = Zj— 24 for 
the fuU-Heusler alloys. Thus in both types of compounds the spin magnetic moment per unit 
cell is strictly integer for the HMF situation. On the other hand, for alloys with non-integer 
site occupancies like the quaternaries X2Yi_2,.Y^Z the moment may become non-integer 
depending on the composition, even for the HMF state. The first-principle calculations 
of the quat e rnary Heusler alloys Co2[Cri_^.Mnj.]Al, Co2Mn[Ali_xSnz] and [Fei_2:Coa;]2MnAl 
(jGalanakid . 1200J) demonstrated the Slater-Pauling behavior and half-metallic properties. 
Moreover, this behavior was investigated theoretically in V-based Heusler alloys Mn2VZ 
(Z=A1, Ga, In, S i, Ge , Sn) which are predicted to demonstrate half-m etallic ferrimagnetism 



(lOzdogan et al. 



2006 



Sasioglu et all l2005a : 



Weht and Pickett 



19991). 



An interesting feature of the half-metallic Heusler alloys is that the Rhodes- Wolfarth 
ratio Pc/Vs ipc are the effective moments, y^ th e saturation moments) can be essentially 
smaller than unity (jlrkhin and Katsnelsonl . Il994l ). Moreover, the effective moment in the 
paramagnetic state, determined from the paramagnetic susceptibility, decreases appreciably 
with temperature. We recall that for the Heisenberg magnets (atomic-like magnetic states) 
we have Pc/Ps = 1 and for weak itinerant ferromagnets such as ZrZn2 Vr.lv^^ ^ 1, i n 
accordance with the concept of "thermally induced local magnetic moments" (JMoriyal . ll985l ). 
In the case of conventional strong itinerant ferromag nets (for example, Fe and Ni) this ratio is 
also greater than unity (for detailed discussion, see (jLichtenstein et al\ . 120011 ) ) . It therefore 
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follows that the inequality pc/Vs < 1 is a striking property of HMF's, which could be used 
in their preliminary experimental identification. 

This behavior may be explained by that the change of electronic structure. The temper- 
ature dependence of magnetic moment in the paramagnetic state may be due to short-range 
magnetic order (local densities of states are similar to those in the ferromagnetic state). The 
numerical calculations demonstrate that the reduction in the moments is a consequence of 
a change in the electron structure as a results of rotation of the magnetic moments as it 



was demonstrated for Fe and Ni in Ref . ( Turzhevskii et al. 



19901 ). One would expect such 



changes to be particularly large in the case of HMF's and they should be of qualitative 
nature (smearing out of the hybridization gap because of spin disorder, for a review, see 
(jlrkhin and Katsnelsoru . Il994j )). From the many-electron model point of view, the decrease 
of the local moment with increasing temperature is connected with the absence of correc- 
tions to ground state magnetization of the type (188|) . However, such corrections do occur at 
high temperatures. 



G. Nuclear magnetic relaxation 

Nuclear magnetic resonance (NMR), which is one of most powerful tools for investigating 
various physical properties, has a number of peculiarities for magnetically ordered materi- 
als and especially for HMF. The localized-spin Heisenberg model is inadequate to describe 
the most interesting systems mentioned above where the role of conduction electrons is 
essential in magnetic properties. Usually the data on the longitudinal nuclear magnetic re- 
laxation rate 1/Ti are discussed within itiner ant-electron models such as Hubbard model or 



phenoraenological spin-fluctuation theories (jlshigaki and Moriya 



Moriyal . 



1985 



1994 



1996 



Millis et al. 



1990 



Ueda and Moriyal . Il975l ). On the other hand, in the s — d{f) exchange 



model (well-separated localized and itinerant subsystems) magnetic properties differ essen- 
tially from those in the paramagnon regime. We discuss the contributions to l/Ti(T) owing 
to electron-magnon interaction for three- and two-dimensional {3D and 2D) metallic fer- 



romagnets with well-de f ined 



(jlrkhin and Katsnelsoru . l200ll ). 



ocal magnetic moments with especial attention to HMF case 



The standard Hamiltonian of the hyperflne interaction, Tihf = hi, {ha = AajsSfs, A is the 
hyperfine interaction matrix) contains the Fermi (contact) and dipole-dipole contributions. 
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Aaf3 = A^5af3 + A^3 ■ According to (JAbragaml . Il96ll ) we have 



where 






(90) 
(91) 



F(°) = ((l-3cos2 0)/r3),F(^) = (sm0cos0exp(-20)/r^), 
F(2) = (sin2^exp(-2z0)/r3),a = --7e7n, 



(92) 



(...) is the average over the electron subsystem states, 7e and 7„ are gyromagnetic ratios 
for electron and nuclear moments, respectively. In the case of the local cubic symmetry we 
have F^""^ = 0. The Fermi hyperfine interaction is proportional to the electron density at 
the nucleus and therefore only s states participate in it, the contribution of core s states 
(which are polarized due to local magnetic moments) being much larger than of conduction 
electrons. It is just the consequence of considerably smaller localization area (and therefore 
higher density on nuclei) for the core states. It is obvious that magnetic f or d electrons 
dominate also in dipole interactions because of large spin polarization. Hence the direct 
interaction of nuclear spins with that of conduction electrons can be neglected in magnets 
with well-defined local magnetic moments. Nevertheless, conduction electrons do effect 
nuclear relaxation via their influence on the local-moment system; besides that, as we shall 
see below, such contributions possess large exchange enhancement factors. 

Using the e xpressions for 1/Ti and linewidth I/T2 in terms of the Green's functions 



flMorival . Il963h 



IT 

— = -— Im^((/i+|/i ))<^J^„, 

^1 ^^ V 



1 IT 

T2 2Ti 27r'^^o ^^^ 1' -'i^^'^' 

{un = {h^) ^ T is the NMR frequency) we derive 

Ti 2 3 

— = ^ + -{(A^--aTW)2ir^^ + a2[2|T«|2i^+-]} 
±2 2,1 1 2 3 

K-^ = -(1/vr) limlmY^{{S^\SZ^))Ju; 



(93) 
(94) 



(95) 
(96) 
(97) 
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Passing to the magnon representation we obtain 



(98) 



where Uf^ = 2S'( Jq — Jq) + uq is the magnon frequency, 7q(co') oc u; is the magnon damping. 
Then we have 

25 ,^ 7q(^n) 



K' 



71U„ 



E 



UJt 



q "q 

The damping owing to the one-magnon decay processes can be represented as 



7^^^^) 



27r/25^^5(tkT)5(ik-qi) 

k 



(99) 



(100) 



where the energy is referred to the Fermi level. The linearity of spin fluctuation damping in u 
is a characteristic property of metals. According to fl95|) this leads to T-linear contributions 
to 1/Ti which is the Korringa law. It is important that the simplest expression for the 
Korringa relaxation 

1/Ti ~ I/T2 ~ A^N^iEF)NiiEF)T, (101) 

{A is an effective hyperfine interaction constant) is practically never applicable for magnetic 
metals: exchange enha i icement facto r s cari change even the order of magnitude of l/Ti 



(jlrkhin and Katsnelson 



1994 



Moriya 



19851 ). Accurate expression for the "Korringa" con- 



tribution in the case under consideration can be derived by the substitution (1991) and (llOOp 
into (I95|) . 

Apart from three-dimensional case we can consider al so two-dimensional HMF keeping in 



mind, e.g., layered CMR compounds like LaSr2Mn207 (Ide Boer and Groot 



1999 



Nagaevl . 



20011 ). According to fl72l) . the damping (llOOp has the threshold value of g, which is determined 
by the spin splitting A = 2|/|S', q* = A/vp {vp is the electron velocity at the Fermi level), 
corresponding characteristic temperature and energy scale being u^ ~ {/S./vpY'Tc- After 
integration for the parabolic electron spectrum the one-magnon damping contribution to 
(!99|l takes the form 



5^^'^K' 






1/4, D = 3 

l/(7rg*), D = 2 



(102) 



with m the electron effective mas. Thus in the 3D case the factor of P is canceled, and the 
factor of I~^ occurs in the 2D case, so that we obtain a strongly enhanced T-linear Korringa- 
type term (remember that V r^ J r^ P/W for the RKKY interaction). This means that the 
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contribution of conduction electrons to T-linear relaxation rate via their interaction with 
localized spins is indeed much more important than the "direct" contribution: perturbation 
theory in the s — d exchange coupling parameter / turns out to be singular. 

In HMF the one-magnon decay processes are absent and electron-magnon (two-magnon) 
scattering processes should be considered (Sect JIII.E|) . Substituting the corresponding damp- 
ing into (199!) yields for D = 3 



5^^^K' 



finTi/2 



1287r' 



:wp^?^'^^^^ 



37ri/2C(|)r, T < 



UJ. 



h.2uj 



T>cu_ 



where C,{z) is the Riemann function, VLq the lattice cell volume. This contribution can 
also modify considerably the temperature dependence of 1/Ti in "usual" ferromagnets, a 
crossover from T^/^ to T^/^ dependence of the correction taking place. 

For D = 2 and T, u;_ ^ uq, small magnon momenta of order of {ujq/T)) ' make the main 
contribution to fl^ . On using the high-temperature expression A^p = T/tUp one gets 

^(2)^+- ^ 1,23- ^°^^ 



STT^^PS/S, 



uo, 



1/2 



T 



(103) 



Thus in the 2D FM case, in contrast with 3D one, the relaxation rate 1/Ti is strongly 
dependent on the anisotropy gap. It is worthwhile to note an important difference between 
relaxation processes via phonons and via magnons. The main difference is due to the gap 
in magnon spectrum. Usually ujq > Un and therefore one-magnon processes contribute 
to the relaxation rate due to magnon darn ping only (cf. discussion of the phonon-induced 
relaxation processes in Ref . (I Abragaml . 1 1 96 ll ) ) . However, the mechanisms of magnon damping 
in magnetic dielectrics (magnon-magno n interactions) are different from tho se in magnetic 



metals and degenerate semiconductors (jAuslender and Irkhin 



1984a 



1985af ) 



H. Thermodynamic properties 

Consider the renormalization of electronic specific heat in an itinerant ferromagnet due 
to interaction with spin fiuctuations. Integration in (JHl), (j9]) at T = gives 



ReY:^{kFa,E) 



— ^^^iV_(E^) 5: a(E - u;„) In ^^-^ 



a=± 



Then the inverse residue of the electron Green's function, 1/Z]^„{E) = 1 — (9/9i?)ReSko-(i?), 
which determines the renormalization of the electron effective mass owing to the electron- 
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magnon interaction, contains a logarithmic factor. We obtain for the coefficient at the hnear 
term in the electronic specific heat at T ^ cj. 



7, = T<»'/Z,(*Fa, Ef) = yJV.lBF) 



"^ -N_,(Ef)U.''- 



UJ+ — UJ_ Co'_ 



(104) 



For weak itinerant ferromagnets we have 



In — ~ -2 \n{UN{EF) - 1) (105) 

SO that the we have a paramagnon enhancement of the specific heat, the numerical factor 
i n (11041) being inexact in this limit because of neglecting longitudinal spin fluctuations (see 



(JMoriyal . Il985l )). On the other hand, we have a considerable enhancement of specific heat 
owing to spin fiuctuations strong ferromagnets which is really observed in a number of 
systems. 

Other thermodynamic properties may be treated by calculating the free energy of the 
system. The spin-wave contribution to the free energy has the form, usual for the Bose 
excitations with the square dispersion law and is proportional to (T/TcY^'^. At low T < co- 
the many-electron (branch cut) contribution reads 

fe, = ^V £ '^'■<«-'.«J> - ^^ £ , ""f 7'"' (106) 

^>J0 q>q. kk' I'T ~ '^k'J. + t^k-k' 

Differentiating ( 1106^ over T one obtains 

= U'^^^^N,{E,)N,{E,f-^T\n f±-— (107) 

uJ-^- — uj- 6 max(u;_, I ) 

Thus at T ^ c<j_ we have instead of (I104p the T In T-dependence of specific heat. 

As one can see from (11041) . the enhancement of effective mass and electronic specific heat 
owing to spin fiuctuations is absent in the half-metallic state. We shall demonstrate that 
the specific heat of a conducting ferromagnet may contain spin-fiuctuation contributions of 
another nature. Write down a general expression for the specific heat in the s — d exchange 
model in terms of the total energy 

C{T) = ^-^ = ^jdEEf{E)N,{E) 

= ^Nt{E)T + J dEEf{E)-^Nt{E,T) (108) 
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where 



Nt{E) 



TX 



Y^\x^G^,,{E) 



kcr 



The first term in the right-hand side of fllOSp yields the standard result of the Fermi- 
liquid theory. The second term is due to the energy dependence of the density of states. 
Such a dependence occurs in the conducting ferromagnet owing to non-quasiparticle (in- 
coherent) states. Using aga in the expressions for the self-energies dH]), © we derive 
(jlrkhin and Katsnelsonl . 119901 ) 






^k+q,-o- 



(109) 



Since at low temperatures /(tk+q,i~^q) = I5 /(^k+q,t~^q) = 0; the non-quasiparticle states 
with a =1 do not contribute to linear specific heat since they are empty at T=0 K. In the 
half-metallic state the non-quasiparticle contributions fllOSp with o =\ are present for J < 
only, and we obtain 



2tx 



SC,{T) = —I'{S^)N,{E^)TJ:^ 



1 



E^Y 



;iio) 



k (^kT 

To avoid misunderstanding, it should be stressed that presence of such contributions to 
specific heat means inapplicability of the Fermi-liquid description in terms of dynamical 
quasiparticles only, which are determined by poles of Green's functions. It may be shown 
rigorously that the entropy of interacting Fermi systems at low T is expressed in terms of 
Landau quasiparticles with the energies, determined as variational derivatives of the total 



energy with respect to occupation numbers (jCarneiro and Pethickl . Il975l ). Thus, even in the 



presence of non-pole contributions to the Green's functions, the descrip tion of thermodynam 



ics in terms of statistical quasiparticles (iCarneiro and Pethickl . Il975l ) holds. (However, the 
quasiparticle description is insufficient for spectral characteristics, e.g., optical and emission 
data.) The anomalous 7T-term is determined by the difference of the spectra of statistical 
and dynamical quasiparticles. 

Similar contributions to specific heat in the Hubbard model with strong correlations 
are discussed in the paper (jlrkhin and Katsnelsonl . Il990t ) too. They dominate in the en- 
hancement of specific heat for half- metallic ferromagnets and may be important, besides the 
effective mass enhancement (11041) . for "usual" magnets with well-defined local moments. 
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I. Transport properties 



Transport properties of HMF are the sub j ect o f numerous exp erimental investig ations 



see, e.g., re cent works for CrOo (JRabe et al 



the reviews ( llrkhin and Katsnelson 



1994; 



2002 



Nagaev 



and NiMnSb ( 



2001 



Ziese 



Borca et a/. , 



2001 



), and 
2OO2I )). At the same time, 



the theoretical interpretation of these results is still a problem. As for electronic scattering 
mechanisms, the most important difference between HMF and "standard" itinerant electron 



ferromagnet s like iron or nickel is t 
former case (llrkhin and Katsnelsonl . 



le ab sence of one-magnon scattering processes in the 
199J). 



Since the states with one spin projection only exist at the Fermi level and one-magnon 
scattering processes are forbidden in the whole spin-wave region, the corresponding T^- 
term in resistivity is absent in the case of a half-metallic ferromagnets. This seems to be 
confirmed by comparing experi mental data on re sistivity of Heusler alloys TMnSb (T = Ni, 
Co, Pt, Cu, Au) and PtMnSn (lOtto et all Il989l ) (see also discussion in Sect. IV. A. II) . The 
T^-contribution from one-magnon processes to resistivity for half-metallic systems (T = Ni, 
Co, Pt) was really not picked out, whereas the dependences piT) for "usual" ferromagnets 
were considerably steeper. 

Two-magnon scattering processes have bee n considered many years ago, the t empe rature 



dependence of resistivity obtained being T^/^ f lHartman-Boutron 



1965 



Roeslerl . 



19651). The 



obtained temperature dependence of the resistivity has the fo rm T^/^. At low enough ten i- 



peratures the first result fails and should be replaced by T^'^ ( jLutovinov and Reizei 



19791); 



the reason is the compensation of the transverse and longitudinal contributions in the long- 



wavelength li mit which is a consequ e nce of the rotationa 



Hamiltonian (jAuslender et al. 



1983 : 



Grigin and Nagaev . 



symmetry of the s — d exchange 



1974 



Nagaev 



19831). We discuss 



effects of interaction of current carriers with local moments are investigated in the standard 
s — d exchange model in the spin-wave region. 



n = 7^o-/(25)'/'E(CkTCk+qi&; + /^.c.) 

kq 



(111) 



kqpcr 
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The zero-order Hamiltonian includes non-interacting electrons and magnons, 

Ho = Yl -^kaC^Cka + Yl ^ci^lP' 



q"q' 



(112) 



kcr 



^ka = tk - a A/2, cuq = 2S{Jo - Jq), 



with A = 21 S being the spin splitting which is included in ?io- In the half-metallic case 
the spin-flip processes do not work in the second order in / since the states with one spin 
projection only are present at the Fermi level. At the same time, we have to consider the 
renormalization of the longitudinal processes in higher orders in / (formally, we have to 
include the terms up to the second order in the quasiclassical small parameter 1/5). To 
this end we can eliminate from the Hamiltonian the terms which are linear in the magnq n 



operators by using the canonical transformation (JGrigin and Nagaev 
Then we obtain the effective Hamiltonian 



1974 



Nagaevl . 



19831). 



1 



^ = ^0+TT E(^kq + ^k- 



-q-p,q/'^kCT'^k+q-pCTOqOp 



(113) 



kqp cr 



Here ^^q is the s — d scattering amplitude which is deflned by (TT2|) {U ^ I). More general 
interpolation expression for the e ffective amplitude, whi ch does not assume the smallness of 



|/| or 1/25 was obtained in Ref . (jAuslender et al. 



19831 ) by a variational approach. 



The most general and rigorous method for calculating t he tra n sport relaxation time 
is t he use of the Kubo formula for the conductivity a^x (IKubol . 119571 ) (see details in 



Ref. dlrkhin and Katsnelson 



2002)): 



/•/3 roo 

crxx = P dX dtexp{-£t){j^{t + iX)jx) 
Jo Jo 



where (3 = 1/T, e ^ 0, 



J = -e X! VkacLcka 

ko- 



(114) 



(115) 



is the current operator, Vko- = dE^^/dk is the electron velocity. Representing the total 
Hamiltonian in t he fo r m 7i = Tip + T^^^the correlator in (I114p may be expanded in the 



perturbation 7i' (JMoril . 

resistivity 



1965 



Nakano 



3 



19571). In the second order we obtain for the electrical 



Pxx = 0"i 



Y POO 



{Px? Jo 

where ?i'(t) is calculated with the Hamiltonian ?io- 



dt{[jx,n'mn',jx]) 



(116) 
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In the HFM situation the band states wit 



at the Fermi level (jlrkhin and Katsnelson 



1 one spin projection only, a = sign/, are present 



199J). We consider the case I > 0, a = + and 



omit the spin indices in the electron spectrum to obtain for the transport relaxation time r 
defined by a^x = e^((f^')^)r 

1 

r 



TV 

4T 



E (^k - <')'(< + 4',q-lc'+k)'^q(l + ^q-k'+k) 



kk'q 



xnk(l -^k') 

X5(tk' - tk - ^q + ^q-k'+k) /E«)^'^(^k) 

Averaging over the angles of the vector k leads to the result 1/r oc /^$ with 
where /pq = 1 for p, g :^ q* and 



;ii7) 



(118) 



/, 



[p X q]' 



pq 



(p-q) (g*)2 



(p, g < g* 



(119) 



As discussed before, the wavevector q* determines the boundary of the region where q- 
dependence of the amplitude become important, so that t(k + q) — t(k) ~ A at g ~ g*. In 
the case q < q* the simple perturbation theory fails and we have to take into account the 
spin splitting by careful collecting the terms of higher orders in /. In the simple one-band 
model of HMF where Ep < A one has q* ~ JA/W {W is the conduction bandwidth) 



(JGrigin and Nagaevl . Il974j : iNagaevl . Il983l ). The quantity g* determines a characteristic tem- 
perature and energy scale T* = V{q*Y oc V{/S./W). 

When estimating temperature dependences of resistivity one has to bear in mind that 
each power of p or g gives T^/^. At very low temperatures T < T* small quasimomenta 
p,q < q* yield main contribution to the integrals and 



p{T) oc {T/T, 



c 



,9/2 



(120) 



Such a dependence was obtained in the narrow-band case (d ouble-exchange model with large 



|/|), where the scale T* is absent (IKubo and O 



in the broad-band case (jLutovinov and Reizer 



lata 



m 



19721 ). and by the diagram approach 



19791). At the same time, at T > T* the 



function /pq in Eq. (11181) can be replaced by unity to obtain 

p{T) oc {T/Tcy" 



;i21) 
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Generally speaking, q* may be sufficiently small provided that the energy gap is much 
smaller than W , which is the case for real HMF systems. We consider the model of HMF 
spectrum fl50p where the majority-spin band is metallic and the minority-spin is semicon- 
ducting. For the temperatures T <^Tc both characteristic q and p are small in comparison 
with the inverse lattice constant and we can put 



AZ^ 



1 



{tki - ^kt) 



k+qT - ^kT 



(122) 



q 25* ' "■" "" tk+qt — tki 

The wavevector q* determines the boundary of the region where q-dependence of the am- 
plitude become important, so that q* = A/vp (the same value as for the spin polarization 
problem). The corresponding characteristic temperature and energy scale is 



T* =V{q*y ocV{A/Wy 



(123) 



which coincides with the case of an usual ferromagnetic metal. The above temperature 
dependences of resistivity are not changed (jlrkhin and Katsnelsonl . 120021 ). 

Now we t reat the two-dimens i onal ( 2D) situation which may be appropriate for layered 



manganites (Ide Boer and Groot 



1999 



Nagaevl . l200ll ). At low temperatures we obtain 



p{T < T*) oc {T/Tc 



c 



>7/2 



(124) 



At the same time, for T > T* we obtain after replacing the scattering amplitude by unity a 
logarithmically divergent integral which should be cut at T*. Thus we get 



p{T > T*) oc {TlTcf \n{T/T*) 



(125) 



To discuss the magnetoresistivity we have to introduce the gap in the magnon spectrum, 
cijq^o = ^Q^ + ^0- Provided that the external magnetic field H is large in comparison with 
the anisotropy gap, ujq is proportional to H. In the 3D case the resistivity at T < T* is 
linear in magnetic field, 

p(T, H) - p{T, 0) oc -uoT'/yr'J' (126) 

The situation at T > T* is more interesting since the quantity d^/dujQ contains a divergence 
which is cut at c^o or T*. We have at T > ujo,T* 



Sp{T, H) oc 



T^o^o 



[max(o;o,T*)]i/2 



(127) 
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(of course, a.tT < ujq the resistivity is ex ponentially s mall) . A negative i7- linear magnetore- 
sistance was observed recently in Cr02 ( iRabe et all |2002| ). 

The discovery of giant magnetoresistance (GMR) lead to tremendou s activity to under- 
stand and develop technology based on high-density magnetic recording fIZutic et al\ . \2004 ). 
Giant magne toresistance for metalli c mult ilayers (superlattices) containing HMF was first 
predicted in (llrkhin and Katsnelsonl . [l99J) . NiMnSb-based spin- valve str uctures using Mo 



space r layers NiMnSb/Mo/NiMnSb/SmCo2 were successfully produced (JHordequin et al. 



9981). The associated GMR exhibits a clear spin- valve contribution about AR/R ~ 1% 



(JHordequin et all Il998l ) . One of the limiting factor for such a small value is the large resis- 



tivity of the Mo layer which determines limited flow of active electrons exchanged between 
the two ferromagnetic layers without being scattered. The giant tunneling magnetoresis- 
tance differs by use of dielectric spacer instead of met allic one. The GMR in t unnel junctions 



based on HMF was considered theoretically in Ref. (Ide Groot et al. 



1983a 



Tkachov et al. 



200lh and recently this issue becanie a subject of intensiye exp e rimental inv estigations 



(jGercsi et al. 



2006 



Rybchenko et al. 



2006: ISakuraba et al. 



2006a 



2007 



2006bh 



It is important that NQP states do not contribute to the temperature dependence of th e 
resistivity for pure HMF. An opposite conclusion was made by Furukawa (JFurukawal . l2000l ) . 
He attempted to calculate low-temperature resistivity of half-metals with account of the 
non-rigid-band behavior of the minority band due to spin fluctuations at finite temperatures 
and derived that the unconventional one-magnon scattering process give T^-dependence in 
resistivity. However, this calculation was not based on a consistent use of the Kubo formula 
and, in our opinion, can be hardly justified. 

On the contrary, impurity contributions to transport pro perties in the presence of po- 



tential scatt e ring are det ermined mainly by the NQP states (llrkhin and Katsnelson 



Irkhin et al. 



1989 



19M)). 



1994; 



To second order in the impurity potential U we derive after 



neglecting vertex corrections and averaging over impurities we obtain for the transport re- 
laxation time 

'5<(^) = -2W^ImEGiS(i?) (128) 

p 

where is the exact Green's function for the ideal crystal. Thus the contributions under 
consideration are determined by the energy dependence of the density of states N{E) for 
the interacting system near the Fermi level. The most nontrivial dependence comes from 
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the non-quasiparticle (incoherent) states with the spin projection —a = —sign/, which are 
present near Ep. Near the Fermi level the NQP contribution is determined by the magnon 
density of states g{uj) and follows a power law, 

rcrE 

6Ni^^oh{E)oc dujg{u)(x\E\''e{aE) [lEl-^uJ). (129) 

Jo 

Here 6{x) is the step function, E is referred to Ep; we have a = 3/2 and a = 1 for 3D and 
2D cases, respectively. The corresponding correction to resistivity reads 

^^^ = -Sa,^AT) (130) 

P 

oc -U'JdE (-^^] mncoi.{E) oc T" 
The contribution of the order of T° with a ~ 1.65 (which is not too far from 3/2) was 



observed in the temperature dependence of the resistivity for NiMnSb (JBorca et all l200ll ) 
above lOOK. The half-metallic properties above lOOK are being challanged, however. The 
incoherent contribution to magnetoresistivity is given by 

(5pi^p(T, H) (X uJodSN,^,oi.{(TT)/dT ex cuoT"-\ (131) 

so that we obtain a temperature-independent term in the 2D case. 

le correction to thermoelectric power, which is similar to (I130p . reads (cf. 



T 



Ref . (Ilrkhin et all Il989l . Il994h ) : 



SQ(T) cx^ f dE{-df{E)/dE)E5N{E) (132) 

Besides that, an account of higher orders in impurity scattering leads to the replacement of 
the impurity potential V by the T-matrix. For the point-like scattering the latter quantity 
is given by 

Expanding 0133p yields also the term 

5Q{T) oc ^ f dE{-df{E)/dE)ERe6n{E) (134) 

with 6TZ{E) being obtained by analytical continuation from 6N{E). Thus we have 6Q{T) oc 

2^3/2 _ 
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J. X-ray absorption and emission spectra. Resonant x-ray scattering 

The NQP contributions in the presence of the potential U, that is induced by the impurity 
at a lattice site, can be treated in the s — d exchange model in the representation fl23l) . The 
impurity potential results in the NQP contribution to this quantity being enhanced for U < 
and suppressed for U > 0. These res ults can be used to cons i der the manifestations of NQP 
states in the core level spectroscopy (llrkhin and Katsnelsonl . l2005al ). 

Various spectroscopy techniques such as x-ray absorption, x-ray emission, and photo- 
electron spectroscopies (xas, xes, and xps, correspondingly) give an important information 
about the electronic structure of HMF and related compounds, i.e. f erromagnetic semicon 



ductors and colossa l magnetoresistance materials (see, e.g., Refs . ([Kurmaev et al. 



Wesselv et al 



2003 



Yablonskikh et al. 



2001 



Yarmoshenko et al. 



2003 : 



1998D). It is well known 



( jMahanl . Il990l ) that many-body effects (e.g., dynamical core hole screening) can be impor- 
tant for the core level spectroscopy even when the system is not strongly correlated in the 
initial state. Therefore it is very interesting to study the interplay of these effects and NQP 
states which are of essentially many-body origin themselves. 

To consider the core level problem in HMF we use the Hamiltonian of the s — d exchange 
model in the presence of the external potential U induced by the core hole. 



kk'cr 



(135) 



where /''', / are core hole operators, U < 0. X-r ay absorption and emission spectra are 
determined by the Green's function (JMahanl . Il990l ) 



Gl^,iE) = {{cU\fcl^))E 



U 



(136) 



As fo: 



lows the investigation of the sequence of equations of motion (llrkhin and Katsnelsonl . 
2005al ). in the ladder approximation the spectral density for two-particle Green's function 
^kk'(-^) ^^ equivalent to the one-particle spectral density in the presence of the core hole 
potential U (note that the ladder approximat ion is inadequ ate to describe the xas edge 
singularity in a close vicinity of the Fermi level (jMahanl . Il990l )). Thus the core hole problem 
is intimately connected with the impurity problem. 
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Since xas probes empty states and xes occupied states, the local DOS 

1. 



KciE) 



TX 



-lu^GUE) 



(137) 



describes the absorption spectrum for E > Ep and emission spectrum for E < Ep. 
To take into account the core level broadening a finite damping 5 should be introduced 
jlrkhin and Katsnelsonl . 



2005al ). For small band filling the "exciton effects" (strong interac- 



tion with the core hole) result in a considerable enhancement of NQP contributions to the 
spectra in comparison with those to DOS. The results for semielliptic bare band are shown 
in FigsHma 

To probe the "spin-polaron" nature of the NQP states more explicitly, it would be desir- 
able to use spin-resolved spectro scopica l meth ods such as x-ray magnetic circular dichroism 
(XMCD, for a review see Ref. ( Ebertl . Il996l )). Owing to interference of electron-magnon 
scattering and "exciton" effects (interaction of electrons with the core hole), the NQP con- 
tributions to x-ray spectra can be considerably enhanced in comparison with those to DOS 
of the ideal crystal. Thus the core level (x-ray absorption, emission and photoelectron) 
spectroscopy might be an efficient tool to investigate the NQP states in electron energy 
spectrum. 

No w we consider NOP eff ects in resonant x-ray scattering processes. It was observed re- 
cently (JKurmaev et a/.l . l2003l ) that the elastic peak of the x-ray scattering in Cr02 is observed 
which is more pronounced than in usual Cr compounds, e.g., in elemental chromium. The 
authors of this work have put forward some qualitative arguments that the NQP states may 
give larger contributions to resonant x-ray scattering than usual itinerant electron states. 
Here we can treat this question quantitatively and estimate explicitly the corresponding en- 
hancement. The intensity of resonant x-ray emission induced by the photon with the energy 
ui and polarization q is given by the Kramers-Heisenberg formula 

(n|Q|0(/|C,|0) 



Iq,q{uj\uj) OC ^ 



y 

^ Eo + uj'-Ei 



iTi 



6{En + uj' -Eo 



00 



(138) 



Here q' and u' are the polarization and energy of the emitted photon, \n), |0) and |/) are 
the final, initial and intermediate states of the scattering system, respectively, Ei are the 
corresponding energies, Cg is the operator of the dipole moment for the transition, which is 
proportional to fc + c^f'^. Assuming for simplicity that P; does not depend on the interme- 
diate state, P; = P, and taking into account only the ma in x-ray scattering c hannel (where 
the hole is ffiled from the conduction band) one obtains (ISokolov et ali 119771 ) 
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luj' oc 



E^ooW 



(139) 



where z = u' ~ Eq + iT. Owing to a jump in the DOS at the Fermi level, the NQP part 
of the Green's function contains a large logarithm \n{W/z) at small z. It means that the 
corresponding contribution to the elastic x-ray scattering intensity {u' = Eq) is enhanced 
by a factor of In^ fH^/F), which makes a quantitative estimation for the qualitative effect 



discussed in Ref. ( JKurmaev et a/.l . l2003l ). Of course, the smearing of the jump in the density 



of NQP states by spin dynamics is irrelevant provided that F ^ u; (cu is a characteristic 
magnon frequency). 

IV. MODERN FIRST-PRINCIPLE CALCULATIONS 

A. Different functional schemes 



In this section we review contemporary approaches to the electronic structure calcula- 
tions with taking into account correlation effects. Model considerations discussed above 
demonstrate relevance of the correlation effects (such as electron-magnon interactions) for 
physics of half-metallic ferromagnets. In order to calculate the electronic structure of real 
materials we have to solve a complicated many-body problem for a crystal, corresponding 
to inhomogeneous gas of interacting electrons in an external periodic potential: 



7i = Tif) + 7i 



int) 



Ho = E/rfr'0.^(r)[-^V2+Kxt(r)]^.(r), 
Hint = iE//^rrfrV+(r)^+(r')\/(r-r')^.'(r')^.(r). 



(140) 



In this section we use the atomic units {h = m = e = 1), ip^j (r) is a field operators for 
electrons, V"cxt(r) describes interaction of electrons with static nuclei which are supposed 
to form the periodic crystal lattice and also may include other external potentials (defects, 
electric fields, etc.) and V{r — r') =l/|r — r'| is the Coulumb interaction between electrons. 
The modern view on various practical schemes for solution of this general many-electron 
problem is based on its f unctional formulation in a framework of s o -calle d effective action 



approach (JFukuda et al 



1994 



Georges et al 



19961: 



Kotliar et al. 



20061 ). The partition 
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function of electronic system within imaginary-time functional integral formalism can be 
expressed as an integral over electronic Grassmann variables: 

a- •' 

S = J2 1 dr iyr^^ir,T)^M^,T)+ jyrnir). (141) 

The free energy of many-electron system F = —T\nZ can be expressed as a function 
of optimally chosen physical variables for a given problem. The most accurate scheme 



corresponds to Baym-Kad anoff, or Luttinger-Ward functional (jBaym and KadanofI 



Luttinger and Ward 



d) 



1961 



19601) of the one-electron Green's function: 



G,(r - r', r - r') = - (T.^.(r, r)V^+(r',r')) (142) 

In this formulation one has to add constraint fields of dual variable S to the action 



r/3 rP 



Sm = S + J2jdrJdr' J^ drj^ rfrU(r, r', r, r')a(r', r, r', r) (143) 

and find the partition function in the presence of the auxiliary source field: 

Z[E] = e-^Pl = f D[tlj+tlj]e-^^''^ (144) 

The corresponding Baym-Kadanoff functional is defined as the Legendre transformation of 
F[S] to the Green's function variable: 

F[G] = F[S] - Tr(SG') (145) 

with further use of the functional derivative G = 6F/6T, to eliminate the constraint fields. 
Using free-electron Green's function corresponding to Tio part of the Hamiltonian, the final 
form of the functional with "Kohn-Sham" decomposition can be written in the following 
form: 

F[G] = -Tr ln(Go^ - S) - Tr(SG') + <^[G]; (146) 

here $[G] is the Luttinger generating functional which can be represented as a sum of all 
irreducible diagrams without legs constructed from the exact Green's function G and bare 
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electron-electron interaction line (bare four-leg vertex) V. The Baym-Kadanoff functional is 
stationary in both G and S and its variation with respect to S leads to the Dyson equation 



G— 1 /^^~i 
— ^0 



(147) 



and the G-extremum gives the variational identity, S = 6^/6G. 

The Baym-Kadanoff functional allows us in principle to calculate not only free energy 
and thus thermodynamic properties of the system, but also the Green's function and thus 
the corresponding excitation spectrum. The main point which makes this scheme rather 
useful for model many-body analysis and preserves its broad practical use in the electronic 
structure calculation is related with difficulties to find an exact representation of F[G] even 
for simple systems. 

In this situations the density-functional schem e of Kohn, Hohenberg, and Sham 



( iHohenberg and Kohn 



196 



I 



Kohn and Shaml . Il965l ) turns out to be the most successful 



scheme for the electronic structure calculations of an electronic systems with not too strong 
correlations. For this purpose, the functional of the static electronic density 



Pa{r) = lmGa{r,T = 0). 



148) 



is constructed. The corresponding constraint fields in the effective action are related to the 
Kohn-Sham interaction potential which is represented as a sum of Hartree and exchange- 
correlation (xc) parts: 

l^nt = Vh + Kc (149) 

Finally, the Kohn-Sham free-energy functional can be written in the following form: 



F[p] = -TrlniG^' - V^^t - l^mt) - Tr(Mntp) + Fh[p] + F,,[p]. 



(150) 



where Kxt is an external potential and Fh is the Hartree potential. Again, there is a similar 
problem: an exact form of the exchange-correlation functional F^^p] is, generally speaking, 
unkno wn, and only a formal ex pression in terms of the integral over the coupling constant 
exists ( JHarris and Jonesl . Il974l ). The practical use of the density-functional theory (DFT) 
is related with the local-density approximation (LDA): 
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F,,[p] ^ J drp(r)e,e[p(r)] (151) 

where exc[p] is the exchange-correlation energy per particle of homogeneous electron gas 
with a given densit y. It can be careful l y par ametrized from the numerically exact Monte- 
Carlo calculations (jCeperley and Aided . Il980l ). Taking into account spin- dependence of the 
DFT through Po-(r) one can study magnetic properties of complex materials. This was the 
method used in the most of electronic structure calculations referred above; in particular, 



the concept itself of the hal 



calculations (Ide Groot et al 



-metal 



ic ferromagnetism was introduced based on this kind of 
1983bl ). In practice, the LDA scheme results sometimes in well- 
known difficulties; in particular, it underestimates usually energy gaps in semiconductors. 
For this reason it can fail to describe pr operly the half-me t allic state, e.g., in the case 
of colossal magnetoresistance manganites (jPickett and Singhl . Il996l ). The DFT scheme is 
formally exact (assuming that an exact E^c is known) to find the energy and electronic 
density of the many-body systems by minimization of density functional. However, the 
excitation spectrum, rigorously speaking, cannot be expressed in terms of the Kohn-Sham 
eigenenergies Si defined by 



;-2 V^ + Kxt + ^nt)^i = Silpi 



(152) 



(see, e.g., the discussion of NQP contributions to thermodynamic properties in Sect JIII.Hl) . 

A reasonable scheme which can overcome the difficulties of the DFT scheme for the 

gap problem uses so-called GW approximation proposed by Hedin ( 



functional approach to the GW scheme has been developed recently (JAlmbladh et a/.l . ll999 : 



Chitra and Kotliar 



Hedin 



1965 L The 



2000 



20011 ) and is related with the free energy functional of both total 
Green's function G and screened Coulomb interactions W = {V^^ — H)~^: 



F[G,W] = -Trln(G'o^-S)-Tr(SG') + -Trln(\/-i-H)+-Tr(HIV)+FH[p]+$[G,W^] (153) 

where H is the polarization operator. Earlier a similar ap proach was used in the theory o f 
phonon-induced superconductivity of disordered systems (JAnokhin and Katsnelsonl . Il996l ). 
In the GW approximation only the lowest-order diagram in the screened interactions is 
included into the generating functional: 
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<^[G,W] = -Tt{GWG) 



(154) 



In this case the polarization operator 11 which serves as a constraint field for the screened 
Coulomb interactions W has the simplest form: 



n 



' 6W 



-GG 



(155) 



and the corresponding electron self-energy reads 



oG 



(156) 



The GW scheme gives an accurate estimation of the screened Coulomb interactions in solids 
and can be used to define the first-principle values of local Hubbard-like multi-orbital energy- 



dependent i nteractions for corre l ated 



metal ions (JAryasetiawan et all 1200J): 



ocal orbitals dAr) which describe d states of transition 



Uijki{uj) = ((l)i(l)j\W{uj)\(j)k(j)i 



(157) 



where W does not take into account the effects oi d — d screening, the latter being ex- 
plicitly taken into account further within an ef fective low-energy Hubbard -like model. The 



numerical estimation of U for metallic nickel JAryasetiawan et al. 



2004 ) shows relatively 



weak energy dependence within the d-band energy width and the static values of the order 
of 2-4 eV, in a good agreement wit h the experimental values of the Hubbard parameters 
( Ivan der Marel and Sawatzkyl . Il988l ). 

The success of GW approximation is closely related to the fact that the bare Coulomb 
interaction V is strongly screened in solids and thus one can use the lowest-order approxi- 
mation for $[G, W]. On the other hand, the spin dependence of self-energy in GW scheme 
comes only from the spin dependence of the Green's function Ga{Y^T) and not from the 
effective interactions W . In the Baym-Kadanoff formalism this means that only the density- 
density channel was taken into account in the screening of the Coulomb interactions. It 
is well k nown that the Hund's intraatomic exchange interactions are weakly screened in 
crystals (Ivan der Marel and Sawatzkyl . Il988l ). and strong spin-flip excitation processes will 
modify electronic self-energy in itinerant electron magnets. In particular, these processes 
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are responsible for electron-magnon interaction and lead to appearance of the NQP states 
in the gap region for half-metallic ferromagnets. 

An accurate treatment of the effects of local screened Coulomb U and exchange J inter- 
actions beyond GW or DFT methods can be carried out within the dynamical mean-field 
theory (DMFT) combined with the GW or LDA/GGA functional. The DMFT scheme 
defines the best local approximation for the self-energy, which uses the mapping of the origi- 
nal many-body system with Hubbard-like interactions onto multi-orbital qua ntum impurity 



mode l in the effective electronic bath under the self-consistency condi tion (JGeorges et al. 



19961 ). The corresp onding GW-I-DMFT scheme (JBiermann et alV l2003l ) or spectral-density 
functional theory (jSavrasov and Kotlian . l2004j ) is probably the best known way to treat 
correlation effects in the electronic structure of real materials. However, it is still very 
cumbersome and computationally expensive; also, methods of work with the frequency- 
dependent effe ctive interaction (11571) are not develope d enough yet (for examples of first 



attempts, see (JRubtsov et al. . 



2005 



Savkin et al. 



20051 )). The only way to consider effects 



of spin-flip processes on the electronic structure of real materials is a sim plified version of 



a general spectral-density functiona 



1997bl : 



Lichtenstein and Katsnelson 



mown as the LDA+DMFT approach (JAnisimov et al\ 



19981 ). In a sense, one can consider the LDA+DMFT 



and LDA as complementary approaches. In both the cases we split a complicated many- 
body problem for a crystal into a one-body problem for the crystal and many-body problem 
for some appropriate auxiliary system where we can hope to calculate the correlation effects 
more or less accurately. For LDA we choose the homogeneous electron gas as this auxiliary 
system. For the LDA+DMFT it is an atom in some effective medium. The latter choice 
is optimal to consider atomic-like features of d or f electrons in solids. The local Green's 
function in magnetic solids is obtained from the effective impurity action with the static 
(frequency-independent) Hubbard-like multi-orbital interactions. 



'S'imp = -Y. ^^ dT'cf{T)gij{T - t')cj{t') + - JI / c?rc+(r)c+(r)f/ijfc/Cfe(r)Q(r), 

(158) 
where Ci{T) are the fermionic Grassmann variables for localized correlated ci-orbitals 0i(r) 
and Qij is so-called bath Green's function which is defined self-consistently within the single- 
particle lattice model. The corresponding interacting local Green's function 
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G,,{T-T') = -{TMrH{r')) 



Sin 



(159) 



can be found, within numerically exact Quantum Monte-Carlo scheme ( JHirschl . Il983l ). or 
some perturbative ap proach which treats accurate l y spi n-flip excitation processes in the 
particle- hole channel (JKatsnelson and Lichtensteinl . Il999h . The corresponding self-energy 
matrix of the impurity model 



can be used in the spectral-density functional: 



(160) 



F[g] = -Tr ln(Go' - S) - rr(S^) + <^[g] 
and satisfies self-consistent equation for the bath Green's function: 



(161) 



Giuj) = ^[G,\k,uj) -j:{^)]-' + ^iuj) 



(162) 



B. LDA-|-DMFT: the Quantum Monte Carlo solution of the impurity problem 



Now we describe the most rigorous way to solve a n effective i r apurit y problem using the 



multi-band Quantum Monte Carlo (QMC) method 



f|Rozenberd. 



1997h. In the 



ramework 



of LDA+DMFT this approach was used first in Ref . (JKatsnelson and Lichtenstein 
the case of ferromagnetic iron. 



2000) for 



We start from the many-body Hamiltonian in the LDA+U form (jAnisimov et a/.l . ll997al ): 



'T_/dc 
fl-LDA 



n = nt^^ + l y: k 

i{a-m} 
ijcr{m} 



mim2m'^m'2 imia im2cr' 'i''n^2"' 



2'j ^im'^c 



(163) 



where {ij) represents different crystal sites, {m\ label different orbitals, {cr} are spin indices, 
and t^^m2 ^^^ ^^^ hopping parameters. The Coulomb matrix elements are defined by 

UL,m2m',..'2 = j j dvdv'mi^^{v)%^y)V,,{v - r')vi>,^; (r) v]/,^. (r') (164) 
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where Vpp(r — r') is the screened Coulomb interaction which remains to be determined. 



We follow again the spirit of the LDA+f/ approach by assuming that within the atomic 
spheres these interactions retain to a large measure their atomic nature. Moreover, the 
values of screened Coulomb (?7) and exchange (J) inte r action s can be calculated within the 
supercell LSDA approach (jAnisimov and Gunnarssoru . Il99ll ): the elements of the density 
matrix n^^/ are to be constrained locally, and the second derivative of the LSDA energy 
with respect to the variation of the density matrix yields the wanted interactions. In a 
spherical approximat ion, the ma trix elements if V^g can be expressed in terms effective 
Slater integrals F^''^ ( I.Tuddl . \l9M) as 



(m, m"\Vee\m', m'") = ^ ak{m, m', m", ■m"')F^^\ 

k 



(165) 



where < k < 21 and 

ak{m, m' , m" , m'") 



Att 



2k + 



- 5:(/m|n,|/m')(/m"|n;|/m'") 



For d electrons one needs F^^\ F^^^ and F'^'^^; they are connected with the Coulomb- and 
Stoner parameters U and J hj U = F(°) and J = (F^^) + F('^))/14, w hile the ratio F^ ^ VF^"^) 



is to a good accuracy a constant, about 0.625 for the 3c? elements ( jAnisimov et al 



de Groot et al 



1993 



19901 ). T^triA represents the LDA Hamiltonian corrected by double counting 



of ave rage static Coulomb interaction that is already presented in LDA (jAnisimov et al . 
1997al ). The index i for the W has a meaning only for the same correlated sites as the 
orbital indices {m}, unlike the LDA term K^^^^ (one-particle Hamiltonian parameters) 
where we have the contribution of all the sites and orbitals in the unit cell. 

The one-particle spin-polarized LDA+DMFT Green's function G(j(k, cj) is related to the 
LDA Green's function and the local self-energy So-(co') via the Dyson equation 



G;^(k,u;)=c. + /x-7^^W,.(k)-S.( 



UJ 



(166) 



where 'H.f^j^p^ ^{k) is the LDA Hamiltonian in local orthogonal basis set depending on the 
Bloch vector k, and /i is the chemical potential. In order to avoid "double counting" , we can 
just subtract the static part of the self-energy, E'^'^ = TrSo-(O). It has been proven that this 
type of "metallic" double - counting is suitable for moderately correlated d electron systems 



( jLichtenstein et al. 



200 ih . 
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The standard QMC scheme for local Coulomb interactions takes int o account only 



density-density like inter actions, although the new continuous-time QMC (JRubtsov et al. 



2005 



Savkin et al. 



20051 ) can overcome this problem and include all the elements of inter- 



action vertex. We use the functional integral formalism and describe the discrete Hubbard- 
Stratonovich transformations for calculating the partition functions and corresponding 
Green's function. In this method the local Green's function is calculated for the imagi- 
nary time interval [0, (3\ with the r aesh r/ = /Ar, / = , ..., L — 1 (Ar = /5/L, (3 = 1/T) by 
using the path-integral formalism (JGeorges et al\ . \l99w . The multi-orbital DMFT problem 
with density-density interactions is described by the following effective impurity action 



S=- fdr f dT'Y.ct{r)g,,{r - r')c,{T') + ^ f dTY^n,{T)U.,n,{T) (167) 

h h ,^^ 2 Jo ,j 

where i = {m, a} labels both orbital and spin indices (we remind that we have no site 
indices since we are now solving the one-site effective impurity problem). Thus we truncate 
the original four-index rotationally invariant vertex and use only two-index approximation 
for it. This is a price we should pay for more exact way to solve the effective impurity 
problem. Without spin-orbital coupling we have Qij = Gm,mi^cTCT'- 

In the auxiliary fields Green-function QMC sch eme one can use the discrete Hubbard- 



Stratonovich transformation introduced by Hirsch (JHirschl . 
exp 



19831) 



'ArUij I niUj — 



-{ni + nj)] =7; Y. exp[AijSi 



[Ui 



^i)] 



(168) 



where Sij{T) are the auxiliary Ising fields for each pair of spins, orbitals and time slices with 
the strength: 



A 



u 



arccosh 



exp \ —Uij 



(169) 



Using Hirsch's transformation (I168p . (I169p we can transform the non- linear action to a 
normal Gaussian one (for a given configuration of the auxiliary Ising fields Sij) and integrate 
out exactly fermionic fields in the functional integr al ( I167p. As a res ult, the partition 
function and Green's function matrix have the form (jGeorges et ali . Il996l ) 

1 



Z 
G 



NfL 



2^/ 
1 1 



E det[G-\s.,)] 



Sij (t) 



Z2^f 



NfL 



J2 G{s,,)det[G-\s, 



(170) 



Sij (t) 
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where Nf is the number of Ising fields, L is the number of time shces, and G{sij) is the 
Green's function of non-interacting fermions for a given configuration of the external Ising 
fields: 



Gi'is) = g;,' + AA,s. 



A, 






1^ 






(171) 



Here we introduce the generalized Pauli matrix 



CTi 




(172) 



To calculate the Green's function G 



V 



can use the Dyson equation fiHirsch 



■s) fo r an arbitrary configuration of the Ising fields one 



19831): 



G' = [1 + (1-G)(e 



v'-v 



1)^0 



:i73) 



where V and G are the potential and Green's function before the Ising spin fiip, and V 
and G' after the fiip. The QMC important sampling scheme allows one to i ntegrate over 



the I sing fields with | det[G ^(5'jj)]| being a stochastic weight (JGeorges et al. 



1996 



Hirschl . 



19831 ). Using the output local Green's function from QMC and input bath Green's functions 
the new self-energy is obtained via Eg. (11601) . the self-consistent loop being closed through 
Eq. (ll62p . The main problem of the multi-band QMC formalism is the large number of the 
auxiliary fields S[^^,. For each time slice / it is equal to M(2M — 1) where M is the total 
number of the orbitals, which gives 45 Ising fields for d states case and 91 fields for / states. 
Analytical continuation of the QMC Green's functions from the im aginary time to the real 
energ y axis can be performed by the maximum entropy method (jjarrell and Gubernatid . 



19961 ). It is important to stress that for the diagonal Green's function Gij = GiSij the 



determinant ratio is always positive. This means that the sign problem, which is the main 
obstacle for the application of the QMC method to fermionic problems ( JTroyer and Wiesd . 
20051 ). does not arise in this case. Real computational experience shows that even for generic 
multiband case the sign problem for the effective impurity calculations is not serious. 

It is worthwhile to illustrate QMC scheme for correlation effects in the electronic struc- 
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tures of HMF by using a simple example. Since solving the full one-band Hubbard model 



T< = -J2 t^ji4aCja + c]^Cia) + UJ2 ^iT^ii 



(174) 



«J,o- 



is difficul t (see discussion in Section IIII.Ap . we treat the Dynamical Mean Field Theory 
(DMFT) (JGeorges et all Il996l ) which is formally exact in the limit of infinite dimensionality. 
Following this approach we will consider the Bethe lattice with coordination z ^ oo and 
nearest neighbor hoping tij = tj ^. In this case a semielliptic density of states is obtained 
as a function of the effective hopping t, A^(e) = (27rt^)~^V4t^ — e^- In order to stabilize the 
"toy" model in the HMF state, we add an external magnetic spin splitting term A which 
mimics the local Hund polarization. This HMF state corresponds to a mean-field (HF) 
solution denoted in Fig. [13] as a dashed line. 

The effective medium Green's function Q^j is connected wit h the local Green's fu nction 
Gf, on the Bethe lattice through the self-consistency condition (JGeorges et a/.l . Il996l ) 



Q-^ =%ujn^[i- t'G, - ll2a^ 



;i75) 



where Un = {'2n + l)7rT [n = 0, ±1,±2...) are the Matsubara frequencies. The 
Green's function corresponding to the DMFT effective action, Eq. (]167p . Gcr{T — t') = 
— {TrCa{T)cl{T'))s^ff, has bccu calculated us ing the Quant u m Mo nte Carlo scheme within 



1 r 



the so-called exact enumeration technique fiGeorges et aii. 
cretization parameter L = 25 (see for details ( 



Chioncel et al. 



996i). b v using the time dis- 



2003al )). We emphasize that 



due to the symmetry of the ferromagnetic state the local G^ and the effective medium 
Green's functions are diagonal in spin space, even in the presence of the interactions which 
enable the spin-flip scattering process. The magnon excitation can be studied through the 
two-particle correlation function 



xtAr-r') = {S^ir)S-{r')) = (T.c|(r)q(r)cl(r')ct(r'))5.,, 



(176) 



which is obtained by using the QMC procedure ( jJarreli Il992h . Being local, this function 
is insufficient to find the q-dependence of the magnon spectrum, but yields only a general 
shape of the magnon density of states. 

The DMFT results are presented in Fig. [131 In comparison with a simple Hartree-Fock 



solution one can see an additional well-pronounced feature appearing in the spin-down gap 
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region, just above the Fermi level: the non-quasiparticle states which are visible in both spin 
channels of DOS around 0.5 eV. In addition, a many-body satellite appears at 3.5 eV. 

The left inset of Fig. [121 represents the imaginary part of the local spin-flip susceptibility. 
One can see a well pronounced shoulder ( :^ 0.5 eV), which is related to a characteristic 



magnon excitation (jlrkhin and Katsnelson 



1985b 



199d . I1994J ) . There is a broad maximum 



at about 1 eV, which corresponds to the Stoner excitation energy. The right inset of Fig. 
T3l represents the imaginary part of self-energy. The spin-up channel can be described by a 
Fermi-liquid type behavior with a parabolic energy dependence — ImE^^ oc {E—EpY, whereas 
in the spin down channel the NQP shoulder at 0.5 eV is visible. Due to the relatively high 
temperature (T = 0.25 eV) in the QMC calculation the NQP tail extends below the Fermi 
level (remember that at zero temperature the tail should vanish exactly at the Fermi level, 
SectHlLBl). 



C. Spin-polarized T- matrix fluctuating exchange approximation 



Most of HMF materials are moderately or weakly correlated systems. This means that one 
can use some perturbative approaches which make computations much less expensive and 
allow one to work with the complete four-index Coulomb interaction matrix. A very efficient 
pertu rbative scheme has been proposed by Bickers and Scalapino (JBickers and Scalapind . 
19891 ) and was called fluctuating exchange approximation (FL EX). This was generalized to 
the multiband case and used in th e cont ext of the DMFT in (JKatsnelson and Lichtensteinl . 



1999 



Lichtenstein and Katsnelson 



19981 ). The latter step means that this approach is used 
not directly for the whole crystal, but for the effective impurity problem, so that the momen- 
tum dependence of the Green's functions is neglected. On the other hand, self-consistency 
of the DMFT procedure makes the description of the local effects in perturbative schemes 
more accurate. For example, in the case of one-band half-filled Hubbard model the second- 
order perturbation expression for the self-energy in the context of the DMFT provides cor- 
rect atomic limit and, actually, very accurate description of the metal-insulator transition 
( jKajueter and Kotliarl . Il996l ). A starting point in the FLEX approximation is the separa- 
tion of different interaction channels. The symmetrization of bare U matrix is done over 
particle-hole and particle-particle channels: 



u: 



mim, m2n^9 



21 p — IP 
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Tjm _ Tji 

Tjs _ _(Tn _|_ Tp \ 

^ miin'^m2m'2 r\\m-im'^m2m'2 ~^ ^ mim'-^m'2m2 1 

Tjt — _(Tp — TP ) n??"! 

^ m\va'-^m'2m2 r\\ni\m'^m'2'm2 miim'^ni2'm'2' \ ) 

The above expressions are the matrix elements of bare interaction which can be obtained 
with the help of the pair operators corresponding to different channels: 



particle-hole density 



particle-hole magnetic 



c?i2 = -7|(cuC2t + cl^C2i) (178) 



m° 



12 - -/|(CitC2T - Ci|C2j) 
^12 = C+TC2i (179) 



particle-particle singlet 



12 - Ci|C2t 



Sl2 = -y|(cuC2t-CitC2i) 



1 

72' 



'12 ~ /K^^li^'i'l '^1T'^2J,) (1"U) 



particle-particle triplet 



^?2 = -y|(cuC2t -CitC2i) 

^-12 — /^V^U''2T ''1T''2J.J 

^12 = Ci|,|C2J.,-f 

^12 = C]^| |C2|^^ (181) 

These operators describe the correlated movements of the electrons and holes below 
and above the Fermi level and therefore play an important role in defining the spin 
dependent effective potentials W^^^^^^^^. 

In the spin polarized T-niatrix fluctuating exchange approximation (SPTF) scheme 



( JKatsnelson and Lichtensteinl . |2002| ) the particle-particle interactions are described in the 
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T-inatrix approach (JGalitskil . Il958l : iKanaraoril . Il963l ) where for the effective interaction the 
sum over the ladder graphs are carried out with the aid of the so called T— matrix which 
obeys the Bethe-Salpether-like integral equation: 

(13|T'"'^'(ifi)|24) = (13|u|24) 

- ^Y.T.Gl%{iuj)G\''s{iuj)G{in-iuj){68\T'"'\in)\24) (182) 

^ oJ 5678 

The corresponding contributions to the self-energy is described by the Hartree and Fock 
diagrams with the formal replacement of the bare interaction by the T— matrix: 



sSi;^(-) 



^'Z\^^) 



^EE(i3|r'^'^'(^fi)|24)Gi;(.fi 



p 



lUJ] 



n 34(t' 



P n SAcr' 



(183) 
(184) 



The "Hartree" contribution dominates for small concentration of electrons and holes; these 
two contributions together contain all second-order in V terms in the self-energy. 

Combining the density and magnetic parts of the particle-h ole channel we can write down 
the expression for the interaction part of the Hamiltonian (JKatsnelson and Lichtensteinl . 



1999 



Lichtenstein and Katsnelson 



19981 ). 

Hu = D*H^ * D- 



1 



rtu = -Tt{D+ *V^UD + m+ *V^*m~ +m- *V^* m+) (185) 

where D is a row vector with elements {d,m^), and D~^ is a column vector with elements 



((i+TTi^), * stands for the matrix multiplication with respect to the pairs of orbital indices. It 
follows from the model consideration presented above that for proper describing the effects 
of electron-magnon interactions it is important to replace the bare spin-flip potential by a 
static limit of the T-matrix (see Eqs.(l6]) and (l35l)). With this replacement, the expression 
of the effective potential is: 

(V;;[)i234 = (13|TTi|42) (186) 



The matrix elements of the effective interaction for z (longitudinal) spin fluctuations are: 

Vdd = iE(E(13|T'^'^'|42)-(13|r-'-'|42)) 
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Vdm = Vmd = iE^((13|T'^l42)-(13|T-'^|24) + (13|T-'-|42)) 

Vmrn = ^ J^lJ] aa'(13|T'^'^' |42) - (13|T'^''^'|42)) (187) 

Further we introduce the expressions for the generahzed longitudinal (x") and transverse 
(X"*") susceptibilities 

X^ii^) = [l + V^T^^{iuj)]-^*T^^{iuj) (188) 

XH^^) = [1 + V^" * xl{^io)]-' * x!(^cu) (189) 

where r[iuj) is the Fourier transform of the empty loop: 

mim2m3m4,\ ) m2»Tt3V /m4miV / ^^lyuj 

and the bare matrix longitudinal susceptibility is: 

1 / r^^ + r^^ r^^ - r^^ \ 

xU^uJ) = -\ ^^ ,, ^^ ,, (191) 

The four matrix elements correspond to the density- density (dd), density-magnetic {dm^), 
magnetic-density {mPd) and magnetic-magnetic channels (m°m°) and couple longitudinal 
magnetic fluctuation with density magnetic fluctuation. In this case the particle hole con- 
tribution to the self-energy is: 



n2a ' 



T 



E W,%Ut)G-,: (192) 



34o-' 



with the particle-hole fluctuation potential matrix 



W^'iiu) =1 " 'M (193) 



and the spin-dependent effective potentials are deflned by 



^v^" * (x" - x!) * v^" 



Wu = ^v^'i * (x" - x!) * v^ii 
Wn = \v^ * ix"- - xt) * Ka 

Wn = \v^ * ix-^ - Xo +) * V.,^ (194) 
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where x"? Xo differ from the values of xK Xo by the replacement r^^^^ ^^ T^^ in equation (11911) . 
The final expression for the SPTF self-energy is given by: 



s = s(^^) + s(^^) + s(p^) 



(195) 



Due to off-diagonal spin structure of the self-energy T,^p^^ this method can be used to 
consider the non-quasiparticle stat es in HMF. A more detailed justi fication and descrip- 
tion of the method is presented in (JKatsnelson and Lichtensteinl . |2002| ). This approach was 
also generalized to the case of strong spin-orbit coupling and used for actinide magnets 



( jPourovskii et al. 



2005 



20061 ). In that case separation into density and magnetic channels 



is not possible and one should work with the 4x4 supermatrices for the effective interaction. 
Recently several LDA+DMFT calculations of different HMF materials have been car- 
ried out. Now we review briefly the results of these calculations focusing mainly on non- 
quasiparticle states resulting from the electron-magnon interactions. LDA+DMFT is the 
only contemporary practical way to consider them in the electronic structure calculations. 

V. ELECTRONIC STRUCTURE OF SPECIFIC HALF-METALLIC COMPOUNDS 



A. Heusler alloys. 

1. NiMnSb: electronic structure and correlations 

The intermetallic compound NiMnSb crystallizes in the cubic structure of MgAgAs type 
(Clb) with the fee Bravais lattice (space group FA3m = Tj). The crystal structure is shown 
in Fig. [TH This structure can be described as three interpenetrating fee lattices of Ni, Mn 
and Sb. A detailed description of the band str ucture of semi-Heusler a ll oys was given by 



using electronic structure calculations analvsis (|de Groot et al 



2002a: 



2003 



Galanakis and Mavropoulos 



Ogiit and Rabe 



2007 



Kulatov and Mazin 



1983b ). f 



2003 



Galanakis et al. 



Nanda and Dasguptal . 



19951 ): we briefly summarize the results. 
To obtain the minority spin gap, not only the Mn— d-Sb— p interactions, but also Mn— d- 
Ni— (i interactions are to be taken into account. Moreover, the loss of inversion symmetry 
produced by CI;, structure (the symmetry lowering from Oh in the L2i structure to Td in the 
Clfe structure at Mn site) is an essential additional ingredient. All the above interactions 
combined with the T^ symmetry lead to a nonzero anticrossing of bands and to gap opening. 



87 



The large exchange sphtting of Mn atom (producing the local Mn magnetic moment of 
about 3.7 fis) is crucial to induce a half-metallic structure. In the spin-polarized calculation 
the position of t2g and Cg Ni states is changed slightly, so that the exchange splitting on Ni 
is not large. The local Ni magnetic moment calculation gives a value around O.SfiB- 

The non-spin-polarized result has a striking resemblance to the majority spin-polarized 
calculations, presented in Fig. [161 Kulatov et al. explained half-met allic ity of NiMnSb by the 
exten ded Stoner factor calculations in the rigid-band approximation (JKulatov and Mazinl . 
20031 ): the minority spin-band gap opens due to the exchange splitting which shifts minority 
bands, so that they become empty. 



There is an excel 



( Ide Groot et ali Il983b[) and recent ^ 



ent agreement between the first band-structure calculation 



gations (JAndersen and Saha-Dasgupta 



-th o rder muffin-tin orbita l s (NMTO) investi- 



2000 



Yamasaki et al. . 



2006 



Zurek et al. 



20051) 



NMTO Wannier Mn d orbitals are shown in Fig. [151 The triple-degenerate manganese 
t2g orbitals are very complicated due to the hybridization with Ni d and Sb p states. The 
dxy orbital at Mn site is deformed by antibonding with Ni d state directed tetrahedrally to 
[111], [ill], [111], and [111]. The same Ni d orbitals couple with Sb p states. The direct 
yhi-dxy-'^h-p vr coupling is not seen since the distance is (i(Mn-Sb):(i(Ni-Sb) = 1 : v^/2. 
Therefore the Ni-rf-Sb-p interactions are more favorable. The dispersion of the Mn t2g bands 
is mainly due to hopping via the tails of Sb p and Ni d orbitals. On the other hand, the 
next-nearest neighbor (n.n.n.) d — d hopping of t2g orbital is small. The Cg orbitals at Mn 
site are much easier to understand: they point towards Sb atoms, and a strong pda coupling 
between Sb p and Mn e^ states is seen. This i nduces large n.n . n. d-d hoppings. 

The Wannier or bital analysis for NiMnSb (j Yamasaki et alV |2006| ) confirms the previous 



conclusions of Ref. dde Groot et al 



the role of d-d hybridization (lGalanakis et al. 



1983br) about t he role of p-d hybridization as well as 
a|). The LDA partial density of states for 



2002£ 



half-metallic NiMnSb and the fat-band structure which marks the main orbital character of 
valence states are presented in the Fig. [T7I1T81 [T9l One can see that the spin-up Mn d and 
Ni d states are located at the same energy region, while spin-down bands are separate due 
to the significant Mn exchange splitting. The top of the valence spin-down bands form by 
the Mn t2g, Sb p and Ni t2g orbitals, while bottom of the conduction bands are due to the 
Mn Cg and t2g states. This means that large Mn spin-splitting and Sb mediated indirect 
Mn-Ni interactions are responsible for formation of the half-metallic gap. One can see this 
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complicated Mn t2g valence orbitals with large contributions of Sb p and Ni t2g states in the 
Fig. flTSj) . The physical picture of the downfolding analysis should not change very much if 
one explicitly includes Sb p and Ni d orbitals. 

We discuss now in detail the prototype half-metallic ferromagnet NiMnSb where the 
gap is situated in the spin-down (minority) channel. The temperature dependence of the 
HMF electronic structure and stability of half-metallicity against different spin-excitations 
are crucial for practi cal applications in spintronic s . A simple attempt to incorp orate finite- 



temperature effects (JDowben and Skomski 



2003 



Skomski and Dowben 



2OO2I), leading to 



static non-collinear spin-configurations, shows a mixture of s pin up and spin down density 
of states that destroy the half-metallic behavior. In Ref. (jChioncel et all [2003aj) a dif- 
ferent more natural approach was used to investigate the proper effect of dynamical spin 
fluctuations on the electronic structure at T < Tc, within the half-metallic ferromagnetic 
state. 



The LDA-FDMFT calculation for NiMnSb (JChioncel et all l2008al ) was the first appli- 
cation of the combined electronic structure and maiiy-body technique to HMF. Aryaseti- 
awan et al pointed out recently ( jAryasetiawan et all 1200J) that a rigorous way to define 
the screened frequency- dependent on-site Coulomb interactions matrix elements for corre- 
lated states related with generalized GW scheme where d — d screening is suppressed to 
preserve the double counting of in the model approach. However, in practice for realistic 
materials, the elimination of degrees of freedom is a very difficult procedure. To find the 
average Coulomb interaction on the d atoms U and corre sponding exchange interactions J , 



a more simple approa c h, the constrained LP A scheme flAnisimoy 



Dederichs et al. 



1988 



1984: 



Gunnarsson et al. 



Norman and Freeman 



1989 



Hybertsen et al. 



and Gunnarssod. 



198 



I 



1991 



McMahan et al. 



19861 ). can be used. In this approach the Hubbard U is cal- 
culated from the variation of the total energy with respect to the occupation number of 
the localized orbitals. In such a scheme, the metallic sc reening is rather inefficient for "i d 



transition metals, and effective U is of the order of 6 eV (lAnisimov and Gunnarssonl . Il99ll ) 



The perfect metallic screening will lead to a smaller value of U . Un 



no re liable scheme to calcul ate U within constrained LDA for metals (jSolovyev and Imadal . 



2OO5I ). and in the works (JChioncel et al. 



2006a 



2003a 



2005 



brtunately, there are 



2006bl ) some intermediate 



values were chosen, f/ = 2 -^ 4 eV and J = 0.9 eV. The recent analysis of angle- resolved 
photoemission and the LDA theory indicates a shift of spectral function of the order of 
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0.5 -r- 1 eV, which can be attributed to corre^ 



ation effects beyon d LDA ( ICorrea et al\ . \200w . 



First angle-resolved photoemission results (JKisker et all 119871 ) in generally agree with the 



LDA band structure, but demonstrate the same shift of quasiparticle dispersion of the or- 
der of 0.5 eV below the Fermi level for spin-down Mn t2g bands. This can easily lead to 
the effective Hubbard interactions in the static mean-field approximation of the order of 
U* = 1 eV, although one should carefully investigate the effects of spin-orbital splitting. 
Since the spherically-average effective Hubbard interactions is equal to U* = U — J and 
the value of intra-atomic exchange interactions is not screened much in solid s being of the 
order of J = 1 eV, we can conclude from the photoemission experiments (JKisker et al . 



19871 ) and resonant x-ray scattering (JYablonskikh et all 1200 ll ) that the Hubbard interac 



tions U = 2 eV is quite reasonable. An account of the correlation effects in the framework 
of the L DA-I-DMFT method also improves the description of magnetooptical properties of 



NiMnSb fIChadov et al. 



20061 ). 



The typical results for density of states using LDA and LDA-I-DMFT is presented in Fig. 
I20l The LDA-I-DMFT density of states shows the existence of non-quasiparticle states in 
the LDA gap of the spin down channel just above the Fermi level. It is important to 
mention that the magnetic moment per formula unit is not sensitive to the U values. For 
a temperature T = 300 K the calculated magnetic moment, /x = 3.96/iB, is close to the 
zero-temperature LDA- value which is integer, /i = 4 yU^. This means that the half- metallic 
state is stable with respect to switching on the correlation effects. The DMFT gap in the 
spin down channel, defined as the distance between the occupied part and starting point of 
NQP tail, is also not very sensitive to U. The total DOS is also weakly [/-dependent due to 
the T-matrix renormalization effects. 

In comparison with the LDA result, a strong spectral weight transfer is present for the 
unoccupied part of the band structure due to appearance of the non-quasiparticle states in 
the energy gap above the Fermi energy (Sect JIH^Bj) . Their spectral weight is not too small 
and has a relatively weak depend ence on the U value fFig. [2T|) , whic h is also a consequence 



2OO2I). 



of the T-matrix renormalization (iKatsnelson and Lichtensteinl . 

For spin-up states we have a normal Fermi-liquid behavior — ImS^(£') (x {E — EpY with 
a typical energy scale of the order of several eV. The spin-down self energy behaves in a 
similar way below the Fermi energy, with a bit smaller energy scale. At the same time, a 
significant increase in 1\ivTj\i{E) with a much smaller energy scale (tenths of eV) is observed 
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right above the Fermi level, which is more pronounced for t2g states (Fig. [201) . The NQP 
states are visible in the spin-down DOS, Fig. [20], at the same energy scale as the imaginary 
part of S^^. A similar behavior is evidenced in the model calculation of Fig. [121 The NQP 
spectral weight in the density of states (Fig. [2U1) is proportional to the imaginary part of 
the self-energy. 

From the general many-body theory, the DMFT approach is neglecting the momentum 
dependence of the electron self-energy. In many cases, such as the Kondo effect and the 
Mott metal-insulator transition the energy dependence of the self-energy is obviously more 
important than the moment um dependence and, therefore, the DMFT scheme is adequate 
to consider these problems (JGeorges et all Il996l ). In the case of itinerant electron ferro- 
magnetism, the situation is much less clear. However, the LDA+DMFT treatment of finite 
temperature magnetism a nd electronic structure in Fe and Ni appeared to be quite successful 



( jLichtenstein et al. 



20011 ). Experimentally, even in itinerant electron paramagnets which are 



close to ferromagnetic instabil ity, such as Pd, t he momentum dependence of the self-energy 
does not seem to be essential (jJoss et all Il984l ). One can expect that in magnets with well 
defined local moments such as half-metallic ferromagnets the local DMFT approximation 
for the self-energy should be even more accurate. In particular, as discussed above, it can 
be used for the calculations of spin-polaronic or non-quasiparticle effects in these materials. 



2. Impurities in HMF: lanthanides in NiMnSb 

Here we discuss the introduction of non-intrinsic deffects in NiMnSb, which preserves 
the half-me tallic properties on one hand while optimizing the magnetic disorder on the 



other hand (lAttema et al 



( lAttema et al. 



2004; 



Chioncel 



2004). Go od candidates are the rare-earth impurities in NiMnSb 



2004J ). The motivation of this choice is connected with the 
substantial spin-orbit interaction in the rare-earth localized / shell. Besides that, the origin 
of the band gap in NiMnSb is closely related to the band gap in III-V semiconductors: it is 
expected that a substitution of some of the tetravalent element in NiMnSb by a lanthanide 
preserves the essential feature of the half-metal: the band gap for one spin direction. 

These systems can be really synthesized, compounds RNiSb (for heavy rare earths R) 
existing with exactly the same crystal structure as NiMnSb. The rare-earth atoms which 
show both a large spin and orbital moment can be expected to introduce a large spin-orbit 
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interaction, in other words around 1/4 (Nd) and 3/4 (Er) in the lanthanides series. Total 
energy calculations allow one to evaluate the coupling between the rare-earth (4/) impurity 
spins and the manganese {3d) conduction electron spins. For a substantial coupling the 
fluctuations of the Mn {3d) conduction electron spins, i.e. the spin wave, might be blocked, 
thus the magnon branch is qualitatively changed. In contrast with the clean limit (pure 
NiMnSb), the magnon spectra of NiMni_a;R^Sb present a fragmented structure with several 
gaps in the Brillouin zone. This fragmentation implies that the finite temperature effects 
are diminished for a suitably chosen rare-earth. 

The ab-initio electronic structure calculations were carried out using the scalar relativistic 



linear muffin-tin orbital (L MTO) met 



flavors: LDA and LDA+f/ (jAndersen 



l od within the atomic-sphere 



1975 



Andersen and Jepsen 



appr o ximation in two 



1984 



Anisimov et al. 



1997al ). To evaluate the coupling between the rare-earth and the Mn sublattices, ferro- 
and ferrimagnetic structures were taken as the initial state of the calculation and they were 
preserved after the self-consistent calculation. 

A simplifled model was used that captures the complex interplay of the Mn {3d) itiner- 
ant conduction band electrons and the localized 4/ electrons, the latter carrying a strong 
magnetic moment. This model deals with a mean-fleld decoupling, in which the Mn 3d and 
the R / states are described by the LDA+?7 whereas the 3d — 4/ interaction is treated as 
perturbation. The corresponding mean-fleld Hamiltonian can be written in the form: 



H-T^LDA+U-^E^f'^i 



+5 



(196) 



where the spin of the conduction electron at site Ri is denoted by af"' and Sfj^^ represents the 
spin of the 4f shell at the Ri+s site. The Mn d local moment fluctuations could be quenched 
by a strong f — d coupling, which affects the magnon excitations. The strength of such a cou- 
pling was evaluated by calculating in an ab-initio fash ion the total energy of NiMui ^ ^R^S b 



compounds for a parallel/antiparallel f — d coupling (jAttema et al. 



2004 



Chioncel 



20041 ). 



Given the geometry of the cell the lanthanide substitution is realized in the /cc-Mn sub- 
lattice, so 12 pairs of R(4/)-Mn(3(i) are formed. The f — d coupling is calculated as the 
i?ll — Ei^ energy corresponding to a pair. Adopting a two sublattices model described by the 
Hamiltonian (I196p . i.e. the sublattice of lathanides R(4/) being antiferro/ferromagnetically 
oriented with respect to the Mn(3ci) sublattice, the J values correspond to the intersublattice 
couplings. In pure NiMnSb the ferromagnetic Curie temperature Tq = 740K is determined 
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by the strength of the Mn-Mn {3d — 3d) sublattice interaction. For a small lanthanide 
content one can expect that this interaction in NiMni_a;Ra;Sb compounds is on the same 
scale as in pure NiMnSb. Therefore the substitution introduces the competition between the 
intra- and intersublatice interactions which are crucial parameters for any practical applica- 
tions. In the case of Nd there is a large 3d — 4/ coupling which should dominate over the 
3d — 3d coupling. For temperatures lower than the Curie temperature corresponding to the 
NigMnyNdSbs compound, T^'^, the 3d — Af coupling interaction will lock the Mn {3d) mag- 
netic moment fluctuations decreasing the available number of magnon states. Above T^^ 
the thermal fluctuations already took away the long range order of the Mn {3d) sublattice, 
and there are no available magnon states at all. As a consequence, the Nd substitution can 
be attractive for the high-temperature applications where the 3c? — 4/ coupling might play 
an important role. In the cubic structure of NiMnSb the uniaxial anisotropy is completely 
missing. Nevertheless, the 4/ impurity spin-orbit coupling lifts the degeneracy in the F 
point. The lowest 3d — 4/ coupling is realized in the case of Ho. 

The LDA+U density of states for HoNi8Mn7Sb8 for 3d — 4/ antiparallel and par- 
allel couplings and f// = 9 eV value of the on site Coulomb interaction is presented 
in Fig. fl22|) . On one hand for the antiparallel coupling the magnetic moments are 
/i^f = — 4.09/iB,/iMni = 3.73/iB,/U^^2 = 3.80/iB and the gap of 0.55 eV is situated in 
the majority spin channel. On the other hand for parallel coupling the magnetic moments 
have almost the same magnitude /^ho = 3.96/xb,//^„]^ = 3.72fiB, fJ'Mn2 = 3.80/iB with a 
similar gap situated in the minority spin channel. It is important to mention that in 
the case of NigMnyHoSbg the Ho(4/) orbitals do not hybridize with the Mn{3d) orbitals 
near the Fermi level. As one can see in Fig. (!22l) . the behavior of DOS near the Fermi 
level is very similar, so that the nature of carriers around Ep is not changed. The elec- 
tronic structure calculation for the a ntiferromagnetic HoNiSb compound was performed in 



(JAttema et al. 



2004 



Chioncel 



20041 ). This is known to be a se miconducting material wi th 



interesting transport pr operties (giant magnetoresistance effect ( JKarla et al. 



Pierre and Karlal. 



down to lOK (JKarla et al 



1999 



1998^ 



2000D). 



Jlbl; 



he i n verse s us ceptibility curves show a Curie- Weiss behavior 



1999 



1998a 



a 



Pierre and Karla 



2000l ). The onset of the an- 



tiferromagnetic ordering was estimated from susceptibility measurements between 1.5 and 
2.5 K, and the neutron diffraction data in dicated an antiferromag netic propagation vector 



;i/2, 1/2, 1/2) flKarla et al. 



1999 



1998al Jbl: iPierre and Karlal . l2000l ) . In order to describe the 
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antiferromagnetic ground state, a rhoinbohedral description of the Ho2Ni2Sb2 unit cell was 
used. The atoms Hoi (0,0,0) and Ho2(l, 1, 1) acquire magnetic moments of fiuoi = 4.0//s 
and /i§02 = — 4.0/iB, respectively. The antiferromagnetic insulating state shows a gap of 
0.29 eV in agreement with the experimental results. There are some qualitative features 
which are similar for the antiferromagnetic Ho2Ni2Sb2 and the HoNi8Mn7Sb8 compounds: 
the almost identical magnetic moment /xho and positions of occupied and unoccupied Ho(4/) 
peaks in DOS. Perhaps the most important observation is that the spin-down channel in 
Ho2Ni2Sb2 is isoelectronic to that in NiMnSb, so that the Ho substitution would preserve 
the half-met allicity of HoNi8Mn7Sb8 in the minority spin channel. 

3. FeMnSb: a ferrimagnetic half-metal 



Early theoretical studies demonstrated that the gap in the minority spin channel is sta- 
ble with respect to change of t he 3d atom X = Fe, Co, Ni in the XMnSb compounds 



(Ide Groot et al. 



1986 



Kubler 



19841 ). Noticeable difference between Ni- and Fe-based 



Heusler alloys is that NiMnSb is a half-metallic ferromagnet with a very small value of 
Ni magnetic moment (0.2/2^), whereas in FeMnSb the antiferromagnetic coupling between 
Fe (— 1/is) and Mn (S/xr) moments stabi lizes the gap and the half-metallic ferrimagnetic 
electronic structure (jde Groot et a/.l . 119861 ). Unfortunately, the ternary compound FeMnSb 
does not exist, but indications of its magnetic and crystallographic p roperties were obtained 
by extrapolating the series of Nii_2;Fea;MnSb (Ide Groot et all Il986l ) to high Fe concentra- 
tion. 

In the non-relativistic approximation there are two essentially different sources for states 
in the gap at finite temperatures. First, there is the simple effect of gap filling due to 
disorder, i.e. due to scattering on static (classical) spin fluctuation or thermal magnons; 
this is symmetric with respect to the Fermi energy. On the contrary, the correlation effects 
result in an asymmetry in the gap filling, spin-down non-quasiparticle states appearing above 
the Fermi level. One has to take into account also spin-orbit coupling effects mixing the 
spin-up and spin-down states due to non-zero elements of spin-orbital interactions V^q (see 
the discussion in Sect JIII.Bl) . To illustrate the differences between the static and dynamic 



effects we plot the DOS of the LDA- 



recent results including SO coupling (Ide Groot et al. 



DMFT calculations ^ yhich should be compared with 



1986 : 



Mavropoulos et a/.l . l2004j ). 
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A relatively weak dependence of the NQP spectral weight on U (Fig. [2^ is evidenced for 
both NiMnSb and FeMnSb compounds. A "saturation" of the spectral weight for FeMnSb 
takes place for almost the same value, f/* ~ 1 eV, as in the case of NiMnSb, which is in 
agreement with exp erimental observati on of relatively weak correlation effects in Heusler 



HMF like PtMnSb flKisker et al. 



19871 ). The small value of effective Hubbard parameter 



can be understood in terms of the 



tions (JKatsnelson and Lichtenstein 



arge T-mat rix renormalization of the Coulomb interac- 



1999 



2OO2I ). The spectral weight values 



br FeMnSb are 



larger in comparison with those obtained for NiMnSb (jChioncel et a/.l . l2003al ). which can be 
attributed to a larger majority spin DOS at the Fermi level. 

Th e spin-orbit coupling prod uces a peak in the minority-spin channel close to the Fermi 



level ( JMavropoulos et al. 



20041 ). which is an order of magnitude smaller than the spectral 



weight of the NQP states. According to the SO results, the polarization at the Fermi level 
for Ni MnSb and FeMnSb are almost the same. In contrast, the calculation ( Chioncel et al\ , 
2006al ) shows that the spectral weight of the NQP states in FeMnSb is almost twice as large 
as the value for NiMnSb. 

To discuss the influence of temperature and local Coulomb interactions on the polarization 
in FeMnSb compound, we present results of LDA+DMFT calculations for T < 400 K and 
different f/'s. Figj25] presents the contour plot of spin polarization P(-E', T) as a function of 
energy E and temperature T for U = 2 and 4 eV. The LDA value, plotted for convenience 
as the T = K result, shows a gap with magnitude 0.8 eV, in agreement with previous 
calculations ( JMavropoulos et a/.l . I2004J ) . 

One can see a peculiar temperature dependence of the spin polarization. The NQP fea- 
tures appear for E — Ep > and are visible in Fig. [2l] for U = 2 eV and T = 300K. 
Increasing the value of U from 2 eV to 4 eV, the NQP contribution in depolarization be- 
comes more significant. When the tail of the NQP states crosses the Fermi level, a drastic 
depolarization at Fermi level takes place, for f/ = 4 eV the NQP contribution being pinned 
to the Fermi level. 

One can see a clear distinction between the finite-temperature behavior of the polarization 
and magnetization, shown in Fig. [26] for different values of U. For U = 4eV, already at 
lOOK there is a strong depolarization about 25%. On contrary, it is interesting to note that 
the reduced magnetization M(T)/M(0) decreases slowly in the temperature range shown in 
Fig. [251 This reduction is a consequence of the finite temperature excitations, i.e. spin- 
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flip processes, affecting both spin channels. In the minority spin channel, NQP states are 
formed, and in the majority channel a spectral weight redistribution around the Fermi level 
(FiglMD contributes to the depolarization. The corresponding depolarization increases with 
the strength of correlations. The density of NQP states displays a rather strong temperature 
dependence (Sect JIII^Uj) and results in an asymmetry that is visible in Figs. [M] and [23 

NQP states dominate in the depolarization of this class of Heusler compounds, while spin- 
orbit contributions are much smaller. In addition, many-body effects are more pronounced 
in FeMnSb than in NiMnSb. This is connected with the larger DOS in the majority spin 
channel in the former material. Therefore, doping of NiMnSb by Fe could be an interesting 
issue to investigate the interplay between alloying and many-body effects. The LDA+DMFT 
calculations for NiMnSb supercell containing 25 % Fe impurities show a half-metallic char- 
acter at the LDA level, with the same strong correlation-induced depolarization effects as 
in pure FeMnSb. Therefore many-body effects for this material are of primary importance 
even in the presence of disorder. Correlation effects on surfaces of half-metals were discussed 
recently and it wa s shown that these states can be probed both directly and via their effect 



on surface states filrkhin and Katsnelsonl . l2006l ). 



4- Co2MnSi: a full- Heusler ferromagnet 

The origin of minority band g a p in full Heuslers was discussed by Galanakis et al 



(JGalanakis and Mavropoulos 



2007 



Galanakis et al. 



20061 ). Basing on the analysis of the 



band structure calculations, it was shown that the 3d orbitals of Co atoms from the two 
different sublattices, Co^(0,0,0) and Co^(l/2, 1/2, 1/2), couple and form bonding hybrids 
Co^{t2g/eg) - Co^{t2g/eg). In other words, the t2g/eg orbitals of one of the Co atom can cou- 
ple only with the t2g/eg orbitals of other Co atom. Further on, the Co-Co hybrid bonding 
orbitals hybridize with Mn{d)-t2g,eg manifold, while the Co-Co hybrid antibonding orbitals 
remain uncoupled owing to their symmetry. The Co- Co hybrid antibonding t2g is situated 
below the Fermi energy Ep and the Co-Co hybrid antibonding Cg is unoccupied and lies just 
above the Fermi level. Thus, due to the missing Mn{d)-t2g,eg and Co-Co hybrid antibonding 
hybridization, the Fermi energy is situated within the minority gap formed by the triply 
degenerate Co-Co antibonding t2g and the double degenerate Co-Co antibonding Cg. 

Fig. [28] shows the results of DOS calculations using the LDA and LDA+DMFT schemes. 
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The inset presents the spin polarization P{Ef) = {N^{Ef) - Ni{Ef))/{N^{Ef) + Ni{Ef)). 
One can see that in the minority spin channel asymmetric NQP states are formed, while in 
the majority a spectral weight red istribution takes place , which contributes to the depolar- 



ization. Contrary to the FeMnSb (iChioncel et all l2006al ). where the density of NQP states 
shows a rather strong temperature dependence, in the full-Heusler Co2MnSi the temperature 
depen dence is not so significant, similar to the result obtained for NiMnSb (IChioncel et al . 
2003af ). 



B. Half-metallic materials with zinc-blende structure 

One of the strongest motivations to investigate magnetic semiconductors and half- metallic 
ferromagnets is the possibility to design and produce novel stable structures on semiconduct- 
ing substrates with new interesting properties. From this point of view first-principle studies 
are an excellent starting point to predict new systems having the desired properties. Using 



full-potential density-functional method, all the 3d transition metal p nictides and cha 



genides with wurtzite structure were investigated systematically in Ref . (IXie et all l2003al ) in 



co- 



order to find half-metallic ferromagnets. These can be fabricated as thin films with thickness 
large for real spintronic applications. Nine of the wurtzite phases, MnSb, CrAs, CrSb, VAs, 
VSb, CrSe, CrTe, VSe, and VTe, were found to be robust half-metallic ferromagnets with 
large half-metallic gaps (0.2-^1 eV). Being compatible with the Ill-V and II- VI semicon- 
ductors, these half-metallic ferromagnetic phases, when realized experimentally, would be 
useful in spintronic and other applications. At the same time, zinc-blende (ZB) phases of 
MnAs, CrAs and CrSb have been fabricated successfully in the form of nanodot s, ul t rathin 



films and multilayers, respectively. A study within density-functional theory (JL 



LhJ, 



20031) 



predicted for the ZB CrSb half-metallic ferromagnetism with a magnetic moment of 3.0 fiB 
per formula. 



1. CrAs: tunable spin transport 



Recently Akinaga et al. ( lAkinaga et all |2000| ) found the possibility to fabricate ZB- 
type CrAs half-metallic ferromagnetic material. Experimental data confirmed that this 
material is ferromagnetic with the magnetic moment of SfiB, in agreement with theoretical 
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predictions (jAkinaga et a/I 120001 ) . According to this calculation, this half-metallic material 
has a gap of about 1.8 eV in the minority spin channel, which has attracted much attention 
to this potential candidate for spintronic applications, keeping in mind also its high Curie 
temperature Tc about 400K. Recent experiments on CrAs epilayers grown on GaAs(OOl) 
evidenced an orthorhombic structure, different from the ZB one, so the structure is rather 



sensitive to the preparation process (JEtgens et all |200J). However, it is highly desirable 
to explore the possibility of existence of half-metallic ferromagnetism in materials which 
are compatible with practically important III-V and II-IV semiconductors. For this purpose 



efforts have been made to inv estigate metastable ZB structures, such as CrAs (JAkinaga et al. 



2000 



Mizuguchi et al. 



2OO2I). 



It is interesting to explore the mechanisms of half-metallic ferromagnetism at fini te tem- 



perat ure from a realistic electronic structure point of view. Theoretical studies (IShirai 



2OO3I ) of the 3d transition metal monoarsenides have shown that the ferromagnetic phase 
of ZB structure CrAs compound should be more stable than the antiferromagnetic one. 
The many-body effects are very sensiti v e to structural properties of the artificially pro- 
duced CrAs compound (jChioncel et al\ . 120051 ). Similar electronic structure calculations 
concerning the stabilit y of the h alf-met allic ferromagnetic state in the ZB structure 



have been carried out (IXie et al. 



2003bl ). The LDA+DMFT calculations were carried 



out for three lattice constants: the GaAs (5.65A), InAs (6. 06 A ) and t 



value (oeq = 5.8A) obtained by density functional calculations fIShirail. 



le "e quilibrium" 



responding LDA computational results agree with previous ones (jAkinaga et all 



Galanakis and Mavropoulos 



2003 



2003h. The cor 



2000; 



Shirail . I2OO3I ). In order to evaluate the average Coulomb 



interaction on the Cr atoms and the correspond ing exchange intera ctions, the cons t raine d 
LDA method (jAnisimov and Gunnarssonl . Il99ll ) was used in Ref (JChioncel et all |2005| ). 
which yielded U = 6.5 eV and J = 0.9 eV. It is important to note that the values of the 
average Coulomb interaction parameter sl ightly decrease going from the GaAs {U = 6.6 eV) 



to InAs {U = 6.25 eV) lattice constants, (jChioncel et all l2005l ) which is in agreement with 



a naive picture of a less effective screening due to increasing of the distances between the 
atoms. 

The typical insulat i ng screening used in the constraint calculation 



(JAnisimov and Gunnarssonl . Il99ll ) should be replaced by a metallic kind of screening. 



The metallic screening will lead to a smaller value of U. Since there are no reliable schemes 
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to calculate U in metals, some intermediate value were chosen, U = 2 eV and J = 0.9 eV. 
It is important to realize that there are no significant changes in the values of the average 
Coulomb interaction for the lattice structures studied, the exchange interaction being 
practically constant. Note that the p hysical results are not very sensitive to the value of U, 



as it was demonstrated for NiMnSb (IChioncel et al. 



2003ah . 



The LSDA and LSDA+DMFT calculation for the density of states, DOS, is presented 
in Fig. [30l Depending on the lattice constant the Cr and As atoms loose electrons and 
this charge is gained by the vacant sites. As a result, the Fermi level is moved from the 
right edge of the gap (for the GaAs substrate) towards the middle of the gap (for the InAs 
substrate). The Cr moments are well localized due t o a mechanism similar to localization 
of the magnetic moment on the Mn atom in NiMnSb (jde Groot et a/.l . ll983bl ). The local Cr 
spin moment is more than S/i^, and the As atom possesses a small induced magnetic moment 
(of order of —0.3/1^) coupled antiparallel to the Cr one. The results are presented in table 
[H Cr magnetic moments calculated in DMFT increase in comparison with the LDA results 
due to the localization tendency of the Cr 3(i states, which is a consequence of correlation 
effects. 



According to the calculations of Ref. (IChioncel et al. 



20051 ). the system remains half- 



metallic with a rather large band gap (about 1.8 eV) for a ll the lattice constants compared 



with the band gap of the NiMnSb which is only 0.75 eV (jde Groot et al 



1983bh . In Fig. 



the non-quasiparticle states are visible for lattice parameter higher than the equilibrium 
one, with a considerable spectral weight in the case of the InAs substrate. This situation is 
very favorable for the experimental investigation of the NQP states. 

Comparing the electronic structure of CrAs growing on InAs or GaAs substrates we 
can conclude that the most significant change in the electronic structure is related with 
the As p states. Having a larger lattice constant in the case of InAs substrates, the Cr 
atom acquires a slightly larger magnetic moment. Nevertheless, in the LDA calculations 
the magnetic moment per unit cell is integer, 3/iB- Expanding the lattice constant from 
GaAs to InAs lattice the Cr states become more "atomic" , and therefore the spin magnetic 
moment increases. This is refiected equally in the charge transfer which is smaller for the 
InAs lattice parameters. A larger Cr moment induces a large spin polarization of the As p 
states, compe nsating the smaller p — d hybridiza tion, the total moment retaining its integer 
value of 3/iB (jGalanakis and Mavropoulod . l2003l ) . 
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The essential difference of the many-body electronic structure for the lattice constants of 
GaAs and InAs is completely due to the difference in the position of the Fermi energy with 
respect to the minority-spin band gap, whereas the self-energy characterizing the correlation 
effects is not changed too much (Fig. [HT]) . The total density of states A^ {E) is rather sensitive 
to the difference between the band edge E^ and the Fermi energy Ep. If this difference is 
very small (i.e., the system is close to the electronic topological transition Ec -^ Ep) one can 
use a simple expression for the singular contribution to the bare density of states, 5Nq {E) oc 
^/W^^Wc {E > Ec). The appearance of the complex self-energy S (E) = Si (E) — 2^2 (E) 
changes the singular contribution as 



5N (E) oc [JZf (E) + Si (E) + Z, (E)] 



1/2 



(197) 



where Zi (E) = E — Ec — T^i {E) (cf. Ref. iKatsnelson and Trefilovl . 119901 ). Assuming that 



the self-energy is small in comparison with E — Ec one can find for the states in the gap 



6N (E) oc S2 (E) I \JEc — E {E < Ec) . One can see that the shift of the gap edge changes 
drastically the density of states for the same S2 (E) . 

A practical use of tunable properties of NQP states in CrAs grown on different substrates 
is possible. For most of the applications room-temperature and the stability of the ferromag- 
netic state are the most important prerequisites. The fer romagnetic CrAs m ight be grown 



on 111-V semiconductors similar to the zinc-blende CrSb 



of the NQP states was obtained in Ref. (jChioncel et al 



Zhao et al 



m 



2OOII ). The presence 



20051 ) for CrAs lattice parame- 



ters larger than 5.8A. It was found experimentally that at 300K, ar ound this value of the 
lattic e parameter, a stable solid solution of Gao.65Ino.35As is formed (JHarland and WooUeyl . 
19661 ). Thus, from the practical point of view, a 65% of gallium in a Ga^Ini_a;As compound 
would constitute the ideal substrate for the CrAs half-metal. This could be a part of the 
epitaxial III-V structure providing an easy way to integrate with the existing semiconductor 
technology. 

To determine possible substrates for growth of layered half-metallic materials, electronic 
structure calculat ions were carried out for lattice parameters in the range 5.60-i- 6.03A 
( jFong et ali I2004J ). According to these calculations, growth with minimal strain might be 
accomplished in a half-metallic multilayer system grown on InAs substrate which would be 
the best choice to evidence the NQP states, since the Fermi energy is situated in this case 
far enough from the bottom of the conduction band. 
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A very high sensitivity of the minority-electron DOS near Ep to the lattice constant 
opens a new interesting opportunity. Suppose we have an antiparallel orientation of the 
magnetizations in the CrAs-based tunnel junction (such as shown in Fig. 12^). then the 
I—V characteristic is determined by the density of the NQP states. Thus if we will influence 
the lattice constant (e.g., using a piezoelectric material) one can modify the differential 
conductivity. This makes CrAs a very promising material with tunable characteristics which 
opens new ways for applications in spintronics. 



2. VAs: correlation-induced half-metallic ferromagnetism? 



Intere s ting materia . 



( JNagaev 



1983 



Ohno 



s for sp intronics applications are ferromagnetic semiconductors 



1998a 



b|). Candidate systems are ordered compou nds such as eu 



ropium chalcogenides (e.g., EuO) and chromium spinels (e.g., C dCr^Se ^ ) (INagaevl . Il983l ). 



as well as diluted magnetic semiconductors (e.g., Gai_2,.Mn2^As) (jOhno 



1998al ]bh. Unfor- 



tunately, all of them have Curie temperatures much lower than room temperature. On 



the other 



land, VAs in the zinc-blende stru cture is, according to density-functional cal- 



culations (JGalanakis and Mavropoulos 



Curie temperature. Un l ike C rAs (jAkinaga et al. 



20031) . a fe r romag netic se miconductor with a high 



2000|), CrSb f Zhao et al. 



200l|), and 



MnAs (jOkabayashi et a/.l . |2004| ). VAs has not yet been experimentally fabricated in the zinc- 
blende structure, but the increasing experimental activity in the field of the (structurally 
metastable) zinc-blende ferromagnetic compounds is promising in this respect. 

The main result including many-body correlation effects is displayed in Fig. [331 While 
this material is expected to be a ferromagnetic semiconductor from density-functional theory 
(LSDA/GGA) or static LSDA+f/ calculations, the inclusion of dynamic Coulomb correla- 
tions within the LSDA+DMFT approach predicts a majority-spin band metallic behavior 
due to the closure of the gap (~50 meV). However, since the minority-spin band gap (~ leV) 
remains finite, the material is found to be a half-metallic ferromagnet. To our knowledge, this 
i s a first examp l e wher e dynamic correlations transform a semiconductor into a half-metal 
( IChioncel et all l2006bl ). This remarkable result demonstrates the relevance of many-body 
effects for spintronic materials. 



The important features of the electronic structure of VAs (JGalanakis and Mavropoulod . 



20031 ) are shown schematically in Fig. [321 The t2g states hybridize with the neighboring As 
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p states, forming wide bonding and antibonding hybrid bands. In contrast, the Cg states 
form mainly non-bonding narrow bands. The Fermi level lies between Cg and antibonding 
t2g states in the majority-spin bands, and between bonding t2g and Cg in the minority-spin 
bands. The spin moment, concentrated mainly at the V atoms, is an integer of exactly 
M = 2 fiB per unit formula, which is obvious from counting the occupied bands for two spin 
directions. 

The exchange constants of VAs were calculated withi n GGA and adiabati c spin dynam- 
ics approach, similar to the one used by Halilov et al. (JHalilov et ali Il998l ). Using these 
exchange parameters in a Monte Carlo simulation of the corresponding class ical Heisenberg 
Hamiltonian yields a Curie temperature Tc = 820 K (Chioncel et all l2006bl ) by the fourth- 
order cu mulant crossing poi nt. This result agrees with the value of Tc = 830 K calculated 
in Ref. (jSanyal et all |2003| ) by using a similar method. The high Curie point is well above 
the room temperature, which makes VAs a very promising candidate for applications in 
spintronics. 

Static correlations were taken into account within the LDA+U method using similar 
values of effective interactions parameters U = 2 eV and J = 0.9 eV as in the case of 
Heusler alloys. It will give just an estimation of correlations effects in VAs, since for exact 
value of U one need t o per f orm a complicated a nalyse of screening effects in this compound 



( JAryasetiawan et al 



2004 



Kotliar et al 



20061 ). The theoretically determined equilibrium 



lattice parameter, a = 5.69A, and a broadening 5 of about 15K, which allows the majority 
spin gap to be clearly resolved, were used. For differ ent lattice parameters (e.g., as for 



InAs) the LDA re sults agree with the previous ones ( JGalanakis and Mavropoulos 



Sanyal et al. 



2003; 



2003h . The GGA DOS is shown in Fig.[3S The main difference in the GGA+U 



spectrum is that the occupied localized majority eg states are expected to shift to even 
lower energy, while the unoccupied minority eg states to higher energy. The semiconducting 
character does not change, since the eg and t2g bands remain separated for both spins; the 
majority-spin gap slightly increases, but remains small. 

In order to investigate dynamic correlation eff ects in VAs the fully self-consisten t in 



spin, charge and self-energy LSDA+DMFT scheme (jChioncel et al 



2003a . 



2005 



2003b|) was 



used. The computational results for the DFT in the GGA approximation and LDA+DMFT 
densities of states are presented in Fig. [331 The non-quasiparticle states in the minority spin 
band are visible just above the Fermi level (inset), predicted also by previous calculations 
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(IChioncel et all l2003al . l2005l ). The weak spectral weight of NQP states is due to that the 
Fermi level is close t o the right edge of the minority spin gap, as discussed for CrAs having 



a similar structure (IChioncel et al. 



20051 ). The local spin moments at V atoms do not 



change significantly (by less than 5%). However, in the case of VAs another correlation 



effect appears: t 



(IChioncel et al. 



le sma ll majority-spin gap at Ep closes, making the material half-metallic 



2006bh . 



In order to investigate the mechanism of the gap closi ng for the majo r ity spin channel, 
the behavior of the electron self-energy was investigated (IChioncel et all l2006bl ). For the 
majority-spin, ImL^^E) ~ {E — Ep)'^ which indicates a Fermi-liquid behavior, as opposed to 
ImS|(_E) which shows a suppression around Ep due to the band gap, as well as a peculiar 
behavior for E > Ep related to the existence of NQP states. 

From the Dyson equation (I166P one can see that the real part ReSo-(-E) causes a shift of 
the LDA energy levels. Therefore, due to the non-zero S!", the Cg orbitals in the close vicinity 
of the Fermi level are pushed closer to Ep. This renormalization is connected with the large 
absolute value of Re{dT,/dE)pp < 0. This causes occupied levels to be shifted to higher 
energy and unoccupied levels to lower energy. Note that this effect is completely opposite 
to the LDA+f/ results discussed above. In addition to this shift, the Cg peak is broadened 
by correlations, its tail reaching over the Fermi level (Fig. [331 inset). Thus the finite- 
temperature LDA+DMFT calculations demonstrate the closure of the narrow gap in the 
spin-up channel, which is produced by the correlation-induced Fermi-liquid renormalization 
and spectral broadening. At the same time, NQP states appear for the minority-spin channel 
just above Ep. 

The slope of the majority-spin self-energy is almost a constant as a function of tempera- 
ture at low T: Re{dT,^/dE)pp ~ —0.4 between 200 K and 500 K. The quasiparticle weight, 
which measures the overlap of the quasiparticle wave function with the original one-electron 
one for the same quantum numbers, is Z = {1 — dUeTj-^/dE)^^ ~ 0.7. As a consequence, 
the closure of the gap in the majority channel is a quantum effect originating from the 
multi-orbital nature of the local Coulomb interaction (energy states are shifted towards Ep) 
rather than an effect of temperature. A similar gap closure is obtained for larger values of 
U, namely, U = 4 and 6 eV, although the latter values should be taken with some caution 
in the FLEX calculation, which is in principle appropriate only in weak to intermediate 
coupling. As a general tendency, increasing U* = U — J produces a stronger Fermi-liquid 
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renorinalization in the majority spin channel, the same effect being evidenced for J = 0. 

On the one hand, density-functional theory calculations within the GGA (Ref. 
( IChioncel et all l2006bl )) predict this material to be a ferromagnetic semiconductor with 
a tiny gap of about 50 meV in the majority-spin states and large gap of the order of 1 
eV for minority-spin states. Quantum effects, such as spin and orbital fluctuations, de- 
scribed by LDA+DMFT destroy the narrow gap and turn the material into an half-metallic 
ferro magnet. 

On the other hand, several other mechanisms could contribute to the band-gap narrowing 
with increasing temperature. A well studied example is electron-phonon interaction. In 
another semiconductor with the zinc-blende structure, GaAs, it amounts to 50meV at 200K 



( jPaesslerl . Il999l ). Also spin-orbit coupling can be essential when considering closing of the 
gap. 



The LDA/GGA calculation su pplemented 



high Curie temperature of 830 K (jSanyal et al. 



j y a M onte Carlo simulation also predicts a 



20031 ) , which makes this material of interest 



for technological applications. One can expect that Tq is not strongly affected by dynamical 
correla tion, for the same reason as the effe ctive exchange interaction parameters (see recent 
works (JKatsnelson and Lichtensteinl . |2002| ) ) . 

The revealed half-metallic (instead of semiconducting) behavior has important conse- 
quences in the potential appli cations of VAs in s pintronics. In contrast to all semiconductor- 



based spin-injection devices ( IZutic et al\ . l2004l ) which avoid the resistivity mismatch prob 



lem, half-metals can b e applied to obtain g i ant rn agnetoresistance or, provided that interface 
states are eliminated ( jMavropoulos et all l2005l ). tunneling magnetoresistance effects. We 
see that correlation effects play a decisive role in the prediction of new spintronic materials. 
The metallic nature of the majority spin channel would be visible in resistivity measure- 
ments. Therefore, the experimental realization of zinc-blende VAs would provide a test of 
this prediction. Further research should address the issue of the stability of the half-metallic 
ferromagnet ic state in a zinc-blende structu re. Some work in this direction has been already 



carried out (jShirai 



2003 



Xie et al. 



2003bh . 
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C. Half-metallic transition metal oxides 



1. Cr02: a rutile structure half-metallic ferromagnet 



Chromium dioxide Cr02 has a rutile str ucture with a = 4 .421 A, c = 2.916 A (c/ 



0.65958) and internal parameter u = 0.3053. (jPorta et a/.l . Il972l ) The Cr atoms form a body- 



center tetragonal lattice and are surrounded by a slightly distorted octahedron of oxygen 
atoms. The space group of this compound is non-symmorphic {PA2/mn'm = Dl"^). The 
Cr ions are in the center of CrOg octahedra, so that the 3d orbitals are split into a t2g 
triplet and an excited Cg doublet. The Cg states with only two valence 3d electrons are 
irrelevant, and only the t2g orbitals are to be considered. The tetragonal symmetry distorts 
the octahedra, which lifts the de generacy of th e t2r, orbitals into a d^y g round state and 



dyz+zx and dy-^-zx excited states (iKorotin et ah . 



1998 



Lewis et al. 



19971 ) (see Fig. [Ml a 



local coordinate system is used for ever y octahedron ) . A d ouble exchange mechanism for 
two electrons per Cr site was proposed (jSchlottmannl . l2003l ). According to this, the strong 
Hund's rule together with the distortion of CrOg octahedra leads to localization of one 
electron in the d^y orbital, while the electrons in the dy^ and d^z are itinerant. 

Measurements of the magnetic susceptibility in the param agnetic phase s how a Curie- 



Weiss-like behavior indicating the presence of local moments ( IChamberland 



19771 ). which 



suggests a mechanism of ferromagnetism beyond the standard b and or Stoner- 



Several recent experimental investigat ions of photoe mission (Tsuiioka et 



tics (ISingley et al. 



x-ra y absorption (Stagarescu et a/.l . l2000l ). resistivity (ISuzuki and Tedrow 



al. 



ike model. 



19971 ). soft 



19981 ). and op- 



19991 ) suggest th at electro n corre lations are essential for the underlying 



physical picture in Cr02. Schwarz ( ISchwarzl . Il986l ) first predicted the half-metallic band 
structure with a spin moment of 2^b per formula unit for Cr02. Lewis fJLewis et all 119971 ) 
used the plane-wave potential method and investigated the energy bands and the transport 
properties, characterizing Cr02 as a "bad metal" (a terminology applied earlier to high 
temperature superconductors and to other transitio n metal oxides, even ferromagnets like 
SrRuOs). A decade later the LSDA+f/ calculation (IKorotin et all Il998l ) treated the con- 
ductivity in the presence of large on-site Coulomb interactions and described Cr02 as a 
negative charge-transfer gap material with a self-doping. Contrary to the on-site strong cor- 
relation description, transport and optical properties obtained within the FLAPW method 
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(LSDA and GGA) (JMazin et all Il999l ) suggest that the electron-magnon scattering is re- 
sponsible for the renormahzation of t he one-electron bands. More recent model calculations 



( ICraco et al. 



2003 



Laad et al. 



2OOII ) propose even orbital correlations. 



Chemical bonding in rutile-type c ompo unds including Cr02 was analyzed by Sorantin 



and Schwarz (ISorantin and Schwarz 



1992). One can see that around the Fermi level the 



bands are primarily chromium 3(i states of t2g manifold, with Cg bands being situated higher 
in energy due to the crystal-field splitting. In the spin polarized case, the exchange splitting 
shifts the minority spin d bands above the Fermi level (Fig. [35]) . For the majority ^23 bands 
the Fermi level lies in a pseudogap. Oxygen p — chromium d hybridization creates both 
bonding and antibonding hybrid orbitals, with the bonding orbital appearing in the occupied 
part and the antibonding hybrid orbital remaining in the Cr ^23 manifold. Half of the dyz and 
dzx components of t2g are pushed upward by antibonding, which explains the dxy dominance 
in the spin density. The non-magnetic DOS shows a sharp peak at the Fermi level, which 
signals the magnetic instability according to the usual Stoner argument. 



Al though there i s a si g nificant difference between the t^g and Cg orb itals ( JKorotin et al. . 



1998 



Mazin et al. 



1999 



Schwarj . 



framework of NMTO technique (jAndersen and Saha-Dasguptal . 



19861: 



Sorantin and Schwarzl. 



9921) . the analysis in the 



2000 



Zurek et al. 



20051) 



shows that their interplay is important not only for the crystal-field splitting of t2g states, 
but also for t he general bonding in the rutile structure. The t2g orbitals form the basis set 
used in Ref. ( lYamasaki et alV l2006l ) to evaluate the effective hopping Hamiltonian matrix 
elements. 

In Cr02 the bands around the Fermi level are primarily chromium "id states of t2g mani- 
fold, eg bands being situated higher in energy owing to crystal-field splitting. The ^29 orbitals 
are further split into single d^y and nearly degenerate dyz±zx bands owing to the orthorhom- 
bical distortion of CrOg octahedra. Despite the differences between Cr t2g — eg orbitals, their 
interaction plays an important role not only in characterizing the crystal-field splitting, but 
also in the general picture of bonding in the rutile structure. 

Concerning the Coulomb interaction U in Cr02, the higher energy eg bands, although 
making no n oticeable cont r ibution at the Fermi level, c ould participate in the screening of the 
t2g orbitals (jPickett et ali 



1998 



Solovyev et al. 



19961 ) ■ thereby giving the values f/ = 3 eV 



and J = 0.87 eV. 



Previous LSDA-|-t/ ( JKorotin et al 



1998 : 



Toropova et ffl/.l . l2005l ) and DMFT fICraco et al. 
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20031 : iLaad et all I2OOII ) studies yielded independently a narrow almost flat band of dxy 
character, which produces ferr omagnetism in Cr02- In contrast to these results, the fully 



self-consistent LSDA+DMFT 



variational cluster approach (jChioncel et al. 



Chioncel et al. 



20071) yields, in agreement with a non-local 



20071 ). an itinerant d^y orbital which crosses 



the Fermi level. Despite the non-localized character of the orbital, a ferromagnetic phase is 
still obtained. 

Results of the LSDA-fDMFT calculation are presented in Fig. [36] for two different values 
of T and compared with the LSDA results. The LSDA Fermi level intersects the majority- 
spin bands near a local minimum and lies in the band gap of the minority spin. Finite 
temperatures and correlation effects close this pseudogap around the Fermi level, as can be 
seen from the LDA+DMFT results in FiglSBl No differences can be observed between the 
two DMFT results at different temperatures, except for the smearing of DOS features at 
larger temperature. For both spin channels the DOS is shifted uniformly to lower energies 
in the energy range —2 -. — 6 eV, where predominantly the 0{p) bands are situated. This is 
due to the Cr(d) bands which affect the 0(p) states through the Ci{d) — 0{p) hybridization, 
so that 0{p) states contribute actively to the ferromagnetic ground sta te formation. 



Fig. 



presents the experimentally measured spin polarization (JHuang et al 



20031) 



in comparison with the theore tical calcu 



ture LDA+DM 



( iKurmaev et al. 



T (T =200 K) (jChioncel 



ations within the LSDA and finite tempera- 



20061 ) . Resonant X-ray emission spectroscopy 



2OO3I ) showed the existence of nearly currentless minority spin states in 



the vicinity of the Fermi level, which can be connected to the non-quasiparticle states Illl.Ji 
As described in previous sections IIV.AI and |V] electronic structure of several half-metallic 
ferromagnets reveal the existence of such NQP states, which are important for the proper 
description of spin polarization near the Fermi level. The LSDA+DMFT description, how- 
ever, is not sufficient to capture the high energy tail of the experimental spin polarization. 
This can be related with improper description of unoccupied Cr Cg orbitals in the LSDA, 
and probably, al so with non-local ex change effects which can be investigated within a cluster 
DMFT scheme ( IKothar et alimm . 



The NQP states occur aro und 0.25 eV in t 



the variational cluster results (IChioncel et al. 



le sp in- down channel. At zero temperatures. 



20071 ) yield the states in the minority channel 



which are actually far from E p, but have a tail vanishing a t Ef. This is in agreement with 



low temperature experiments (jJi et al. 



2001 



Soulen et all Il998l ) which support very high 
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polarization of Cr02. 



2003) 



The effects of local and non-local electronic correlations in Cr02 (iChioncel et al 
change considerably the mean-field LSDA+f/ picture, despite the interaction is not too 
strong. In particular, in LSDA+f/ the single occupancy of the Cr d^y orbital is determined 
by the exchange and crystal-field splitting. On the other hand, the competition of the latter 
with correlation effects, which is taken into account in the DMFT and variational cluster 
perturbation theory calculations, induces a ferromagnetic state with itinerant-type dxy or- 
bitals possessing a large effective mass rather than with localized rnoments, in contrast to 



previ ous results ( jCraco et al. 



2003 



Korotin et al. 



19981 : 



Laad et al. 



2001 



Toropova et al. 



20051 ) . In the minority-spin channel, correlations induce NQP states which are crucial for 



the occurrence of spin depolarization in Cr02. However, a quantitative analysis of the de- 
polarization would require the inclusion of additional effects, e.g., disorder or phonons (sect. 
UlLCj) . 

VI. EXCHANGE INTERACTIONS AND CRITICAL TEMPERATURES IN HALF- 
METALLIC COMPOUNDS 



Owing to the strong interest in the half-metallic ferromagnetism, the number of theoreti- 
cal studies of exchange interactions and calculations of Curie tem peratures in Heusle r alloys 



has been drastically increased last time (see, e.g., recent review ( jZhang et al. 



2003)). The 



first investigation s of the exchange in teractions in half-metals within the DFT formalism was 
made by Kiibler (JKiibler et a/.l . Il983l ). The mechanisms of ferromagnetism in Heusler alloys 
were discussed on the basis of total-energy calculations for the ferro- and antiferromagnetic 
configurations. Since the antiferromagnetic state of the system is not half-metallic, such 
an estimation gives only crude values of exchange parameters in HMF. Therefore, a more 
precise evaluation of the exchange interactions from the first-principle theory is required. 
In this section we present the real-space Green's function and frozen-magnon techniques to 
calculate the exchange parameters, and their applications to fuU-Heusler, semi-Heusler, and 
zinc-blende half-metals. 
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A. The Green's function formalism 

Within the first-principle Green's function approach the exchange parameters are ob- 
tained by mapping the second variation of electronic band energy to the classical Heisen- 
berg Hamiltonian and making use of magnetic version of the Andersen's local force theorem 



(JMachintosh and Andersenl . Il980l ). 



^e// = -^EE«'«R«R' (198) 

fJ-,!^ RR' 

In Eq. (11981) . the indices fi and u mark different sublattices, R and R' are the lattice 
vectors specifying the atoms within a sublattice, s^ is the unit vector in the direction of the 
magnetic moment. By introducing the generalized notation for site (i = /i, R), a simple and 
trans parent expression for the e xchange interaction parameters is obtained in the following 



form (JLiechtenstein et al\ . 119871 ) : 



Jij = ^ / deImTrL{A,Gl^AjGlj}, (199) 

where G^- is the real space Green's function and Aj is the local exchange splitting for site i. 

B. The frozen-magnon approach and DFT calculations of spin spirals 

The approach is based on the evaluation of the energy of the frozen-magnon configurations 
defined by the following atomic polar and azimuthal angles: 

e^ = e;<P^ = ci-R + cP>^ (200) 

The angle 6 defines the cone of the spin spiral, and the constant phase 0'^ is normally chosen 
to be zero. The magnetic moments of other sublattices are kept parallel to the z-axis. Within 
the classical Heisenberg model (11981) the energy of the spin-spiral configuration is 

E'^'ie, q) = E^^{9) + J'^^(q) ■ sin^O (201) 

where Eq^{9) does not depend on q and the Fourier transform J'^^(q) is defined by: 

J'^^(q) = E^Rexp(*q-R) (202) 

R 

In the case where u = fi the sum in Eq. (I202p does not include the local interaction with 
R = 0. Calculating E^^{6,q) for a regular g-mesh in the Brillouin zone of the crystal 
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and performing the inverse Fourier transformation one gets exchange parameters Jg^ for 
sublattice fi. 

The Curie temperature can be estimated within the mean- field (MFA), 



kBT^'' 



E^^R 



M 1 



R 



6/iB A^^ 



E^(q) 



(203) 



and random phase approximation (RPA) 

1 6/iB 1 



E 



(204) 



with uj{q) the spin-wave dispersion, N the number of q points in the first Brillouin zone, 
and M the atomic magnetic moment. In the MFA the Curie temperature is determined by 
the arithmetic average of the magnon energies, while in the RPA Tc is determined by the 
harmonic average. Therefore the value of Tq within MFA is larger than the RPA one, the 
two values being equal only provided that the magnon spectrum is dispersionless. 

C. First-principle calculations 

1. Semi-Heusler Clj, alloys 

First-principle studies of exchange interactions and magnetic phase tra nsitions for NiCrZ 



Z=P,Se,Te), NiVAs, NiMnSb, CoMnSb were carried out by many authors (JRusz et al 



Sasioglu et al 



2006 



2005b[) In fig. [38] the f rozen- magnon dispersion and the exchange interactions 



are presented (ISasioglu et all i2005b) . The exchange interactions are Fourier transforms of 



the frozen-magnon spectrum. 

A remarkable feature of the exchange interactions is their short-range character, the Curie 
temperature being determined by the interaction within first two coordination spheres. The 
deviation from th e coUinear a li gnment of magnetic moments can b e also characterized by 



magnon energies (jSandratskii 



1998 



Sandratskii and Bruno 



200 



The deviation leads 



to the mixing of the majority and minority spin states, which makes half-metallicity less 
favorable. A detailed discussion of different depolarization mechanisms is given in section 
HTLUl 
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2. FuU-Heusler L2i alloys 



Recently the studies of the interatomic ex change interact i ons in several fuU-Heusler com 



pounds were reported by Kurtulus et al. (JKurtulus et al 



20051 ) Sasioglu et al. studied 



the exchange interactions in non-half-metallic Ni2MnZ (Z=Ga, In, Sn, Sb) and half-metallic 
Mn2VZ (Z=Al,Ge). The importance of the intersublattice exchange interactions has been 
demonstrated. For Mn2VZ (Z=Al,Ge) it was shown that the ferrimagnetic coupling between 
the V and Mn moments stabilizes the ferromagnetic alignment of the Mn moments. 

In Co^^iMn Z (Z=Ga, Si, Ge, Sn) the presence of Co atoms makes the interaction more 
complicated (JKurtulus et all |2005| ). The interaction between the Co atoms in the same 
sublattice ,Col(2)-Col(2), and between Co atoms at different sublattices, Col(2)-Co2(l) has 
to be take n into account . This approach gives results which goes beyond the initial approach 



of Kiibler (IKiibler et al. 



19831 ) . From the Jij values it is clear that the exchange interactions 



are relatively short ranged and do not exceed the first four neighbors in each sublattice. The 
main exchange parameter corresponds to the nearest neighbor Co(l)-Mn interaction. This 
already gives 70% of the total c ontribution to J and is about ten times larger than the 



Co-Co and Mn-Mn interactions ( Kurtulus et al. 



20051 ). Thus, it was concluded that the 



Co-Mn interactions are responsible for the stability of ferromagnetis m fsee Fig. [39 



The interaction changes weakly with the contraction of the lattice (jSasioglu et al 



2005ch 



interactions between Co at different sublattices favors ferromagnetism and is stronger than 
the ferromagnetic interaction between the Mn atoms. The intrasublattice Mn-Mn ad Co- 
Co interactions have different signs for different distances between atoms, which means a 
RKKY-type interaction. 

3. Zinc-blende half-metals 



Besides the large number of results on the electronic properties of zinc-blende half-metals, 
the exchange interactions constitute an i mportant aspe ct to understand the stability of 



half-metallicity in these structures. Shirai (jShirai 



20031 ) obtained that in VAs, CrAs and 



MnAs ferromagnetism is energetically favorable in comparison with the a ntiferrom a gnetic 
state, unlike FeAs where an opposite effect was demonstrated. Sakuma (jSakumal . |2002| ) 
predicted ferromagnetism in the isoelectronic MnSi, MnGe and MnSn compounds. Similarly 
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to CrAs, the ferroinagnetism in these systems is stabihzed by short-range interactio ns (direct 
Mn-M n and indirect through sp atoms), giving a Curie temperature of lO OOK ( 



Sakumal . 



20021). Tc for the VAs, CrAs and MnAs was calculated also by Kiibler flKiiblerl . l2003h 
yielding the same ra nge of magnitudes. Using GaAs and InAs lattice consta nts, Sanyal et al 
(ISanval et ali l2003l ) calculated Tc for VAs, CrAs and MnAs. Sasioglu et al (ISasioglu et al. . 
2005al ) calculated the exchange parameters for a large number of pnictides. 

MnC presents an interesting situation since its half-metallic gap is situated in majority 
spin channel. The Curie temperature was found to be 500K. FigJlQ] represents the frozen- 
magnon energies for a selected direction in the Brillouin zone and the exchange constants. 
Ferromagnetism is stabilized by the direct Mn-Mn interactions and Mn-C ferromagnetic 
coupling. A remarkable feature o f MnC is small difference between the MFA and RPA 



values of Tr-. (ISasioglu et al. 



2005ah . 



VII. CONCLUSIONS 



The idea of half- metal 



culations (Ide Groot et al. 



ic ferr omagnetism appeared as a result of band-structure cal- 



1983bl ). During a long time, dominant activity in this field was 



theoretical. Conceptually, HMF are of interest from the point of vie w they provide an oppor- 
tunit y to probe in a clear form some essentially many-particle effects (jlrkhin and Katsnelson . 
1990[) . Whereas for generic me t allic system the Landau Fermi-liquid theory (jNoziered . 
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I 



Vonsovsky and Katsnelson 



19891 ) works, most correlation effects being hidden in a 



parameter renormalization (such as effective mass, magnetic moment etc.), in HMF the 
spin-polaronic effects lead to a qualitatively new feature: occurrence of non-quasiparticle 
(incoherent) states in the energy gap for one of spin projections near the Fermi level. On 
the contrary, similar effects of electron-magnon interactions in traditional itinerant-electron 
ferromagnets are mixed with other kind renormalizations (e.g., electron-phonon, electron- 
electron ones etc.). The NQP states occur only above (below) Ep for minority spin (majority 
spin) ga p. Therefore HMF are i deal o bjects to investigate effects of electron-magnon inter- 
actions (jlrkhin and Katsnelsonl . Il990l ). 

Even stronger motivation to study HMF is connected wit h perspectives to use 
them in giant magnetoresistance and t unnel magnetoresistance ( Ide Groot et al. . 



Irkhin and Katsnelson 



1994 



Prinz 



1983al : 



devices. This initiated great theoretical activity in 
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the fi eld of heterostructures containing HMF (jlrkhin and Katsnelsonl . l2002l : iTkachov et all 
20011 ). At the same time, the interest in the search and prediction of new HMF was growing 
on the basis of band-structure calculations, as well as attempts to understand better the 
features of electronic structure and chemical bonding which are relevant for half-metallicity 
(see SectJV]). 

Recently, numerous at tempts have been performed to build hetero s tructures wit h HMF, 



such as Heusler a lloys (IGercsi et al 



( Miao et al. 



2006 



Sakuraba et al. 



20061 ). Fe304 (IRvbchenko et al. 



2006 



2006a 



Zhao et al. 



2007 



2006bh . CrOs 



20061 ). Therefore the half- 



metallicity predicted by electronic structure calculations becomes practically applicable. 
Nevertheless, direct experimental evidences of half-metallic structure for specific compounds 
are still rather poor. Perhaps, the unique method of testing genuine, bulk, half-metallic 
properties remains the spin-resolved positron annihilation. This underexposed technique 
enables the direct measurement of the spin-polarization in the bulk. Advanced techniques 
borrowed from semiconductor technologies which access spatially resolved spin polarization 
at the Fermi level would be interesting alternatives for positron annihilation. Although being 
extensively used to characterize semiconductors, they are poorly known for the spintronics 
community. 

Among other expe r iment s expected to advance the fi eld we mention the STM 



(jlrkhin and Katsnelson 



dreev refiection (jSoulen et al. 



2006) 



spin- polarized photoemission (jPark et all Il998l ). and An 



19991 ). Investigations of nuclear magnetic relaxation rate 



should be specially m entioned since the absence o: 



of the half-metallicity ( jlrkhin and Katsnelsonl . l200ll ) and NMR gives is not surface sensitive 



the Korringa relaxation is a clear sign 



and gives direct information on bulk properties. Also data of core-level spectroscopies, espe- 
cially XMCD would be very useful (see Sect JllLJ]) . As for the band-structure calculations, 
application of state-of-art methods taking into account correlation effects, such as GW or 
DMFT (see Sect ilV.Al) looks very promising. Another perspective direction is the use of 
electron-structure calculations to search essentially new types o f HMF like sp electron (or 



anionic) magnets (jAttema et al. 



2005 



Edwards and Katsnelsonl . 



20061 ). 
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Figures 




FIG. 1 (color online) Density of states of non-quasiparticles, for half-metallic ferromagnets, pos- 
sessing the gap in the minority spin channel. NQP states are the dominant many-body feature 
around Ep in comparison other mean- field effects, such as spin-orbit or non-collinearity as it will 
be discussed in the following sections. 
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FIG. 2 Atom-resolved DOS for the Co2MnZ compounds with Z=A1, Si, Ge, Sn ([Galanakis et all 
20M). 
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FIG. 3 The density of states (DOS) of Sr2FeMo06 (IKobavashi e.t d\ . Il998l ) 
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FIG. 4 Density of states for the semielliptic DOS at concentration of carriers current 6 = 0.35 
{6c < 6 < 6'^) and 6 = 0.02 < 6^ (inset) in the self-consistent approximation. 
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conduction band is shown by dashed line; the short-dashed line corresponds to the mean-field 
approximation . 
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FIG. 7 Density of states in a half-metallic ferromagnet with / = —0.3 < 0, other parameters being 
the same as in Fig.l. The spin-down subband (lower half of the figure) nearly coincides with the 
bare band shifted by IS. Non-quasiparticle states in the spin-up subbands (upper half of the figure) 
occur below the Fermi level; the short-dashed line corresponds to the mean-field approximation. 
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FIG. 11 The local density of states N^^JyE) (solid line) for a half-metallic ferromagnet with S = 
1/2, / = 0.3, 5 = 0.01 in the presence of the core hole potential U = —0.2. The dashed line shows 
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0.3- 


N,(E) 


^ 




gO.2- 




/ 








1 


J 






1 1 


// 






/ 1 


7 


0.1- 




/ 1 

/ 1 

/ 1 

/ / 


I 


0.0- 


1 1 


' 1 ' 


1 1 1 ' 



-0.3 



-0.2 



-0.1 



0.0 
Energy 



0.1 



0.2 



FIG. 12 The local density of states N^^JJii) (solid line) for a half-metallic ferromagnet with S = 
1/2, / = —0.3, 5 = 0.01 in the presence of the core hole potential U = —0.3. The dashed line shows 
the DOS N-^{E) for the ideal crystal. The value of Ep calculated from the band bottom is 0.15. 
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FIG. 13 Density of states for HMF in the Hartree-Fock (HF) approximation (dashed hne) and 
in QMC solution of DMFT problem for semi-circular model (solid line) with the band-width 
W=2 eV, Coulomb interaction U=2 eV, A=0.5 eV, chemical potential ^=-1.5 eV and tempera- 
ture r=0.25 eV. Ins ets: imaginary part of t he the spin-flip susceptibility (left) and imaginary part 
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FIG. 14 (Color online) Cl^ structure with the fee Bravais lattice (space group F43m). Mn (green) 
and Sb (purple) atoms are located at (0, 0, 0) and (^, ^, ^) forming the rock salt structure 
arrangement. Ni (orange) atom is located in the octahedrally coordinated pocket, at one of the 
cube center posit i ons (4 , f > f ) leaving the other ( j, j, j) empty. This creates voids in the structure. 
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FIG. 15 (Color online) NMTO Mn-d Wannier orbitals of NiMnSb. Ni is orange, Mn is green, and 
Sb is purple. Red/blue indicates a positive/negative sign. Upper panel: t2g orbitals; dxy (left), 
dyz (middle), dzx (right). The triply-degenerate t2g orbitals can be obtained by the permutation 
of axes. Lo wer panel: Cp orbitals: d^^^^i (left), d^-^-y^ (middle). These eg orbitals are doubly 



degenerated. (JYamasaki et al. 
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FIG. 16 Pull basis set spin-polarized (ferromagnetic) bands for NiMnSb; majority spin (left) 
and minority spin (right). The high-symmetry points are; W{^,1,0), L(^,^,^), r(0,0, 0), 
X(0 1,0), K(^^ 0) in W-L-r-X-W-K-r line, and X(0,0, 1), UU, i 1), L(i i i), in X-U-L 



line . ( Yamasaki et al. 



20061) 
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FIG. 17 Decorated (fat) bands for the spin-polarized (ferromagnetic) NiMnSb; Mn majority spin 
(up) and minority spin (down). The high-symmetry points are; W{^,1,0), L(i, i, i), r(0,0, 0), 
X(0, 1,0), K(\l,0) in W~L-r~X-W-K-r line, and X(0,0, 1), [/(i, i, 1), L(i, i, i), in X-f/-L 



line ( Yamasaki 



20061 ). 
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FIG. 18 Decorated (fat) bands for the sp i n-pol arized (ferromagnetic) NiMnSb; Ni majority spin 



(up) and minority spin (down) (JYamasakil . 
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FIG. 19 Upper figure: FLAPW calculation of spin-polarized NiMnSb (iLezaid . 




20061 ) ■ Lower figures: 



Decorated (fat) bands for the sp i n-pol arized (ferromagnetic) NiMnSb; Sb majority spin (up) and 



minority spin (down) ([Yamasaki 



2001 
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FIG. 20 Density of states for HMF NiMnSb in LSDA scheme (dashed hne) and in LDA+DMFT 
scheme (sohd hne) with effective Coulomb interaction U=3 eV, exchange parameter J=0.9 eV and 
temperature T=300 K. The non-quasiparticle state is evidenced just above the Fermi level. The 
imaginary part of self-energies ImSj^ for t2g (solid line) and Cg (dotted line), ImS^ for t2g (dashed 
line) and eg (dashed dotted line) respectively. 
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FIG. 21 Spectral weight of the non-quasiparticl e state, calculated as a function of average on-site 



Coulomb repulsion U at temperature T=300 K. ljChioncel et al. 



2003a|) 



144 



CO 

o 

Q 




-4 -2 

E-Ep (eV) 




E-Ep {eV) 



FIG. 22 Density of states for the half-metallic ferromagnet NigMnyHoSbg in the case of anti-parallel 
3d — 4/ coupling left, and parallel coupling right. 
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FIG. 23 Density of states of semiconducting HoNiSb . (|Chioncell . 



200J) 
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FIG. 24 (Color online) Left: Density of states of half-metallic FeMnSb, LSDA (black line) and 
LSDA+DMFT (red line), for the efective Coulomb interaction U=2 eV exchange parameter 
J=0.9 eV and temperature r=300 K. Lower right panel: zoom around Ep for different values 
of U. Upper right panel: Sp ectral weight of the NQ P states calculated as a fun ction of U. The 



values obtained for NiMnSb (IChioncel et al 
2006al ) 



2003al ) are plotted for comparison. ([Chioncel et al 
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FIG. 25 (Color online) Contour plots of polarization as a function of energy and temperature for 
different values of local Coulomb interaction U. Left U=2 eV, right [7=4 eV. The LSDA polarization 
is plotted as the r=0 K temperatur e result. The asymmetry of the NQP states, is clearly visible 



for [7=4 eV. (Chioncel et al. 
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FIG. 26 (Color online) Temperature dependent polarization at the Fermi level, P{E = Ep,T) 
(so lid line) and magneti zation (dashed line) for different values of local Coulomb interaction 
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■ ^^IG. 27 Schematic illu stration of the gap formation in Co2MnZ compounds with Z=A1, Si, Ge,Sn 



(Galanakis et al. 
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FIG. 28 Total density of states of Co2MnSi full-Heusler alloy. Notice that the gap is formed be- 
tween the occupied Co-Co antibonding t2g orbitals and empty Co-Co antibonding eg orbitals. The 
LDA+DMFT results are presented as well for U = 3eV, and J = 0.9eV, and different temperatures. 
The inset shows the finite temperature spin-polarization. 



150 



(a) 




(b) 

tt 



ft tt 



tc5^t tt<5^ 



tt 

(c) 

lent 



tt tt 



tt 



a 



tt 



u 
It 
u 



tt 

'A 

tt 



t tit t^ 
ttjtcf 
ttjtt 

uitt 

Itltd" 

Uitt 



FIG. 29 The tunneling transport between strongly correlated ferromagnets. The density of states 
in the lower Hubbard band (a) is provided by standard current states for majority-spin electrons 
(above) and by non-quasiparticle states for minority-spin electrons (below), the latter contribution 



being non-zero only above the Fermi energy (occupied states are shadowed) . E 



is possible both for parallel (b) and antiparallel (c) magnetization directions. ([Chioncel et al. 
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FIG. 30 (Color online) Cr density of states calculated in LSDA (dashed line) and LSDA+DMFT 
(solid line) methods corresponding to a temperature of T=200K, average Coulomb interaction 
parameter U=2 eV and exchange J=0.9 eV. The non-quasiparticle states are clearly visible for 
lattice parameters larger than a^q = 5 .8^, in the unoccupie d part for minority spin channel just 



above the Fermi level, around 0.5 eV. (jChioncel e t alf, 2005) 
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FIG. 31 (Color online) Energy dependences of imaginary parts of the electron self-energy E2(i?), 
for lattice constants of GaAs (a), equilibrium one (b), and of InAs (c): e^ down solid line, t^.n down 



decorated solid line, Cg up dashed line, t2g up decorated dashed line. (JChioncel et al. 
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FIG. 32 Schematic r epresen tation of the p-d hybridization and bonding-antibonding sphtting in 



VAs. dChioncel et al. 
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FIG. 33 (color online) DOS of VAs within the GGA (dashed blue hne) and GGA+DMFT (solid 
red line) for a temperature of T=200 K, U=2 eV and J=0.9 eV. Inset: Focus around Ep showing 



the semiconducting gap within the GGA. To i 



larger scale for the spin down channel is used. (|Chioncel et al. 



lustrate the rn i nority- spin NQP states a ten times 
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FIG. 34 (Color online) Cr02 (rutile) structure. Crl (green) and Cr2 (orange) are located at (0, 
0, 0) and (2,2)2)' ^^ atoms are octahedrally coordinated by oxygen atoms (purple). The local 



coordinate system is used for each Cr atom; x 



-ia+ib 



1 



V2 



c, yi 



a + ib + ^c, 



1 



^1 = 75^+ 71^' ^"^^2 



1. 



ib - -i=c, y2 



ib + ^c, Z2 = -^a + ^b. xi,2 and 



2 2" v^^' J -^ 2" 2' 

yi,2 are approximately point to O atom, and zi.2 are exactly point to O atom. The local axes a re 



transformed into each other by a 90° rotation around the crystal c axis. (JYamasaki et al. 



200a) 
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FIG. 35 Full basis set spin-polarized (ferromagnetic) bands for Cr02; majority spin (left) and 
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FIG. 36 (Color online) Density of states obtained within the LSDA and LSDA+DMFT calculations 
for different temperatures. The in set shows the results for a smaller energy window around the 



Fermi level. ( Chioncel et al 



2003) 
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FIG. 37 (Color online) Energy dependence of spin-polarization of Cr02. Experimental measure- 



ment blue scattered line (JHuang et al. 



20031), LSDA calculation ( 



jlack dot-dashed), DM FT (green 



dashed) and non-local variational approach, red continuous line (jChioncel et al. 



2003) 
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FIG. 38 Left pannel(a): Frozen- magnon dispersion for NiCrZ (Z=P,Se,Te). Right pannel(b): 



Interatomic exchange interactions. (jSasioglu et al. 



2005b|). 
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FIG. 39 Exchange constants for Co2MnSi as a function of interatomic distance. (ISasioglu et al. 
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FIG. 40 Upper pannel: Frozen-magnon energies for MnC and MnSi in the [00*7] and [Oqq] directions. 
Vertical broken line shows the B rillouin zone boundary . Lower panel: Exchange constants as a 



function of interatomic distance. ([Sasioglu et al 
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(i^b) 


(K) (eV) (eV) 


,,GaAs 


3.191 


-0.270 


-0.009 0.089 


3.00 
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,,GaAs 
PDMFT 


3.224 
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-0.023 0.067 


3.00 
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A'lda 


3.284 
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/^DMFT 
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TABLE I Summary of the results of the calculations (iChioncel et all 

moments corresponding to the GaAs, InAs and the equilibrium lattice constant Ueq- For the latter 

one the value aeq = 5.8 A wa s used. Parameters of the DMFT calculations are presented in the last 



three columns of the table. (Chioncel et al 



2005|) 
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